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1. Introduction. The equation under discussion in this paper is the following: 


Ou Ou au 
— t—— =P 5 
Ot Ox ox’ 


(1) 


where u = u(x, t) in some domain and » is a parameter. The occurrence of the first 
derivative in ¢ and the second in zx clearly indicates the equation is parabolic, similar 
to the heat equation, while the interesting additional feature is the occurrence of the 
non-linear term u du/dx. The equation thus shows a structure roughly similar to that 
of the Navier-Stokes equations and has actually appeared in two separate problems in 
aerodynamics. An equation simply related to (1) appears in the approximate theory of 
a weak non-stationary shock wave in a real fluid. This is discussed in Ref. 1 (pp. 146-154) 
where a general solution of (1) is given. The equation is also given in J. Burgers’ theory 
of a model of turbulence (Ref. 2) where he notes the relationship between the model 
theory and the shock wave. Historically, the equation (1) first appears in a paper by 
H. Bateman (Ref. 3) in 1915 when he mentioned it as worthy of study and gave a special 
solution. Eq. (1) is of some mathematical interest in itself and may have applications 
in the theory of stochastic processes. The aim of this paper is to study the general 
properties of (1) and relate the various applications. I wish to thank Professor P. A. 
Lagerstrom and F. K. Chuang for helpful collaboration. 

2. Relationship of (1) to Shock Wave Theory. The solutions to Eq. (1) can approxi- 
mately describe the flow through a shock wave in a viscous fluid. They can be related 
to the shock wave in several ways. In Ref. 1 an approximation based on the Navier- 
Stokes equations for one-dimensional non-stationary flow of a compressible viscous fluid 


gives 


dw ow _ 4,8 
at + én 3 Ox (2) 


for w = excess of flow velocity over a sonic velocity where 8 = (y + 1)/2, »* = kine- 
matic viscosity at sonic conditions. Eq. (2) is reduced to Eq. (1) by Bw = u; 4/3 »* = », 
In this paper a different discussion, intended to illustrate the production and mainte- 
nance of a shock, will be given. In the cases in which we are interested we can say that 
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v is a small parameter, a statement which will be made more precise later. Thus, we 
study the case vy = 0 
ou du 2 
4 - (3) 
in order to see the underlying mechanism of propagation. Eq. (3) is similar to the non- 
linear equations for propagation of waves of finite amplitude in one dimension (Ref. 4, 
p. 482) 
0 


‘2 + (u+c) oh + u) 
(4) 
‘2 + (u—c) 2 \ —u) = 0, 


I 


where w = [°, (dp/dp)'”” dp/p; c = (dp/dp)'”* = velocity of sound 

p = pressure, p = density. 
Eq. (4) reduces to (3) exactly if c = 0. Thus we have as a model a fluid in which w is 
transported by the fluid motion itself (i.e. with a velocity wu). According to the theory 
of characteristics for first order equations, the projection of the characteristics of (3) on 
the (z, 4) plane are straight lines whose slope is 


a =u. (5) 
In addition Eq. (3) interpreted geometrically states that u is constant in the character- 
istic direction. Thus uw corresponds to the Riemann invariants u + w of (4). It follows 
that each characteristic is straight over its length and carries a definite value uv. There- 
fore knowing the values of u at any particular instant the solution for the future may 
be found by following the characteristics. However this process may terminate after a 
finite time when two characteristics intersect. For example this must happen if for 
2, < 22, u(x, , t) > u(x, , t). This phenomenon is the steepening of the wave front for 
waves of finite amplitude known from the study of (4). 

We can now regard the viscosity as a mechanism for preventing the formation of 
discontinuities. The viscous stresses depend on changes of rates of strain so that the 
viscosity » appears in Eq. (1) multiplying a term of higher order 0°u/dx*. This gives 
(1) the nature of a diffusion equation where velocity (actually momentum) is the quantity 
diffused. It is, of course, very much different from solid friction as expressed in first 
order damping terms. The characteristics of (1) are different from those of (3) for any 
vy > O and they are given in the (z, t) plane by the curves y/(z, t) = const. where 


vy, = 0 or ¢ = const. (6) 


These characteristics occur as a double set which can be regarded as the limit of a pair 
of characteristics indicating a high signal speed. Thus here the speed of signals is infinite. 

Hence Eq. (1) shows the typical features of shock wave theory: (i) A non-linear 
term tending to steepen the wave fronts and produce complete dissipation, (11) A viscous 
term of higher order which prevents formation of actual discontinuities and which 
tends to diffuse any differences in velocity. 

3. Relationship of (1) to Turbulence Theory. Eq. (1) is related to turbulence theory 
as a mathematical model. The similarity of the Navier-Stokes equation to Eq. (1) is 


1951] A QUASI-LINEAR PARABOLIC EQUATION IN AERODYNAMICS 227 
responsible. Both contain non-linear terms of the type: unknown function times a first 
derivative; and both contain higher order terms multiplied by a small parameter. The 
problems in turbulence are not very well defined but in most theories one is interested 
in some kind of spectrum, the feeding of energy through the spectrum and dissipation 
of energy. The model equation contains the non-linear terms and viscous terms vital 
to a study of those topics. Following Burgers, we regard Eq. (1) as a model for decaying 
free turbulence; he discusses other cases in much detail. 

The mathematical (rather tian physical) aspect of the model can be emphasized as 
follows. We study as before the underlying wave propagation for vy = 0. For the usual 
turbulence theory we are dealing with flows in two or three dimensions so that in addi- 
(ion to equations like (1) a kinematical or continuity restriction must be added. For 


an incompressible fluid in two dimensions we have 


u, + uu, +vu, = — 1 Pe , (7a) 
1 

v. + w, + w, = — » Pe : (7b) 

uz tv, = 0. (7c) 


The continuity equation (7c) is a statement that only transversal waves are present in 
the flow field. However as discussed in Sec. 2 a typical feature of (1) is the longitudinal 
steepening effect (i.e. steepening in the flow direction). It does not seem possible to have 
waves of that type for a system like (7) but the steepening must be transverse to the 
local flow direction. The underlying structure of (7) is given by the characteristic sur- 


faces V(x, y, t) = const. which now satisfy the equation 


w+ W(v, + wv, + v¥,) = 0. (8) 


> | 


The vanishing of the first factor gives a double set of surfaces ¢ = const. corresponding 
to the characteristic cones which have degenerated into planes, and propagation of 
pressure signals with infinite speed. The vanishing of the second factor gives stream 
surfaces. If any steepening effect occurs it will have to be related to these surfaces. In 
this case an invariant on these stream surfaces is the vorticity &(x, y, t) which satisfies 
the equation 


&, + ué, + vt, = 0. (9) 


Thus, as pointed out by Burgers, the situation in the actual case is complicated very 
much by the kinematical restrictions. 

The quantity wu in the above model is, of course, a measure of the turbulence. The 
model is completed by some kind of statistical analysis based on Eq. (1) or else Eq. (1) 
is considered as a stochastic differential equation subject to random boundary values. 

4. General Properties of Eq. (1). It may be expected that Eq. (1) is similar to the heat 
equation 


06 0 
af on (10a) 
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“ad 


and to the heat equation in a moving medium 


06 00 
Uf stig s 
ot Ox Ox 


06 


(10b) 


with respect to the type of boundary value problems which are sensible. In this section 
some comparisons will be made and other general properties of (1) will be studied. 

An energy equation for (1) can be found by multiplying (1) by wu and integrating 
over a spatial domain, for example (2, < x < 2,). This gives 


pte 94 


| 9, o l ; Si ) 
2). ay ) da +3 iu (te, t) — ula, , db} 


(11) 
Ou Ou \ r=* fau\? 
= W\ulze , 2) — (22,0 — ule, ,Oo=—(%,,07 —»v | \2. dx. 
t Ox , " C2 1 J : Ox 
The various terms in Eq. (11) have the following meaning: 


2 4,,2 
ou 


l , Sere : 
= | 5, ax = total rate of change of kinetic energy in system, 
a vr, Cl 
Po, ae Lal 
g (w(t » t) — u(x, , t)} = net flux of kinetic energy out of system across boundaries, 
e 
au) Ou | ; : 
u- ‘) —\u- > = rate of work done on system at boundaries, 
Ox) ma Ox oa ’ 
ate \2 
Ou Bios : , : : P 
v rn dx = total dissipation of energy by viscosity in system. 
Cu 


The non-linear term in (1) provides a means of feeding energy into the system across 
the the boundaries. We can have a steady-state solution to (i) in an infinite 


domain (— © < x <o@) with an energy balance. 
0 \o 
EP oe F ou \” 
~(ui — Ww) =p | —] dx, (12) 
3 Jaa Ox 
where vu; = u(—©™, t), uw = u(+o, t) and u, > uw. A linear equation like (10a) or 


(10b) can have no such (bounded, non-zero) solution in an infinite domain. For such a 
steady state the total dissipation, as given by (12), is independent of the value of ». 
The steady solution is actually 


u(ae — 2) P 
um —y, tanh ————_, (13) 
2v 
where uw, and 2, are constants. Thus u. = —wu, . This gives the steady flow through a 


shock wave and shows how the non-linear terms are responsible for a change from u > 0 
(supersonic) to u < 0 (subsonic). As v > 0 the steep front of (13) at x = x, approaches 
a discontinuity which corresponds to the shock wave in a fluid where v = 0. This shows one 
reason why the conditions on continuity of solutions for vy = 0 have to be relaxed if 
the solutions are to correspond to reality. Using (13) it is also possible to give steady 
state solutions for finite domains (0 < x < I) which have regions of rapid transition 


either in the interior or adjacent to the boundary. 
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A translation property of (1) is also of interest. If we consider (1) in a coordinate 
system moving in the positive x-direction with a constant velocity U, defined by 





z=2-— Ut 
14 
ge (14) 
we obtain 
au - OU ou 
=+(u—-U) ==, (15) 
ot Ox Ox 
so that w = u — U satisfies the same equation in (z, t) as u in (2, t): 
dw - dw aw 
ot Ox Ox* 


In this sense (1) is invariant under a Galilean transformation. Bateman used the steady- 
state solution (13) as w(z, t) in order to show a shock u(x, t) progressing with a velocity U 
u(x — z, — Ub) 


u= U — u, tanh ——_— - 


2v 


Laws of similarity are also of importance in understanding the joint effects of non- 
linearity and viscosity in (1). For clarity, consider a solution u of (1) depending on the 


following parameters: 


l = significant length; e.g. size of domain, 
Uo = significant initial amplitude, 
y = viscosity. 


Then it is possible to express any solution of (1) in terms of these non-dimensional 


variables: 


» _ ul _ Uol _w _2 
et. ee 7 ae fot 
as 
R = F(R, , 7, 2). (16) 


This relationship is, of course, general. However it should be compared with the corre- 
sponding linear case where it is possible to express the solution as 

u ' 
= F(r, é). (17) 


: $ 
Uo 


These results are easily derived by introducing the dimensionless variables 7 and & in 
the corresponding equations and seeing what is required of u to make the equations 
dimensionless. As one example, compare two solutions with the same / but different 
values of viscosity v, , v2 . In the linear case we can say that the ratio [u2/uo, = Ui/Uo,]e 
In the non-linear case however we can only say that wo, must be 
adjusted so that wu», = Uo,v2/v; if the F is to have the same value, and in that case u. = 
As might be expacted the non-linear equation cannot give much information 


when ft, = 1,,/vo . 


U;V;/Vo 
under linear transformations. 
Various writers (Refs. 5 and 6) have studied the existence and uniqueness of solution 
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of different types of problems for quasi-linear parabolic equations. For initial value 
problems it is clear that only one condition is needed at ¢ = 0, u(a, 0) = u(x). Using 
this and treating Eq. (1) as an inhomogeneous heat equation, they reduce the problem 
to an integral equation. Picard iteration procedures can be used to prove existence 
and uniqueness of solution in the neighborhood of the initial line. The boundary condi- 
tions always involve constant values of u(or u,) on lines x = const. The problem of 
radiation where conditions like u(z, t) = f(t) have to be considered, has not been treated. 
It is much harder to prove existence and uniqueness, or to discover necessary and 
sufficient conditions for this, in the radiation case. Some further remarks about unique- 
ness will be made in Sec. 5. ; 

5. General Solution of Initial Value Problem. The general solution developed here 
applies directly to the case when the initial values are known in some domain 


u(x, 0) = u(x), (18) 
and the boundary conditions are of a simple type. For example, we may have 
u(x, ,t) = 0, u(t. ,t) = 0. (19) 
The solution is supposed to be bounded function having the necessary derivatives. 
The general result is: If @(2, t) is any solution to the heat equation 
00 0°60 
a ess, (20) 
Ot Ox 
then 
: 6, 
u(x, t) = —2Qv— (21) 
] 
is a solution to 
Ou ou eu 
— > = 73° (1) 
ot Ox Ox 


The result can be checked directly by differentiation. However it was derived as follows. 


Let 


fel) ‘ 
c=, 6 = 42,0 (22) 
Ox 
and substitute in (1). Integrating with respect to x we obtain 
ee (2) es a (23) 
ot 2 \dx, Ox” 


where the function of integration is omitted. Like the heat equation (23) is invariant 
under the transformation x — az, t > a’t (a = constant). This suggests finding solutions 
of the form 


o(z, t) = Fie(2, 0} (24) 
where 6 satisfies (20). Introducing (24) into (23) we have 


F'-0, + 5 FPO = vf FG + F700). (25) 
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Hence we obtain the ordinary differential equation for F'(6 
; 1 





1 ary __aF 
2 (4 ~ "a (26) 
which has the solution 

F(0) = —2v log (@ — c:) + ce. (27) 


Therefore, F is the log of a solution to the heat equation and u(z, t) can be expressed 


as in (21). 
Integrating (21) with respect to x we obtain the equivalent relation 


6(x, t) = C(t) exp f. 2 [ u(é, t) at), (28) 
2v Jb 
where b is an arbitrary constant and C = 6(b, t). Without loss of generality b may be 
normalized to the value zero. The initial values are simply related. If 
u(x, 0) = uo(x), (29) 
then 


6(x, 0) = @(x) = Cy exp { x [ Uo(€) at). (30) 


Various representations of the solution to the heat equation may be used. A representa- 
tion suitable for an infinite domain (— © < 2 < @) is given by 
c (x — a) 
Hc, = ———35 | “xX E <_< | @ . 31 
r, t) 2(mvt)'? - exp dot O(E) dé (31) 
Thus given 6, from (30), @ is found from (31) and u(z, ¢) from (21). If an integration 
by parts is carried out in the expression for 6, (21) becomes 
[22 exp [—(a — £)°/40t] exp [—(20)"* S§ wo(n) dnluol@) dé 


a a , 82) 
f22 exp [—(x — £)*/40t] exp [—(2s)"* J§ uo(n) dn] de 





u(z, t) = 


an expression for u(a, ¢) in terms of its initial values. 

The uniqueness of the solution to (1) in the domain (—° < 2 <@) under the 
initial conditions (29) can be discussed as follows. Any solution u(x, t) of (1) defines 
a function @(x, t) according to (28), where C(t) can be found such that @(2, t) satisfies 
the heat equation (20). For each u, C(t) is uniquely determined within a multiplicative 
constant. The initial values of @(x, t) depend only on the initial values of u(z, t) ac- 
cording to (30). Now assume there are two solutions wu, v of (1) having the same initial 
values. If w(x) is suitably restricted @(x, t) is uniquely determined by its initial values. 
But u, v as solutions of (1) are computed from their corresponding @ according to (21). 
In this formula the factor C(t) cancels out. Since @ is uniquely determined, u and v 
are identical for all (x, t) and the solution to (1) is unique. 

6. Examples of Solutions. The first example is that of a shock wave approaching a 
steady state. We choose a step-function as the initial condition where the strength of 
jump is chosen so that the center of the jump will remain fixed. Then we see how vis- 
cosity smoothes even a sharp discontinuity. 
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The initial conditions are 


U(x) = uz <0, 


(33) 
= —y, 2 > 0. 
From (28), putting b = 0 
6.(x) = Cy exp (u,2/2r) x > 0, 
(34) 
= Cy) exp (—u,27/2v) x < 0. 
Upon substituting (34) in (31) we obtain 
iil Co wie (itV\04 Ux a“ (1) f x+ Ut 
(x, t) = Z ex — § 2ch exp |- 2 ae 
\ QP \4y ]\" P\ 2 (vt)? 
(35) 


: Ux C—- Ul 
— exp — 9,7) et ogy 
Hence from (21) the solution is 


u(x, t) 


Qsh(u,x/2v)+ { explux 2verf[(a+1, t) 2(vt)'*|+ exp[—wx/2v]erf[(e—u, b) 2(vt)'”*}} 





=~ "1 O¢h(u,2/2v) + { explu,x/Qverf | (au, t)/2(vt)”] — exp[—wx/2r]erf[(e—wt) /2t)7]} 


(36) 

It is easily verified that as t — 0 the initial conditions are satisfied if we use 
erf(o) = 1, erf(—o) = —1., (37) 
For large values of t/v (and | x | # ut) we can substitute in (36) the asymptotic formulas 

for erf, 
c ~ l 2 2 l 90 
erf z= 1 — —pm exp (—z) + exp (—z)0| 3], z>0. (38) 
wz Zz 


This shows that the approach to the steady state given by Eq. (13) is very rapid. The 
deviations from the steady state die out like exp (—wujt/4v). 
The passage to the limit » — 0 in the solution of (1) is important for determining 


the behavior of discontinuities in the solution to the equation (3) with v = 0 from the 
start. Putting v = 0 in (36) we have, if u, > 0 


u(x,t) = —(sign x)u, , 


so that the initial conditions are preserved and we have a stationary shock wave. The 
invariant quantity across this shock wave is the kinetic energy uj . This is to be taken 
as the rule for treating discontinuities in the solution to (3), from the viewpoint of an 
observer at rest relative to the discontinuity. In addition the flow velocity must de- 
crease in passing through the shock. In the case the initial velocities are not chosen as 

(33) the situation is the same in a moving coordinate system. The solution for a shock 
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tends to a quasi-steady state progressing with a certain velocity U, as indicated in 


Sec. 4. For example if 


u(¢, 0) = ula) = % ,; z <0 
(39) 


= i, , z>0 where u; > Us, , 


we may define a moving coordinate system so that the solution following (33) applies. 
. . t=] 


Let 
u, = (u, — U) = —(u, — UV) (40) 
O 
. I 
U = 5 (us + wu). (41) 
Then according to (15), the solution (36) applies to (39) if x is replaced by x — Ut, 
and u by u — U. The speed of propagation of the shock U is the average of the velocities 


on both sides. 
The second example shows the decay of an arbitrary periodic initial disturbance. 
This corresponds in the turbulence model theory to the decay of free turbulence in a 


box. The initial and boundary conditions are: 


u(x, 0) = u(z), 0O<2<l (41a) 
u(O, t) = u(l, t) = 0, (>. (41b) 
From (21) and (30) these conditions induce the following conditions on 6(2, t) 
: l ar 
(2, 0) = &(z) = C, exp] — = | Uo(é) dé |, (43a) 
2 J 
6.(0, ) = 6,(l, t) = 0 (43b) 


The problem for the heat equation specified by conditions (43) has a unique (bounded) 
solution. We can represent the solution to the heat equation in a standard way by a 
Fourier series in x whose coefficients are exponentials in ¢ 


6(x, t) = Ao + bs exp | -» as tla. cos 7 , (44) 
so that (43b) is satisfied. The coefficients A, , A, are determined at ¢ = 0 as 
1 al C. al 1 z 
A. = | 6.(x) dx = i | exp E ra [ Uo(E) at| dx (45a) 
rr nwa -— 7 ew 
A, = | 6,(x) cos P dz = 4 | exp| — > i uo(é) dé | dx (45b) 


Hence from (21) the solution is 


Pe ur _ exp [—mis"t/F]nA , sin (nra/l) , (46) 
| A+ } , exp [—vn'at/U]A, cos (nrx/l) 
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For large values of time ¢ only the first term in the numerator remains so that 


Qvr A, | x |. xz 
u(r, t) SS — — exp | —y = t| gn (47 
\ / l Ap [- l \ 7) 
This may be contrasted with the solution to corresponding linear problems at large 
times 
r - WL 
6(x, t) = B, exp | -» F | sin, (48) 
where 
Ff Tx 
B, => | U(x) sin — dz. 
. Ll Jo l 


The solutions (47) and (48), are seen to have the same form in dependence on (z, ¢), 
but different amplitudes. The similarity of the solutions is an expression of the fact that 
when the amplitudes are small the non-linear equation behaves like the linear one. 
However the decay process over intermediate ranges of time is considerably different. 
We can find out something about this process by considering special cases. For example, 


consider a simple sine wave 





U(x) = Up Sin =. (49) 


U 


The coefficients are explicitly evaluated in this case as 


l 
C,. f : Uol Wx ee Ul Ul 2 
Ao = & [exo ~ $3 (1 ~ cos) | ae = cs exo | ~ 55 |r (58), (50) 


l 2rv 
9c, f' Uol ( =) | nex | Uy L] (u t) 
= — 4 — — — cos — 708 —— dx = 2C, ex — — I, (|— f 
A, l i, exp | ony 1 cos 7 cos = dx Cy exp a I on)” (51) 


so that (46) becomes 


Avr 2, exp (—vn?xt/P)nI,(uol/2xv) sin (nxx/l) (52) 


u(x, t) = 








9 2 


LU To(ul/2mv) + 2 S52, exp (—v’n?x?t/P)I,(nol/2nv) cos (nxx/1) 


The conditions at ¢ = 0 are satisfied by (52) for 


I,(z) + 2 > I,(z) cos nz = exp (z cos 2). (53) 


n=1 

The significant parameter occurring in the solution is the Reynolds number R, based 
on the size of the box 

Uol - 

R, =—- (54) 

Vv 
The complexity of (52) is due to excitation of higher frequencies because of the non- 
linearity of (1). The same excitation of higher frequencies was noted by Burgers by 
directly splitting the solution of (1) into the form 


u(r, 1) = >> E,(t) sin nx 


n=1 
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and observing the coupling between the EZ, . This can be contrasted with solutions to 
the heat equation under the same initial conditions 

O(x, t) = Up exp | -» a | sin T (55) 
which, for all time, shows only the fundamental frequency. This is typical of linear 
equations as is the linear dependence on the amplitude of the initial disturbance. The 
solution (52) emphasizes the non-linear dependence on the initial conditions. The initial 
amplitude enters through R, . As Ry — 0, u(z, t) = 0(2, t) + 0(R,). This estimate shows 
to what approximation the non-linear terms can be neglected and emphasizes that the 
dimensionless parameter R, (not merely wu.) should be small (Ry < 1). For large values 
of R, all the J,(R,/27r) are almost equal so that large changes in R, (increases in wp) 
produce relatively little effect on the solution. Asymptotically we have 


Ro) ~ exp (Ro (ae) § _ iW} ( L)} ‘ 
I (2) = | ~ “8R, 2m 0 R: (56) 


As the first approximation for large Ry the J,,(Ro/2r) may be cancelled in (52) to give 








4m a exp (—vnin't/P)n sin (nwx/D (57) 


L442 $02, exp (—vn’2’t/l’) cos (nrx/l) 

Equation (57) is an approximation which should describe what happens for ¢ greater 
than some ¢; > 0. Equation (57) can be put in a simpler form by using an identity for 
theta functions (Ref. 7, p. 489) 





u(x, t) = 


; _ _ 2m ys (—)"sin (nrx/l) 
u(x, ) = — TT 2 sinh olna’t/P) (58) 





For large R, , Eq. (58) gives the spectrum, involving all higher frequencies, explicitly. 
The solution is independent of the initial amplitude and the spectrum damps expo- 
nentially with the first power of the wave number for large n. For small n and vn*t/l? < 1 
the coefficients depend on 1/n. The dissipation is proportional to u; and is thus in- 
dependent of n for the first few n. The exponential cut off for large n assures a finite 
total dissipation. For 0 < x < / (58) can be written in another form, 


mn Sh(max/vt) (59) 


x ly 
ataet Pc 
u(x, t) i +2 t ae sh(ml/vt) 


For x/vt large, the series can be approximately summed as 


u(x, t) = 1 anh (! ==) — (1 _ «\ (60) 


a form which shows a steep front near x = l. The general picture presented by the 
above considerations is the following. The initial sine wave (49) shows after the first 
instant a tendency to develop a steep front near x = l, if Ry is sufficiently large. After 
a while this steep front broadens and dies out until at the end only a sine wave remains. 
This sine wave has an amplitude which is smaller than that of the corresponding linear 
problem because of the increased dissipation over the intermediate ranges of ¢. It is 
clear that similar considerations apply to any initial distribution of the same general 





form as a sine wave. 
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7. Concluding remarks. The simple examples which have been worked out are in- 
tended to illustrate some general features of the interaction of non-linearity and vis- 
cosity. The main effects are always the same, namely, steepening of the velocity profiles 
by non-linearity and prevention of discontinuities, diffusion of momentum and dissipa- 
tion of energy by viscosity. Under different interpretations these effects are considered 
responsible for the formation of steep but continuous shock wave fronts and for the 
finite dissipation and feeding of energy through the spectrum. As an example of different 
interpretations consider a velocity distribution with a steep front. The discontinuous 
front of non-viscous flow has infinite dissipation when considered in viscous flow and it 
contributes terms like 1/n to the spectrum. When viscosity is considered from the outset 
the front is steep but continuous, the dissipation is finite and independent:of v; while 
the front is steep there are some terms like 1/n in the spectrum but the spectrum dies 
out like exp (—7) for large n. Another important general feature is the non-linear de- 
pendence of the solution on a characteristic Reynolds number Ry = uol/v, uox/v. For 
low R,(R, «K 1) the non-linearity is not important and the solution behaves like the 
solution to the corresponding heat equation but as /, increases the solution changes 
very much. For large R, it is typical that there are ranges of (x, é) for which the solution 
depends very little on the variations in R, . Part of the problem for the future is a more 
precise determination of the ranges in which the various approximations are valid. 

The same type of result applies to some special solutions in higher dimensions. For 


the equation in two or three dimensions which is analogous to (1) 
qa. +q:-Vqa=vrV'q (61) 
it can be verified that the method of Sec. 5 applies. If @(2, ¢) is any solution of 
,=yp Vv" 0, ,; (62) 
then 
q = —2y V (log 4) (63) 


is a solution to (61). It should be noted that (63) gives an irrotational flow field. 
For future work it seems worthwhile to investigate simple solutions further and study 
the radiation problems. Then the three-dimensional cases may also be studied. 


{EFERENCES 


1. P. Lagerstrom, J. D. Cole, and L. Trilling, Problems in the theory of viscous compressible fluids, Mono- 
graph, California Institute of Technology (1949). 

2. J. Burgers, A mathematical model illustrating the theory of turbulence, Advances in Applied Mechanics, 
Academic Press (1948 

3. H. Bateman, Some recent researches on the motion of fluids, Monthly Weather Review, 43, 163-170 
(1915). 

4. H. Lamb, Hydrodynamics, Dover Publications, New York (1945). 

5. S. Sundaram, On non-linear partial differential equations of the parabolic type, Proc. Indian Acad. Sci. 
(A) 9, 479-494 (1939). 

6. D. Karimov, Sur les solutions périodiques des équations différentielles non-linéaires du type parabolique, 

C. R., Acad. Sci. URSS, 25 (1939). 
Whittaker and Watson, Modern analysis, Cambridge University Press (1940). 


“J 


bo 
wo 
“i 


COMPRESSIBLE FLOWS WITH DEGENERATE HODOGRAPHS* 


BY 
J. H. GIESE 


Ballistic Research Laboratories, Aberdeen Proving Ground 


1. Introduction. The theory of compressible perfect fluids has developed slowly 

because the basic equations are non-linear. Thus it has been profitable to consider 
special examples, such as will be studied here. The present problem originates in the 
study of steady, two-dimensional, isentropic, irrotational flow. If there is a biuni~ue 
mapping of the physical plane onto the hodograph plane, then the equation for the 
velocity potential function can be linearized by a Legendre transformation [5].** This 
draws attention to the case in which the transformation may fail because the mapping 
is nowhere biunique. This suggests the problem to investigate all three-dimensional 
flows whose images in the hodograph space, for Cartesian coordinates and velocity 
components, are curves or surfaces. Such flows are sometimes said to be “‘lost”’ [10] by 
contrast with the nomenclature used here. By analogy with the usage in [5], flows with 
one- or two-dimensional hodographs will be called simple or double waves. The hodograph 
of a flow will be called degenerate when it has fewer dimensions than the original physical 
space. 
The problem can also be motivated as follows. Among the most familiar compressible 
flows are Prandtl-Meyer expansion around a corner or curved wall [8]; Busemann’s 
cylindrical or “swept-back” flow produced by superposition of plane flow and uniform 
flow normal to that plane [2]; Taylor-Maccoll axisymmetric flow about a cone [1, 3, 7, 11]; 
and Busemann’s general conical flows [2]. In these examples the loci of particles of equal 
velocity are planes or straight lines, so their hodographs are degenerate. The question 
arises, Whether this enumeration is exhaustive. 

In this paper the flow will be assumed to be steady, isentropic, and irrotational. 
Characterizations of one- and two-dimensional hodographs will be developed, and 
generalizations will be found for the properties of the examples mentioned above. As 
an example the construction of flows with axisymmetric degenerate hodographs will be 
considered. 

Some aspects of this problem have been considered by Germain [6]. M. H. Martin 
has also made an unpublished investigation along these lines. The construction of all 
axisymmetric flows with degenerate hodographs was studied by Bateman and later by 
Stewart [10]. Opatowski [9] has discussed very concisely the more general problem to 
determine those flows for which the covariant velocity components in some curvilinear 
coordinate system depend only on two coordinates. 

2. Fundamental equations. Compressible perfect flow obeys the equations of motion 


u,; Ou,/dx' = —p ' Op ax" (2.1) 


and the equation of continuity 


O(pu,)/dx’ = 0. (2.2) 


*Received Aug. 15, 1950. Presented to the American Physical Society, Feb. 1, 1947. 
**Numbers in brackets designate papers listed at the end of this note. 
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The x’ (¢ = 1, 2, 3) denote Cartesian coordinates in the physical space, u; velocity com- 


ponents, p density, and p pressure. The convention that every pair of repeated sub- or 
superscripts implies summation over their range has been adopted. For irrotational flow 


Ou; zi = Ou; Ox’. (2.3) 

For isentropic flow 
P/Po = (p/po)” (2.4) 
for certain reference values po and p, , and y = ¢,/c, , the ratio of the specific heats at 


constant pressure and volume. Equations (2.1), (2.3), and (2.4) imply Bernoulli’s 


equation 


suj;u; + a /(y — 1) = ie, (2.5) 
where 
a’ = dp/dp = yp/p (2.6) 


is the square of the speed of sound, and the constant c is the limiting speed of flow. By 


(2.3) there exists a velocity potential function ¢ such that 


Uu; = dp/d2’. (2.7) 

By (2.1, (2.3), (2.4), and (2.6) 
a’ dp éx* = — pu; Ou; dx", (2.8) 

and by (2.2), (2.7), and (2.8) 
(a°5;; — u,u;) 0e/dx* dx’ = 0, (2.9) 

where Kronecker’s delta, 6;; = 1 (0) if 7 = (¥) 7, and where by (2.5) 

a’ = 2 (y — 1)(c’ — u,u,). (2.10) 
x — u; maps a 


3. Degenerate Legendre transformations. The transformation 
three-dimensional region of the physical space onto an n-dimensional region of the 


hodograph space if and only if 


0g az = u; = u,(u),; (3.1) 
the functions u“(x) (a = 1, , n) are independent, and n of the functions u,(u) are 
also independent. Hence the 

rank of || du*/da’ || = rank of || du;/du" || = n. (3.2) 


Disregard n = 0 (uniform flow). Then n = 1 or 2 for degenerate hodographs. Let 


k=o—2'U; (3.3) 
By (3.1) and (3.3) dk/dx’ = —2'(du;/du")(du*/dx’). Accordingly, the Jacobian matrix 
of k and wu” has the same rank as that of u* alone, so k = k(u). (3.3) becomes 
g = x'u,(u) + k(x), (3.4) 
and by (3.1) and (3.4) (a du;/du* + dk/du*) du*/ax’ = 0. By (8.2) 
x’ du;/du* + dk/dp* = 0. (3.5) 
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By (3.5) 
(a 0°u;/dp* du” + 07k/dp" Ou”) Ou"/dx' = —du;/dpu". (3.6) 
By (8.2) this implies 
rank of || a’ Ou; /du* du + dk/dp* ou” | = n. (3.7) 
Hereafter assume that u,(u) and k(u) have been chosen to satisfy (3.7). Then (3.5) 
can be inverted to yield u*(z). 

By (3.1) a point on the hodograph is determined by setting u* = wo . The set of 
points in the physical space which is mapped onto u,(uo) will be called its prototype. 
(3.5) implies 
THEOREM 3.1: If the coordinate axes of the physical and hodograph spaces are parallel, 


the prototype in the physical space of a point, P, of a one (two) dimensional degenerate 
hodograph, H, is contained in a plane (line) parallel to the plane (line) normal to H at P. 


So far @ has only been compelled to yield a degenerate map. For a compressible 
flow (2.9) and (2.10) must also be satisfied. In (2.9) d°g/dx* dx’ is required. By (3.1) 


d°e/dx' Ax’ = (du,/du*)(du"/dz’'), = 


where 0u°/dx’ must be obtained from (3.6). 
4. Simple waves. When n = 1, (3.6) to (3.8) yield 


¢/dx'dx’ = —ulul/(x"u', + k’’), (4.1) 
where primes denote ordinary derivatives with respect to u’. By (2.9), (2.10), and (4.1) 
a’ulul = (u,u’)’. (4.2) 
If s is arc-length measured from some point of the hodograph curve and q is speed of flow 
(s’)* = ulect , (4.3) 
q = uu; . (4.4) 
Now (4.2) implies 
a® = 3(y — 1)(c — q’) = G(dq/ds)’. (4.5) 


Construct a cone, K, with vertex, V, at the origin of the hodograph space and passing 
through the hodograph curve C. When K is developed onto a plane, C will be deformed 
into a plane curve C’ to which (4.5) also applies. Hence C’ is the familiar epicycloid of 
the Prandtl-Meyer flow around a corner. Accordingly, C will be called a conically de- 
formed Prandil-Meyer epicycloid. 


THEOREM 4.1: The hodograph of a simple wave consists of arcs of conically deformed Prandtl- 
Meyer epicycloids. Conversely, a sufficiently small arc of a conically deformed Prandil- 
Meyer epicycloid, on which the direction of the tangent vector varies continuously, is the 
hodograph of a simple wave. 

For the converse, construct a velocity field with the prescribed hodograph. Suppose 
that for A < pw’ < B, u; = u,(u') is an are of a conically deformed Prandtl-Meyer 
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epicycloid, the u, being of class C' on AB. Let f(u') be an arbitrary function continuous 


on AB, and let 


a(n’) = 2, + | Ui(u) f(u) dy, 
vA 

where the x4 are constants. To prevent the curve 2’ = 2'(u') from intersecting itself, 
decrease the interval AB, if necessary. As suggested by Theorem 3.1, through each 
point z'(u') construct a plane normal to u/(u'), and assign u;(u') to every point of this 
plane. By making the interval AB small enough, and by considering only a region close 
enough to the curve x = 2'(u'), a continuous single-valued velocity vector field can 
be obtained. Finally, by constructing in this vector field a family of streamlines close 
to the streamline x z'(u'), astream tube, and hence a flow with the desired hodograph 
will be produced. 

If (4.5) is interpreted as an equation of a plane curve, it is clear that in a simpl 
wave the flow must be supersonic. Discontinuities in the second or higher order derivatives 
of u;(u ) are propagated along prototype planes. Thus the Mach cone at any point, P, 
of a simple wave must be tangent to the prototype plane, I, through P, and the stream- 
line, S, through P intersects II at the Mach angle. 

The reader may verify the following assertions. (1) Sufficiently small ares of any 
curve with continuous curvature can be arcs of streamlines of simple waves. (2) For a 
sufficiently small range of values of u' any one parameter family of planes A,,(u')2”" + 
B(u') = 0 can be chosen to be the prototype planes of a simple wave, provided A,,(u') 
and B ms are of class C’, and proy ided that not all of these planes are parallel. 

As an example for this section, consider a simple wave, W, in which the envelope 
of the prototype planes is a cylinder, S. By Theorem 3.1 the hodograph, H, of W is a 
plane curve. Orient axes so H lies in u; = constant, and let Q° = ¢ — uz. (4.5) becomes 
L(y — 1)[(e — u3) — Q*] = (QdQ/ds)°. This defines an epicycloid obtained by shrinking 
the generating circles of the usual Prandtl-Meyer epicycloid by a factor (1 — u3/c°)'”’. 
As indicated in Sec. 5, W is a swept-back version of Prandtl-Meyer flow. 

5. Double waves. When n = 2 let A = det || 2” d°u,,/dp* Op” + Ok/du" dy” |I, 
where A + 0 by (3.7). By (3.6) 


3p 


Adu’ /dx = (-1 ‘(x” 0°u Ou ' Ou ‘4+. Ok On ' Ou” ) du, On 


where a is summed, but not 8, and where a + 1 and 8 + 1 are reduced mod 2. For fixed 
nu” the solutions x’ of (3.5) lie on a line. Let v'(u) be parallel to this line, so 
vy ou Ou” = (), (5.2) 
and let xj(u) be a particular solution of (3.5). The general solution is 
x = X(u) + rv (p), (5.3) 
where the parameter r is independent of u*. Now (5.1) becomes 


: ‘Ou?’ ‘lau, Ou” (5.4) 


Adu’ /da* = (—1)77P (rve™ + 25) 0'un,/ Op ‘Our + 8’k/du" 


for 6 not summed. By (2.9) and (3.7) 
21 . . x r » bys x 
(a°5"’ — uju;)(du;/du*)(du;/dOp w=) 


3 


-[(rv™ + 275)0° Um Ou” Ou ' +- O°k Ou" * Ou" **] —— 
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with both a and 8 summed. Since r is independent of uw“ 
= e — @ e a & a+ 2 a+la B+ ~ 
(a°5"' — uju;)(du;/du*)(du; /dp")(— 1) ay" u,/Op du*' = 0, (5.6) 
(ad — u,u;)(Ou Ou )(du, Ou”)(—1)**" (250? Um Ou" * aye? 
(5.7) 
+ dk/dp**'du"*') = 0. 
Now let gs, be the covariant metric tensor of the hodograph surface and 6, its 
second fundamental tensor. By definition 


Jas = (du,/du")(du;/dp"), (5.8) 


bus = v" (O° Up /Iu* dy’), (5.9) 


where v” is a unit normal to the surface, i.e. 


vy = 1, (5.10) 

Also 
u, Ou,/On" = ¢ 0qg/dp" = aM dq/du", (5.11) 
where 1 = q/a. Write dqg/du" = q,.. , Where the subscript ,, denotes the covariant 


derivative with respect to uw“ and based on g,g . Tnen (5.6) becomes 
(gas — M°q.09.8)(—1)° "basi gs) = O. (5.12) 
A particular solution of (3.5) is 
x, = —(dk Ou”)g*’ (du, Ou’), (5.13) 


where g“” is the inverse of g,g . Since the Christoffel symbols of the first kind, based on 


Jag , are [aB, y] = (du;/dp")(d°u;/dp" dp”), the second covariant derivative of k becomes 
kop Ok /du" Op” — (dk du”)g? (du, Ou’) (07u, Ou” On’). (5.14) 
Hence (5.7), (5.8), (5.11), and (5.14) imply 


(gap — M°q,a9.8)(—1)° "hk asigs+1 = 0. (5.15) 


(5.12) is a second order quasilinear partial differential equation for three functions. 
To determine u,(u) requires two more equations, which may be obtained by assigning 
a special form to the coefficient tensor gaz — 4/°q,.9,. The resulting systems are classified 
according to the nature of the characteristic curves of their integral surfaces. 

A characteristic is a curve on which the coordinate functions, their first partial 
derivatives, and hence the metric tensor are continuous, while the components of the 
second fundamental tensor may have discontinuities. Suppose du;/du* are known along 
uw = p(t) on u; = u,(u). By (5.9) the strip conditions d(du,,/du*)/dt = 
(d°u,,/On" Ou") du” /dt imply 


bas du’ /dt = v™ d(du,,/du*) /dt. (5.16) 
Then b,, fails to be uniquely determined along n* = u*(t) by (5.12) and (5.16) only if 


(gas — M*q.29.8)(du*/dt)(du*/dt) = 0. (5.17) 
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This defines the characteristic directions du“/dt. By (5.14), if u,(u) are known, then 
(5.15) is a linear partial differential equation for k which also has the characteristic 
directions (5.17). Equations (5.12) and (5.15) will be said to be of hyperbolic, parabolic, 
or elliptic type wherever 2 = det || gas — M°q.a9.e || < 0, = 0, > O. Hereafter (5.12) 
and (5.15) will be assumed to be hyperbolic. Cohn [4] has constructed a double wave 
of hyperbolic type and has also given simple canonical forms for the hodograph equa- 
tions for both the hyperbolic and elliptic cases. 

If s. is the are-length of a characteristic, then (5.17) becomes 

q (dq ds.)° = a’, 

which is identical with (4.5). Hence 
THEOREM 5.1: The characteristics of the hodographs of double waves are composed of arcs 
of conically deforme d Prandtl-M eye r epicycloids, 1.€. of one-dimensional hodograph S. 

By (5.18) the component of g,, along either characteristic is q,. du“/ds, = +a/q = 
+1/M. Hence 
THEOREM 5.2: On the hodographs of double waves the curves of constant speed and their 
orthogonal trajectories bisect the angles between the characteristics. 

Curves x = z'(t) (other than prototype lines) in the physical space are mapped 
onto curves u* = u(t) on the hodograph. It is convenient to know the relation between 
tangent vectors of a pair of corresponding curves. By (3.5), (5.3), (5.13), and (5.14) 


(dx*/dt)(du;/dp*) + [r(bas + k.ag](du’/dt) = 0 (5.19) 
for some r(/). Unless det || rbag + k,2s|| = 0, this determines dy”/dt. Conversely, if 
the curve u* = u*(t) is given on a hodograph surface, its prototype is the ruled surface 

x(r, t) = xo(t) + [r — A(d]r' (0, (5.20) 


where v’ is a unit normal to the hodograph, 2(¢) is defined by (5.13), and A(t) by dA/dt = 
y’ dx}/dt (to make the curves r = constant orthogonal to the rulings). For (5.20) an 
analog of (5.19) is 

(dx /dt)(du;/dn*) + [(r — A)bag + k.as|(du"/dt) = 0. (5.21) 
Since v’ dx'/dt = 0, (5.21) implies 

da‘ /dt = —(du;/du*)g""[(r — A)bas + kas] du’/dt. (5.22) 


In general, the direction of dx°/dt will vary with r along a ruling, so (5.20) need not 
be developable. This raises the question, what curves on the hodograph have developable 


prototypes? Since d°x'/dr° = 0, (5.20) will be developable if and only if 
det || dx°/dr, ax’ /dt, d°x' /drat || = 0. (5.23) 
Since v* and du;/du" are linearly independent, by (5.20) and (5.22) (5.23) is equivalent to 
(Chas + Dk,as) du’/dt = 0 (5.24) 
for some C(t) and D(t) not both zero. By (5.12) and (5.15) 


~ or 


(gap — M°q,a9,8)(—1)**"(Cbasise: + Dk, asipei) = 0. (5.25 
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Hereafter assume Cb. + Dk,.g # 0 for some a and 8. Then (5.24) and (5.25) imply 


(5.17), i.e. u* = u*(é) is a characteristic. 
Next, show that the prototypes of both families of characteristics are developable. 
Let u* = u(t) define one characteristic from each family through a given point P of 


the hodograph. At P, by (5.17) 
2(gae — M°q,09.8) = (—1)***{(dus*"/dt)(duz**/dt) + (dut**/dt)(dus*"/dt)]f 
for some f ~ 0. Then by (5.12) and (5.15) 
bas(dus.,/dt)(du8/dt) = k,ae(dus,,/dt)(du’/dt) = 0, (5.26) 


where ¢ is not summed. Since these have non-trivial solutions dut,,/dt, (5.26) implies 
“and some C = C, and D = D, not both zero. 


(5.24) for uw“ = pw, 
Now investigate the relation between tangents to characteristics and unit normals 
ni to prototypes x° = xi(r, d) of characteristics 4" = u°(t). By (5.20) niv’ = ni axi/ar = 0. 
For some A 4 
n. = A% du;/du*. (5.27) 
Since ni dxi/at = 0, then by (5.25) A*[(r — A)bas + kas] du’/dt = 0. Since z* = 


azi(r, t) is developable, it must be possible to choose A% so that ni does not vary on a 


ruling. Thus A* may be assumed independent of r, so 


A°bagdue/dt = Ak, apdu’/dt = 0. (5.28) 
If for some c(t) and d(Z) 
det || cbag + dk,as || ¥ 0, (5.29) 
then by (5.26) and (5.28) A% = gdut,,/dt for some g. By (5.27) 
ne = g(du;/du*)(dut.,/dt), (5.30) 
i.e. the tangent to u* = yu%,, is normal to the prototype of u* = u* . These considerations 


and elementary calculation yield 


TueoreM 5.3: If for a double wave det || cbag + dk,as || # 0 for some c and d, and tf bap 
and k. 4 are linearly independent: 

(1) The characteristics are the only curves on the hodograph with developable prototypes, 

(2) The tangent at any point of a characteristic is normal, at points of the corresponding 
ruling, to the prototype of the other characteristic through those points. 

(3) The Mach cone at any point of the prototype of a characteristic is tangent to the 
prototype. 

(4) The streamlines intersect the prototypes of characteristics at the Mach angle. 

For the omitted cases, first suppose det || b., || = 0, which includes the case det 
| cbag + dk,.g|| = 0. Then the hodograph is developable, so one family of lines of 
curvature consists of rulings. There exist b, such that b,, = b,bg . Suppose b, is non- 
null. Let n* be a non-trivial solution of b,n* = 0. Then b,n® = 0, so n° is tangent to 
a line of curvature of curvature zero, i.e. a ruling. Since some b,, ¥ 0, the lines of curva- 
ture of the hodograph are uniquely determined. On the other hand, by (5.12) 
(gas — M’q,09.8)(—1)°**basibpa1 = 0, 80 (Gas — M?q,29,3)n"n* = 0, and n° is a char- 
acteristic vector. Hence one family of characteristics must consist of rulings. By Theorem 
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5.1 a plane characteristic is a Prandtl-Meyer epicycloid, not a straight line. Hence 


b, = 0, so bg, = 0, and the hodograph must be in a plane, which may be assumed to 
be u; = constant. Thus prototypes of curves on the hodograph are cylinders, with rulings 
parallel to the x°-axis. These are Busemann’s cylindrical flows. For u* = u, (3.2) and 


(3.4) define the familiar Legendre transformation from the physical to the hodograph 
plane for plane flow, and (5.13) takes the usual form 
Q°bag — Uats)(—1)**" O°k/dUg+; OUgi, = O. (5.31) 


l 


Finally, suppose det || b.3 || # 0, but b,, and k,,, are linearly dependent. By (5.9) 


and (5.14), for u. = pu", d°u;/du, dug and d°k/d u, Jug are also linearly dependent. 
Hence for some f(u; , U2) 
Ok/du, 0Ug = f dO uz/du, Iusz . (5.32) 
By (5.32)(0°u;/du, dus) (dO) du,) = (O°u;/du,0u,)(Of /dug). Since det O° u;/OU, Ug || ¥ 
0, then df/dus = 0. (5.32) yields 
k= B-A'u 5.33) 


where A‘ and B are constants. By (3.5) and (5.33) all prototype lines pass through 
xz = A’. For each streamline S pass straight lines through z' = A’ and each point of 
S. The cone so constructed will be a stream sheet covered by streamlines similar to S. 
Accordingly such flows are said to be conical, a type considered by Busemann. The 


most familiar example is Taylor-Maccoll flow. 


Tueorem 5.4: If for a double wave det || cbag + dk,ag || = 0 for all c and d, or if ba, and 
koe are linearly dependent: 

(1) The flow is conical or cylindrical. 

(2) The prototype of any curve on the hodograph is developable. 


(3) Conclusions (2) to (4) of Theorem 5.3 apply to these flows. 


To every double wave that is neither cylindrical nor conical there corresponds a 
conical flow with the same hodograph. Such general double waves will be called skewed 
conical flows. 

Reconsider the conditions for hyperbolic, parabolic, or elliptic type for (5.12) and 
(5.15). They are sin x >, =, < 1/M, where x is the angle between the velocity and 
the direction of a prototype line. For subsonic flow (5.12) and (5.15) must be elliptic. 
In sharp contrast with plane flow, they need not be hyperbolic for supersonic double 
waves. To see this, consider a supersonic cylindrical flow based on a subsonic plane flow. 
Then (5.15) or (5.31) is elliptic. 

6. Double waves with axisymmetric hodographs. An important class of examples can 
be constructed as follows. Assume u; = wu;(u) is axisymmetric about the u;-axis. The 


hodograph may be represented by 


u, = u(t) cos 8, U, = u(t) sin 8, uz; = w(t) (6.1) 
for some u(t) and w(?). If (6.1) is a curve, two possibilities arise. If wu = 0, (4.5) implies 
w=qg =a = (y — Le’/(y + JD, so (6.1) reduces to two points. If vu and w are con- 
stant, (4.5) implies a” = 0, i.e. (6.1) is the circle g = c, w = constant. This is a singular 


case of a velocity field with constant speed g = ¢ in a vacuum. 
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Next, suppose (6.1) is two-dimensional. Let u' = ¢, uw? = 6. Then by (5.8) and (5.9) 


$i = 4!” + w'. Ji = 0, Joo =U, 
bi = (ww — u’’w')(u? + w’?)~ - (6.2) 
bis = 6; bo. = uw'(u”? + w"”)'” 


where primes denote derivatives with respect to ¢. (5.12) implies 

aulw'u’ — wu’) + w'[a’(u? + w”’) — (uu’ + ww’)*] = 0. (6.3) 
This has the singular solution g = c¢, i.e. a spherical hodograph. It also has the solution 
w = constant, i.e. the hodograph lies in a plane and the corresponding flow is cylindrical. 
Hereafter, suppose w is not constant. With no loss of generality, set / = w. Then (6.3) 
becomes 


a*(uu’”’ — u” — 1) + (ww + w)* = O. (6.4) 


This is a form of the differential equation for the hodographs of axisymmetric conical 
flows, of which Taylor-Maccoll flow or a convergent flow considered by Busemann [3] 
are particular examples. As stated at the end of Sec. 5, to these flows there correspond 
skewed conical flows with the same hodographs. To construct examples, find the function 


k(w, 6). By (5.14) and (6.2) 


ki, = &k/dw’ — 3[log (1 + u’”*)ak/dw, 


kop = Ok/00’ + [uu’/(1 + wu”) ldk/dw. 
By (5.12), (5.15), and (6.2) b,k22 — bk, = 0, so 
ak /ae" + (u u"’)d°k/dw” = 0. (6.5) 


A skewed Taylor-Maccoll flow can be constructed, in the small, by solving the ordinary 
differential equation (6.4), the linear partial differential equation (6.5), and finding the 
prototype lines (3.5). Thus, by relatively elementary processes a class of three-dimen- 
sional solutions of the non-linearized equation (2.10) can be constructed. 
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ON TORSION OF PRISMS WITH LONGITUDINAL HOLES* 


BY 
CHIH-BING LING 


Aeronautical Research Laboratory, Taiwan, China 


Abstract. This paper presents a method of solution, called the method of images, for 
the torsion of prisms having one or more longitudinal holes. The method is applicable 
to prisms of the following four, and only four, sections: a rectangle, an equilateral tri- 
angle, an isosceles triangle and a 30°-60°-90° triangle. These four sections form a group 
by themselves. 

The solution is obtained by adding to the known solution of a corresponding solid 
prism without holes a system of harmonic functions which vanish on the entire external 
boundary of the given section, and besides possess a singularity at the centre of each 
hole. Such a system of functions may be constructed from Weierstrass’ Sigma function 
and its allied functions. 

The solution is illustrated by applying it in detail to a rectangular prism having a 
central longitudinal hole. Numerical results are shown for the special case of a square 
prism. 

Introduction. The torsion of a circular cylinder having longitudinal circular holes, 
with or without a central hole, has been investigated by Kondo’ and by the present 
writer.” In the present paper, the investigation will be extended to a prism, which is 
also pierced by such longitudinal holes. 

Both problems in fact belong to the same general class of torsion problems dealing 
with cylinders of multi-connected sections. Analytic solutions of such problems are 
generally difficult except in some simple cases, and indeed very few solutions have ever 
been found. It appears, however, that certain prisms of this nature can be solved by 
adapting to them the method of images. There are altogether four such prisms, the cross 


sections of which are as follows: 


1) a rectangle, including square as a special case, 
2) an equilateral triangle, 

(3) an isosceles right triangle, 

1) a 30°-60°-90° triangle. 


It is not difficult to show that by reflection about the edges each of the above four 
sections forms a doubly infinite set of images. Furthermore, it can be shown conclu- 
sively by theory of groups’ that these four sections are the only ones which form such 
images. 

The images formed in each case are shown in Figs. 144 respectively. In each figure 
the shaded area represents the fundamental region, or the given section of the prism; 


*Received August 9, 1950. 
1M. Kondo, The stresses in twisted circular cylinder having circular holes, Phil. Mag. (7) 22, 1089-1108 


(1936). 
2C. B. Ling, Torsion of a circular tube with longitudinal circular holes, Q. App]. Math., 5, 168-181 
(1947). 


3W. Burnside, Theory of groups, 2nd ed., Cambridge University Press, 1911, 410-418. 
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The regions marked by positive signs represent the images which are formed by an even 
number of reflections, while the regions marked by negative signs represent those which 
are formed by an odd number of reflections. Note that any two adjacent regions must 
have alternate positive and negative signs. The set of points in each figure represents 
the images due to a given point in the fundamental region. 
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Fic. 1. Rectangle. 




















® 


Fia. 2. Equilateral triangle. 


It may be noted that in Fig. 1 the rectangles may be regarded as grouping them- 
selves into identical rectangles each composed of four adjacent rectangles. Similarly, in 
Figs. 2, 3 and 4 the triangles may be regarded as grouping themselves into identical 
regular hexagons each composed of six adjacent triangles, into identical squares each 
composed of eight adjacent triangles, and into identical regular hexagons each com- 
posed of twelve adjacent triangles respectively. Such identical regions, which are indi- 
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cated by heavy lines, form a doubly infinite set in each case. Obviously, such groupings 
are not unique, but they are immaterial in the present treatment. 
Whereas the present solution is restricted to the four sections as mentioned above, 
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Fia. 3. Isosceles right triangle. 
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Fia. 4. 30°-60°-90° triangle. 


there is however no further restriction theoretically as to the manner of distribution of 
the holes within the section. Naturally, the simplest is the case of a single circular hole, 
especially when the hole is symmetrically located within the section. The next to the 
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simplest is a group of similar circular holes. The latter implies that all the holes are of 
equal radii and symmetrically located within the section so that the properties of all 
the holes are alike, and consequently if the boundary conditions on any one of the holes 
are satisfied, the boundary conditions on all the other holes will be automatically satisfied. 
For convenience, the former will be referred to as the one-hole case and the latter as 
the invariant case. Various invariant cases are shown in Fig. 5. For the 30°-60°-90° 
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Fic. 5. Invariant cases (of more than one hole). 


triangle, it appears that no such invariant distribution (or more than one hole) is 
possible. 

Method of solution. In general the solution of the torsion problem requires a har- 
monic function, say ¥, whose conjugate function is single-valued, such that the function 
W defined by the relation 

VW=y - (2? + y’) (1) 
becomes a constant on each boundary of the given section; x and y being the rectangular 
coordinates in the plane of the section. As no generality is lost, the particular constant 


on the external boundary will henceforth be taken as zero. 
Suppose that the function y is composed of two parts as follows: 


Y=evtnh (2) 
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where ¥ represents the solution of a corresponding solid section without holes. This 
implies that Y> is a harmonic function which possesses a single-valued conjugate function 
and is equal to }(x* + y’) on the external boundary of the section. Consequently, the 
function y, must possess the following properties: 


(1) it is also a harmonic function, 
(2) its conjugate function is also single-valued, 
(3) it vanishes on the entire external boundary. 


Now, the functions Y for a solid prism of the four particular sections mentioned 
above are well known.* In particular, a solution has been given in a unified manner by 
Hay.” His method of solution is also described as a method of images, but is essentially 
based on a different consideration. 

The present problem is thus reduced to find a function ¥, which possesses these 
properties, or more precisely, to find a complete system of such functions for y, so that 
the group of parametric coefficients attached to them can be further adjusted to satisfy 
the remaining boundary conditions at the internal boundaries of the given section. In 
addition, since the function y, for a solid section possesses no singularity inside the 
boundary, therefore the function y, must be of a different system such that it possesses 


Y 








4 + —i xX 
rt 


O 
|-2a-+| 


Fria. 6. 




















singularities inside the external boundary. Such singularities will eventually be excluded 
from the material of the section by the internal boundaries or holes. Thus, if the holes 
are circular, they are preferably placed at the centres of the holes. Naturally, the system 
of functions will be considerably simplified whenever the single hole or the group of 
holes are symmetrically located within the given section. 

It appears that the system of functions for the four sections mentioned above may 
be constructed from Weierstrass’ Sigma function and its allied functions. However, as 
the sections are so different in nature, it seems that no general expression can be put 
forward for all the four sections as a whole. In the following, the case of a rectangular 
prism having a central longitudinal circular hole will be chosen as an illustration. The system 


‘St. Venant, Mém. des Savants étrangers 14, 233-560 (1856), for rectangle and equilateral triangle. 

B. G. Galerkin, Bull. de l’Aecad. des Sci. de Russie 13, 111-118 (1919) for isosceles right triangle. 

G. Kolossoff, Comptes Rendus 178, 2057-2060 (1924), for isosceles right triangle. 

B. R. Seth, Q. J. of Math. 5, 161-171 (1934), for equilateral triangle, isosceles right triangle and 
30°-60°-90° triangle. 

5G. E. Hay, Proc. of London Math. Soc. 45, 382-397 (1939). 
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of functions will first be constructed. Subsequently, formulas for shear stresses and 
torsional modulus will be derived. Methods will also be given for evaluating the co- 
efficients or functions involved in the solution. Finally, numerical results will be shown 
for the special case of a square prism. 

Rectangular prism having a central longitudinal hole. A class of doubly periodic 
functions. Consider in the z-plane an identical rectangle composed of four adjacent 
rectangles, including the fundamental rectangle, as shown previously in Fig. 1. Let the 
dimensions of each rectangle be equal to 2a X 2b and the origin of the coordinates be 
at the centre of the fundamental rectangle as shown in Fig. 6. Now consider the doubly 
infinite set of image points due to a point at the origin. Since the identical rectangle is 
of double periods 4a and 47b, the affixes of the set of image points in the positive regions, 


including the given point in the fundamental rectangle, are given by 


Pin = 4ma + 4nib | 


and 


Qn = 2(2m + 1a + 2(2n + 1)ib = P,, + 2a + 2ib, ) 


(3) 


where m and n are extended to all positive and negative integers including zero. Again, 
the affixes of the set of image points in the negative regions are given by 


P*, = 2(2m + 1)a + 4nib = P,,, + 2a 
and | (4) 
Q*, = 4ma + 2(2n + 1)tb = Pian + 2b, 


where m and n are extended to the same range of values as before. 
A function with a logarithmic singularity at each such point will be defined by 





g,(2)o2(z) 
7(2=> — _— — 5 
Wole) log = o*(z)o%(z) ’ (5) 
where oc; , o2 and o*%, o% are functions of z defined by the following doubly infinite 
products: : 
o,(z) = an. (1 = = exp (i +- se), 
— IT. (1 ~ i) o— (- ™ fr). 
r (6) 
a= TI (1-ze) (5 ae) 


Q*, 202 : 


in each of which the double multiplication is extended to all positive and negative 
integers of m and n including zero, except that in the first the pair of simultaneous zeros 
is omitted. This omission is distinguished by adding to the product sign an accent as 
indicated. It may be noted that the function o, is Weierstrass’ Sigma function while 
the remaining three functions’ are closely allied to o, . Thus, we have 


‘Cf. E. T. Copson, Functions of a complex variable, Oxford University Press, 1935, p. 378. 
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Wo) = —loge — QU” {log (: ~ #5 a p.. + + spr} 


= flog (1 - ) a 
7 a lls = 0. ” ig 2 20, 
(7) 
+ od {log ( _ x) + pe spa) 
+ ES {toe - e) + é + abe 


where the accent on the summation sign indicates likewise the omission of the pair of 
simultaneous zeros of m and n from the double summation. By expanding the logarithmic 
terms in the neighborhood of the origin, this leads to 


W(z) = —logz — - 5 2", (8) 


k=3 


where, for k W 3, 


c-* 
~~ 


oo 1 « (1 1 1) 
-_ -- 9 
-- e+ ze (A+a-2 (9) 
It is readily shown by symmetry that the coefficient 2, vanishes identically when k is 
odd, and is real when k is even. 

A class of analytic functions with poles of integral orders at the foregoing sets of 


image points will be defined by 





puis (-—1)° dad 
2) = 7 0 ’ 1 
W.2) = Gig We (10) 





where s is the order of poles of the function. When s > 2, the class of functions are doubly 
periodic or elliptic functions. 

The form of the functions thus derived is different according as s is odd or even. 
The results are as follows. The initial function W, is also rewritten for the sake of 


uniformity. 


Woe) = loge — as", 
/ 1 iin! ae, 
W.(z) = 3 > 2”, 
“a n=1 
and | - 
eo 5 = "dee, (8 > 2) 
W ocai(2) = ty > a aa (s = 0) 
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where P ste] 
a, = Qr+e 
7 
and in particular (12) 
Bw 
‘a =-Q,, ‘a, = 0. 
7 
Note that here the coefficient °a, is not defined. It will be reserved for later usage. 
Now split these functions into real and imaginary parts as follows: For s > 0, 
W.(z) = S,(z, y) — 1T.(a, y). (13) 
Also, define a pair of polar coordinates (r, 6) by 
z2=2+ iy = re’. (14) 
We then have, since "a, is real, 
S, = —logr — Zz. ae cos 2n6, } 
= 
| | 
T =6+ 2. “"aor” sin 2n6, 
— 
| 
oe cos 26 = on ~ 
S, = 7, = b “aor COS 2n8, 
7 a= | 
908 . | 
Y COS 2& } 2n on ¢ 9\ { on 
S., = eS > Q>2," cos 2n8, (s > 2) (15) 
r n=0 
in 3 ~ 
sin 280 —\ on —— 
T.. = — a zz: “Qo,” sin 2n8@, (s > 1) 
r n=1 | 
cos (2s + 1)6 =. ond on+ / 
Bay, = a a + p "esair | cos (2n + 1)8, (s > 0) | 


| 
T n=0 | 
| 


sin (2s ais 5 eae ee 
Pon = a. — > "0,47"! sin (Qn + 1)0, (> 0) 


r n=0 


The preceding expressions give the expansions of the harmonic functions S, and T, 
in the neighborhood of the origin. Note that S., is even in both z and y, T>, is odd in 
both z and y, S,,4; is odd in x but even in y, and T,,,, is even in x but odd in y. 

It can be shown from Cauchy-Riemann differential equations that the following 


relations exist. 


OS, OT, OT, OSp ss 
oo _ _STo__g,, He M__p,, 
Ox oy Ox oy 
(16) 
os, oT. a OT, Os, 
i 2 Btn et,,,, Sta Sree Oy. 
Ox oY Ox OY 








1951] ON FORSION OF PRISMS WITH LONGITUDINAL HOLES 255 


Furthermore, it can be shown that for s > 1: 


(1) when z = 0, +2a, +4a,---, 
T,,= 0, Sus, = 0; 
(2) when az = a, +3a, +5a,---, 
S. = const S., = 0, T's041 = 0; 
(17) 
(3) when y = 0, +2b, +4b, ---, 


(4) when y = +b, +3b, +5b,---, 


S. = const., Ba, = %, Sea, = 0. 
The two real constants for S, are in fact identical and will be denoted by °a, . It 
is found that 
a. = 20(2a + 2bi), (18) 
where g(z) is Weisrstrass’ elliptic function of double periods 4a and 407. It is seen that 
the functions S,, all vanish on the external boundary of the given rectangle, i.e., at 
x = +a and y = +b, except S, which becomes a constant. Therefore, the system of 
functions is useful in constructing the function y, . 
The solution. For a solid rectangular prism of cross section 2a X 2b, the function 


> is known as’ 


(19) 





_ gpo( 2 > S _(—1)"_ cosh (2n + 1)rx/2 (2n + 1)ry 
46 ( p> (2n + 1)* cosh (2n + 1)za, ts sana 2b . 


the origin being at the centre of the rectangular section. With the aid of the expansion 


Tv 





= 1 9 
cosh px cos py = p> (2m)! (pr)"” cos 2mé (20) 
the above function may be expressed in polar coordinates as follows: 
Yo = fo t+ Do fear” cos 2n6 (21) 
n=1 
where, for n > 0, 
, ] (2m + Ebel" (2m + 1)ra 
> — h* A = —- = 22 
Fon don + 5 51 on man d¢ 1) {Om sech oy} (22) 
in which 6mn = 1 or 0, according as m = n, or m # n. 


Now, we construct 
7A. E. H. Love, Mathematical theory of elasticity, 4th ed., Dover Publications, 1944, 317-318. Note 
that the constant is modified so that Yo = }(x? + y*) on the boundary. 
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vi = >> Ao, S2,(x, y) — “arAs , (23) 


I 
where A>, are arbitrary constants to be determined. The initial function S, is rejected 
on the ground that its conjugate function is not single-valued. The function y, thus 
constructed evidently meets all the requirements as outlined previously. 

To adjust the remaining condition on the internal boundary or the central hole, we 
have in terms of polar coordinates, 


l 2 
y= — = @wt+y)twt 


+ > (soar + A,r" — rr" > "Alo As.) cos 2né. 


Hence the boundary condition on the rim of the circular hole, where r = X say, is satisfied 
provided that for n > 1, 


or 
Big ee KR YO FT Me As, « (25) 


s=1 
The value of Y on the rim of hole then becomes a constant, say Vo , as follows: 
Y= L y2 + fo — > a, 1, (26) 
0 = fo oaAo, - 2 
2 7 
The system of linear equations in (25) may be solved by successive approximations 


as follows. Write 
A; = 2, A? (27) 


where 


and, by iteration, 


Ax =X" Ya, As”. (28) 


Naturally, the validity of the solution depends upon the convergence of the series 
(27). From physical considerations alone, it seems likely that there will be convergence 
as long as the rim of hole does not touch the external boundary, i.e., when 


A < min (a, b). (29) 


Torsional modulus. The torsional modulus’ of the prism is given by 
H=2 | |v de dy + Qed , (30) 


®R. V. Southwell, Theory of elasticity, 2nd ed., Oxford University Press, 1941, p. 323. 
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where the double integration is extended over the entire rectangular section, excluding 
the hole. 
There is no difficulty in evaluating the following integrals: 





+ »b pa eh pT 
| y, dx dy = 4 | | WY dx dy — 2 | | Yor dé dr 
(31) 
2 es ale ‘ 1/4. = 1 (2n + 1)ra 
— f ij 3 e F — = =» — é 
3 ub(a” + 5bd*) Tr fo 5 (4) b p> (n+) tanh ob 3 
ep ab pa r ae 
I (2? + y) dx dy = 4 | [ (x? + y’) dx dy — 2 [ [ r dé dr 
Ja “0 ~0 “0 “0 
(32) 


= ; ab(a’ + b’) — 5m. 


But in evaluating the integral [f y, dx dy the foregoing method fails owing to the fact 
that ¥, possesses a singularity at the origin. However, it may be evaluated by means 


ms - 


}> tx 














= 
Fig. 7. 


of contour integration based upon a corollary from Green’s theorem. Suppose that F 
is a function of z regular in a domain S which is enclosed by a contour C. Then we find 


[[ af 48 = sé ff F d, (33) 


dd 


where Z is the conjugate of z; the contour being taken in a counter-clockwise direction; 
Now, put F = W,(z) where s > 1. Then 


[[ wn@as=- 5 [ woe. (34) 


Thus by referring to the contour C in Fig. 7, for s > 2, 


e 


[ | S., dx dy 


RP. |] W2.@ a8 = - Fas iy RP. | Wee) ae 


II 


(35) 
4 


=— (Qs = Des 2) T0500; b) + 722,r°. 
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In particular, for s = 1, 
If tien « 5 RP. [ iW.(2) 2 = —4T0(a, b). (36) 
Consequently, we have 
// y, dx dy = —A,{4T,(a, b) + (4ab — md’)"ar} 
(37) 


4 "a a | 
_ yA Wo 12s — 2) 8 b) —_ 7Qs,r if 


Hence the torsional modulus is equal to 


6 fn 
H = I > ab? — : mr* — (4) b* | — wae” a = oe ws 
3 2 n=0 (2n + 1)° 2b 
(38) 
— 2A,{4T,(a, b) + (4ab — rd*)°a.} — 8 Ps Da(Ds - re i) oe 2s 
The resulting twisting couple is given by 
T = urH (39) 


where uz is the modulus of rigidity of the material and 7 is the angle of twist per unit 
length of the prism. 

Stress components. The non-vanishing stress components of the prism are two shear 
stress components given in terms of rectangular coordinates by 


OV ov 
Z, = wr — Z, = —ar— 
ao e » ur (40) 
or in terms of polar coordinates by 
ov ov 
Z,== Zy = —ur 
r 06’ ‘ OY ar (41) 


It is now a straightforward matter to calculate the stress at any point in the prism. 
In particular the shear stresses on the external and internal boundaries are as follows: 


f 8b cosh (2n + 1)rx/2b , ) 
, aes 4 ij = : 2 ee ee ee - otio a { r } 
[Zz],~s 2u y T a=0 (2n + 1)° cosh (2n + 1)ra 2b ot 8A a Da ilt, o} 
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T\-3 r 2 , Cos 
mw 425 (2n + Then 2b 2b 
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N 
= 
nd 
Lal 
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: ae 1)” Qn + 1)ra (2n + 1)xy 
aul gg 20+ Ley 
; r (42) 
-+- > 8Ao,So041(4, , 
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The maximum shear stress occurs on the boundary at those points which are nearest 
to the next boundary. Due to the presence of the holes, the greatest shear stress does 
not necessarily occur on the external boundary. 

Evaluation of coefficients. Numerical results will be given for the case of a square 
prism of cross section 2a X 2a (i.e., a = b) having a central hole of two different radii 
\ = a/2 and a/3 respectively. 

When a = b, it can be further shown that for k > 1, 








Q4n4+2 = 0, 
1 2(—1)* (43) 
Q4 = (2a)** 1 pi 97 Tsk » 
where 
Cu = ees 9 (44) 


min=—= (m + ni)” 


rT . . . . 9 
There exists an algebraic relation between the coefficients’ o,,. For k > 2, 


(2k — 3)(4k + 1Cu = CiCug + CsCur-s + Ci2Car-12 + °° + Cue (45) 


Ol— 


where, for k > 1, 
Cu = (4k — low (46) 


The values of o,, have been computed by the present writer some time ago.’° They 


TABLE I. CoEFFICIENTS FOR A SQUARE PRISM 























Am a Aon a 
n San Qin(2a)™ Sona”? mt for\ = 3 7) for\ = 3 
0 - -- 8.93704 XK 107? — —- 
1 3.151212 4.726818 0 0 0 
2 4.255773 3.723801 —9.11851 XK 10°7| 3.56902 x 10-*| 1.38992 x 107% 
3 3.938849 4.061938 0 0 0 
4 4.015695 3.984322 1.92418 X 10-* | —2.84698 & 10-* |—4.46471 x 10™ 
5 3.996097 4.003902 0 0 0 
6 4.000977 3.999023 —1.20581 K 10-4] 7.77531 X 10-"| 4.28303 x 107% 
7 3.999756 4.000244 0 0 0 
8 4.000061 3.999939 1.1841 XX 10-°|—2.4495 xX 10-4|/—6.3624 x 10-*! 
9 3.999985 4.000015 0 0 0 
10 4.000004 3.999996 —1.45438 xX 10-*| 1.4594 x 10-% 1.2010 xX 10 
11 | 3.999999 4.000001 0 0 0 
2 | 4.000000 4.000000 2.0472 xX 10-7|—6.7229 xX 10-*|/—2.5589 x 10-” 

















are reproduced in Table I where values of 2,, and f,, are also tabulated. It is noted that 
by symmetry, for n > 0, 
Sans2 _ 0, Anse al 0. (47) 


9Cf. E. T. Copson, loc. cit., p. 360. 
Doctoral thesis by C. B. Ling presented to London University, England (1937). 
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Values of *“a,, are shown in Table II. With these values, the coefficients A,, can 
now be computed from (27) by successive approximations. The results for two different 
radii \ = a/2 and a/3 are shown together in Table I. 


TABLE II. 4"a,,a*"*** FoR A SQUARE PRISM 


4s in = 4 in = 8 in = 12 in = 16 in=20 | 4n = 24 


10-4 |4.35886 * 10-6 
10~? |2.44895 X 1073 


76622 * 10°-7|3.70005 * 10-3 |4.22136 
92402 * 107! |5.84357 & 107? |1.382335 
22282 x 107! |1.94292 107! |7.88560 10-2 |2.42858 x 107? 
cayenne < 107! |2.79896 K 107! /1.89056 10! 9.14617 x 107? 
31 < 1071! |2.36320 & 107 |2.50741 XK 107 |1.82006 x 107? 
8: x 10-2 1.37193 & 1071 |2.18408 « 10-1 |2.29133 x 107 


415.09114 X& 107 |1.63628 x 107° |2. 
8 13.27256 X 107 |3.91222 « 10/1. 
12 |8.29864 X 10-7? |2.88604 < 107! /3. 
16 {1.48002 * 107? |1.16871 xX 107! /2. 
20 |2.11068 *K 10-3 |3.30837 x 107? /1. 

}. 


24 |2.61531 X 10-4 |7.54684 X 10 ced 


To proceed further in computing the torsional modulus and shear stresses on the 
boundaries, values of “a, , 7';,..(a, a) and 7',,,,(2, a) are required. 
It is readily shown that when a = 8, 


a, = 0. (53) 
Now, define similarly a class of analytic functions W%(z) whose initial function is 
W*(z) = —log {o,(z)o.(z2)}, (54) 


where c, and o> are defined before in (6). The expansion of W*(z) in the neighborhood 


of the origin, for s > 1, is 
W%*(z) = aS = W*(z) = : + (—1)’ : "at2”, (55) 
where, for a = J, 
"at = (1) nenr/a(n + : ~ ; ss gn (56) 


¢,+, is defined in (44), which vanishes unless the suffix (n + s) is an integral multiple 
of 4. An algebraic relation between the functions is as follows: For s > 5, 


] ~ — 
= (§ — 3)(8 — 2)B, = B,B,-2 + BsB.-s + BsBsa + +++ + BiaB, (57) 

where, for s > 2, 
B, = (s — 1)W*(z) = (s — 1){ S*(z, y) — iT*(2, y)} (58) 


When a = bB, it can be further shown that for s > 1, 


We .(a + az) = 0, 
We(a + ai) = Sié(a, a), (59) 


Ws_(a + ai) = —iT#_.(a, a) = —}17,,-.(a, a). 
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The two initial functions may be found from 


i(150,)'” 


he ° 


W*(a + at) = 
(60) 

5o4 

~ 32a‘ 


W*(a + az) 


With the aid of (57) for successive functions, values of 7',,,2(a, a) can then be found 
without difficulty. 
To evaluate 7',,.,(2, a), it is noted that for a = b, 
T ...,(x, a) = 2T%,,,(z, a) = R.P.{[QW?,.,(a + ta)} (61) 


In general, the three initial functions of W*(z) may be found from 


{W%(x) — we) . 7 
= lia + Wes ~ M@ + WM, 


W *(z) 
(62) 

ie eee 15 - ‘aie 

Wz) = i s(z)}° + a W x09 | . 


504 
24.4 
64a 


W*(z) = {W%(2)}? + . 
With also the aid of (55) and (57), values of T,,.,(z, a) can then be computed for any 
particular value of x. 


The results are tabulated in Tables III, IV and V. 


Taste III. 7;,-2(a, a) and 7'4.4:(z, a) FOR A SQUARE PRISM 


Tss41(2, a) 











T 1, a a 
r=0 x=02a | r=04a | tc =06a | x = 08a | =a 
| 1.71880 | 1.93656 0.46903 | --0.25063 | —0.43000 | —0.25138 0 
2 —0.50778 | 2.00239 | —0.16957 | —0.48920 | 0.04270! 0.10612 | 0 
3 0.12501 | 





Tar_Le IV. Cross Sectrionat AREA A AND TorstonaL Mopu.us H or a SquaReE PRIsM 





— ——— — = = l 








v/a A/a? H/a‘ | % reduction of A | % reduction of H 
1 3.2146 2.1502 | 19.63 4.40 
3.6509 2.2298 | 8.73 0.36 

Solid 1.0000 2.2492 | 0 0 
| 
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TABLE V. SHEAR STRESSES ON BOUNDARIES OF A SQUARE PRISM 
Z./u7a ony =a | Zo/ura on r = X 
z/a A = fa A = ja Solid | 6° A = fa A = ta 
oe ere er = = 
0 | 1.858 | 1.351 | —1.351 0 0.524 | 0.360 
0.2 | —1.316 —1.315 —1.315 5 0.523 0.359 
0.4 | 1.210 —1.211 —1.211 10 0.518 0.354 
0.6 | 1.014 | —1.014 —1.014 20 0.504 | 0.338 
0.8 | —0.683 | —0.684 —0.684 30 0.488 0.320 
1 0 0 0 45 0.476 0.306 








The writer wishes to thank Mr. T. C. Lee of the Aeronautical Research Laboratory, 
Taiwan, for his assistance in preparing the manuscript. 
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THE TORSION AND STRETCHING OF SPIRAL RODS (I)* 


By H. OKUBO 
(Institute of High Speed Mechanics, 
Tohoku University, Sendai, Japan) 


SYNOPSIS 

In this paper the torsion or the stretching problem for a spiral rod is treated theo- 
retically. The equations of equilibrium expressed in terms of displacements are reduced 
to forms which are independent of one co-ordinate. They are readily integrated for the 
particular case where the helix angle is small, and the corresponding displacements and 
stresses can be expressed in forms which contain three arbitrary plane harmonic func- 
tions, determination of which is dependent upon the shape of the section. As an appli- 
cation of the general solution, the problem for an elliptic section is solved explicitly. 

Two-dimensional problems in elasticity have been studied extensively from early 
times on account of their simplicity in stress analysis and their useful applications in 
many engineering problems. For a similar reason, various problems of axially sym- 
metrical stress distribution have been investigated by many writers. 

In this paper we shall treat the torsion or the stretching problem for a spiral rod. 
The stress distribution for this case is neither two-dimensional nor axially symmetrical, 
and each stress does not vanish in general and consequently the analysis becomes some- 
what complicated. But the problem is not a three dimensional one without any re- 
striction, since if we rotate the co-ordinate axes about the axis of the helix so as to cc- 
incide with the fixed directions with respect to a section which is perpendicular to the 
axis of helix, then the stress distribution referred to the rotating axes is the same in any 
section. 

Starting from the equations of equilibrium expressed in terms of displacements, we 
shall introduce equations which are independent of the position of the section. The 
differential equations of displacements are readily integrated for the particular case 
where the helix angle is small. The corresponding displacements and stresses are ex- 
pressed in forms which contain three arbitrary plane harmonic functions, determination 
of which is dependent upon the shape of the section, and thus we can considerably 
simplify the problem. 

We shall take the axis of the helix as the z-axis, and shall denote the displacements 
in the x, y, z directions by u, v and w, respectively. Then the equations of equilibrium 


. . 1 
can be expressed in the forms 


(A + ») o- + uV*u = 0, 
x 
OA 2 
(A + 2) ay +2Vv =0, (1) 


aA 
(A + ») 5 + uV'w = 0, 


where A, » are Lamé constants and A is the cubical dilatation. 


*Received November 16, 1950. 
1H. Love, The mathematical theory of elasticity, 4th ed., 1927, p. 133. 
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We shall make the following transformation 
a’ + ty’ = e*(x + iy), (2) 


where k is a constant which is related to the obliquity of the helix. 
Let u’, v’ be the displacements in the x’ and y’ directions respectively. It follows that 


ul + wv’ = e'“*(u + w). (3) 
Let us put the displacements in the form 

ul = u(x’, y’) — ay’z, 

vy’ = v(x’, y’) + azz, (4) 


w= w(x’, y’) + Bz, 
where a, 8 are constants. Using the expressions for the displacements in Eq. (4), we 


obtain the cubical dilatation 


yn av jw 8 av 
si ou = ow = OU} ica es kD.(w,) + B, 


_ - - - ” / 
OX oY Oz Ox Oy 
and 
OA AA OA . 
= — > cos ke + —>sin ke. 
oy 


ax Oa 
Remembering the relation in Eq. (3), we have 
Vu = {Viu + k’D,(u) — 2k°D.(v,)} cos kz 
+ {Viv, + k’D,v,) + 2k? D.(u,)} sin kz, 


where the operators D, , D, and V; are 





9 oO re a 0 0 
D, = 9? aa — 22’! a tt” ta — 2 ao wl 1, 
= “84 Ox OY Oy Ox oy 
a rf) a” a 


Ray 5-5, Viera 


__- > =? 
OX oY 


To satisfy the first equation in (1), we have 


, OA ; rs a ; 
(A + yu) ay! + wf Vim + kD,(u,) — 2k D.v,)} = 0, 
(5) 
. OA : a Sy rae 
A+ uaa + wl[Virn + k Die) + 2k Dr(u)} = 0. 
OY 


The displacements which satisfy the above conditions, also satisfy the second equation 
in (1). From the third equation in (1) we have the relation 


—k(\ + w)D(A) + u{Viw, + k’D,(w,) + kw} = 0. (6) 
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(5) and (6) are the equations of equilibrium for this case. Inte- 


Accordingly, Eqs. 
and w, , which are independent 


erating these differential equations, we obtain U, , V; 


Substituting them in the following equations 


ol 
— Ou 
Al = KA + 2c —;, 
Ox 
. Ov 
Pie SRA A Be 
OY 


Z, = A+ (8B — kD{w,)), 


(7) 
= Ou, , OW 
xi, = (2 4%) 
\OY Ox 
- Ow, , 
X’'’ =4 — * kD.(u,) + kv, — ay’), 
Or 
- Ow, 
Yr. = —; — kD,(v,) — ku, + az’), 
\ OY . 
obtain the stresses, which are independent of z 
If we express the equation representing the bounding curve of the section by 
Fiz’, 7) = 0, (8) 


the condition that the bounding surface of the rod is free from traction is satisfied if 


the equations 


oP oF 
x3. 4 + X!, = — kDAF)X! = 0, 
OX OY 
x!. oF » Fe oF — kDAF)Y! = 0, (9) 
OX oy 
saa or 
x Fy: 8 _ ppymz, =0, 
Ox Oy 


hold at all points of the bounding curve of the section. 
Let us consider the equilibrium of a portion of the rod cut by two parallel planes 


perpendicular to the axis of the helix. Since the stresses in Eq. (7) are independent of 
z, the resultant of the shearing stress on each plane is the same in magnitude but gener- 
ally different in direction; and from the condition of equilibrium of tractions, the resultant 
stress vanishes. The effect of the normal traction on each plane is static- 


of the shearing 
to a single force and a couple. The former can be cancelled by taking 


ally equivalent 
the constant 8 so as to satisfy the condition 


[| Z, dx’ dy’ = 0, (10) 


and by a similar consideration, we see that the latter vanishes. Accordingly, the effect 
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of all tractions acting on each plane is equivalent to a couple due to the shearing force, 
and the solution which satisfies Eqs. (5) and (6) with the boundary conditions (9) and 
(10) is the one for the torsion problem of the rod. If we determine a so as to satisfy the 
condition that the moment of the couple due to the shearing force vanishes on each 
plane, instead of employing condition (10), then we have a solution for the problem of 
the stretching of a rod 

First, we shall consider the torsion problem. When k is small, u, , v; and A are small 
quantities since they all vanish for a straight rod, and so if we neglect the smaller quan- 
tities of the second order, the equations of equilibrium (5) and (6) can be written in 


simpler forms as 


OA ML 2 

ag? tp Vi = 0 

OA m ° 

oo +“ vm, =0, 

ay + X , 171 (11) 


Viw, = 0, 


and the third equation of the boundary condition (9) becomes 


eg oF 
x' = cae = G, 
25,7 + } ay (12) 
and the stresses X/ and Y! become 
X'= (5 om ‘ Y’= (2 + aa’) (13) 
A; Wo y'}; s=H ay’ x]. 


Accordingly, when k is small, the shearing stresses X/ and Y{ and hence the torque 
acting on the rod are the same as those for a straight cylinder, but the other stresses 
do not vanish as in the latter case. 

Let us put 

w, = Uf; — fs), (14) 
where f; is an arbitrary function and f;, = f;(¢), fs = f3(0), £ = 2’ + ty’ and ¢ = 2’ — iy’. 
w, satisfies the third equation in (11). Substituting this expression for w, into the first 


and second equations in (11), we have 


A+ Qu dy, yy pe a peg EH 
Xu on? + ae ay titpay? Nt ht th + th; 
(15) 


ik ff — fe + off’ — tf;'). 


I 





_# om , Ou, A+ WIM 
N+ dar? * dx’ dy’ A+ 420 


The particular solution of Eq. (15) is 
mak feactk | Red, 1 =0. 
Setting the right sides of Eqs. (15) equal to zero, we obtain the relations 


Viu, = 0, or Viv, = 0, (16) 








| 
] 
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from which we can put 
h=fit hi + 2'(fe + fo); 


where f, and f, are arbitrary functions of ¢. Inserting this expression for u, into (15), 


v, = if, — fr) + it'(fe — fr.) +i n+B4l fp ar — [ Rat. 


Accordingly, the displacements which satisfy the equations of equilibrium (11) can be 


expressed in the forms 


W=fththtf +a f nears [Rede | — av’s, 


we have 


v’ = if, — fr) + info — fr) + (2p + | f f.dt — [i if | + ar’z, (17) 


w= ifs — fs) + Bz, 


and the corresponding expressions for the stresses are given in the forms 


Xt = WK AKA Dh + fd teh +f) + MEK + TID} +, 
Y= —wi f+ hit pt Wh + fr) +2'(fi + J} +, 
Z, = —2Qu{(1 — p)fo + fo) + WES + Efs)} + A + WB, 








(18) 
Xj, = Din f ~F+e(h-D+@t+ Mh -D+ECKh- ‘P}, 
X! = wli(fi — fa) — ay’}, 
Y? = —u{fs+ fs — ax’} 
where 
, = a. 
A+ u 
As an example, we shall consider an elliptic spiral rod whose section is given by 
z+ ye =1. (19) 
Let us put 
2ff= Art Ai’, 2 =Bot+ Bs’, fs = hf, (20) 
where 
a’ — b° 
h = +B a. 








H. OKUBO 


ing the expressions in (20) into | 18), we obtain 


X!, = Qu{Ao + pBo + 2Bix” + (A, + pB, + kh)(v” — y’”)} +, 


Yi, = —2Qu{A, + (p + 2)Bo + 2B,z” + [A, + (p+ 2)B,\(2’”” — yy} + AB, 
Z, = —2Qu}(1 = B + (1 = p)B, + kh\(2’” aes y’)} + (A + 2u)B, 
(21) 
Xj, = —2yu{2A, + 2p + 2)B, + kh}x’y’, 
x = —plh a a)y’, 


Y2 = —ul(h — a)z’. 


From the boundary conditions (9) and (10) we obtain a system of simultaneous equa- 
tions for the determination of the unknown constants A, , By , A; , B; and B as 


Ay 
£04 2B, + A+ @+ 2B, +5°5 + kh =0, 
a a “ea 
1, 
fe 4 OTE pA, — + DB, — as = 0, 
D b 2ub 
1 3 3 6 | : 
( 2+ JA, + (2, + p+ 8)p, + ( 2+ ay = 0, (22 
a b a b a b 
3 é 6 
(3,44 )a. + (2@$S +244), ~0, 
a b a b 
. = b? : 2 
i= gee Sly - pe 4 ay ~- Ey we. 
4 Qu 
As a numerical example, we shall take the dimensions of the section as a = 2cm, 
b = 1 cm and shall assume the values of the elastic constants as \ = 8.66 X 10,° 4 = 


8.20 X 10° (unit Kg weight per cm’). Substituting these values into (22) and solving 
the simultaneous equations, we find (unit ka) 


Ao => — 0.5786, A, = — 1.6793, B, = 0.1156, B, = 0.4627, 


B = 0.4500. 


In our calculation, the constant 6 does not vanish. Hence we see that when a spiral 
rod is twisted, an axial elongation (or contraction) takes place. Substituting these 
numerical values of constants into Eq. (21), we calculate the stresses on the 2x’ and 
y’ axes, which are plotted in Figs. 1 and 2. In this case, the predominating stress which 
is concerned in torque is X/ and its maximum amount is 1.6 wa. While the predomi- 
nating stress among the other stresses is the normal stress Z, , which attains its maximum 
amount of 5.44 kua at the points 7’ = +2, y’ = 0, as is shown in the figures. Accordingly, 
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the latter can not be ignored compared with the former except in cases when k is ex- 


tremely small. 
rhe shape of the cross section is not deformed by twist in a straight rod, but in the 


XAKE 











Stresses 
Ns 





Fic. 1. Stresses on 2’ axis. 











SNCSSES 


~ 8S ~~ N W 





Fia. 2. Stresses on 7’ axis. 


case of a spiral rod the distortion of the section, as shown in Fig. 3, is caused by twist. 
Next, we shall treat the stretching problem. For convenience of calculations, let us 


put the displacements in the forms 
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ul = u(x’, y’) — yx’ — ay’z, 


v’ = 9,(2’, y’) — vy’ + ox’, (23) 
w= w,(zx’, y’) + bz, 
where 
_ 
: ie 9 & ™ p)B. 


When k is small, neglecting the smaller quantities of the order of k’*, the further 
calculations become quite similar to those of the previous case, viz., Eq. (11) is equally 


k>0 











Fic. 3. Distortion of the section. 


valid for this case and the stresses are represented in the same forms as in Eq. (18), 
except the following stresses 


XL =wWf+tftolh + hytxr(nt MEME + SSDI, 

Y= —Miflt+ ft (p+ Wht fy t2(nrt fo}; (24) 
f — so \ 

Z, = —2u) (1 = p)( fo a Je) + kK(¢fs + fs) acs 2 (3 tc P)B 


For an elliptic section, we shall put 
2fi = Ao t+ A:f’, 2f2 = Bo + Bie’, 2f3 = Cog. (25) 
The corresponding stresses become 


X!, = Qu{ Ao + pBo + 2Byx”? + (A, + pB, + kCo2” — y"4}, 


Yi, = —2Qu{Ay + (p + 2)Bo + 2Bix” + [A, + (p+ Bile” — y”)}, 


Z, 


2 1 
~244(1 — p)Bo + [1 — pBi + kCo|(2” — y”) - 9 (3 — net, 
(26) 
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1951] 
X!, = —2u{2A, + 2p + 2)B, + kCp}x'y’, 
X? = —pn(Cy + a)y’, 
y? = «pil, ~ det. 


From the boundary condition (9) with the condition that the moment of couple due 
to the shearing stresses vanishes, we have a system of simultaneous equations for the 
B, , B, , Co and a, as 


determination of the unknown constants A» , A, , 


(a’ — b°)C, = (a + DB’)a, 


(a’ + b°*)C, = (a — Bb’) {a + (3 — p)kB}, 


Be 4 F By + A, + (p + 2)B, + kC, = 0, 


a 
we | D we + §) kCo (1, 5) koe (1, a 1) “ (27) 
(4, + 3A, + (2, + b? B, + 2 a + b? + 2 b? a’ —_ 0, 
Ao r 
a Be B, — A, — (p + 2B, = 0, 
oe ot ) (w+s e+ 4) kCo (4 1) _ ka (1, a 1) a 
(3, + bh? A, T a’ + b? B, + 2 a’ + b 2 b? a’ — 0. 


If g is the mean value of the axial normal stress over the section, then 


q= — || Z, dx’ dy’ 
(28) 
Ef, rE 
- 3 oa D 3 ai p)B = 2(1 Fie p)Bo = 2 (a Pane b){d ae p)B, + KC] ? 
where EF is Young’s modulus. From Eq. (28), 8 can be expressed as a multiple of q. 
We shall consider an elliptic section of the same dimensions as in the previous 


example. From Eqs. (27) and (28) we have 


Cy = 2.3565k8, a = 1.4139kp, Ay = —1.8851k’, 
B, = 0, A, = —4.8148k8, B, = 1.1783k'8, 
and 
q 





8 = Fi — 1.7673k} 


Accordingly a twist arises when a spiral rod is pulled axially, since a does not vanish 
in our calculations. The normal stresses X/, , Y{, attain their maximum value of 1.5 
k*°EB at the center of the section, and so they are very small quantities compared with 
Z, . The normal stress Z, is distributed almost uniformly over the section and it attains 
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+ 2.356k")E8 and (1 — 9.425 k°)Ep 


its maximum and minimum values of (1 
ends of the minor and major axes of the ellipse respectively. The predominating 


tres 
Stress 


among the shearing stresses is Y’ and it attains its maximum value of 1.5 FEB at both 
ends of the minor axis 
In conclusion, the writer wishes to express his thanks to Miss E. Itagaki, his assistant, 


for her earnest he lp In this study. 








THE EFFECT OF INITIAL DEFORMATIONS ON THE BEHAVIOUR OF A 
CYLINDRICAL SHELL UNDER AXIAL COMPRESSION* 


BY 
P. CICALA 


E cuela SS pervor de Aerotécnica and University of Cordoba, Arge ntina 


The theoretical investigation presented here concerns the effect of certain initial de- 
] 
the buckling under axial compression of thin plate structures. A general 


formations on 
theory is developed and used to make approximate analyses of the complete cylinder and 
urved panel. Two extreme cases of the curved panel are treated: (a) uniform longi- 
tudinal shortening and, (b) uniform longitudinal stress. Particular attention is paid to 
the possibility of inducing buckling in a favorable mode, and thus increasing the load 


which a structure can withstand. 


INTRODUCTION 


The buckling of a cylindrical panel under axial compression has recently been the sub- 
ject of a number of investigations using the non-linear theory of thin plates. This 
theory is valid for deflections of the same order of magnitude as the thickness of the 
pane ] 

The non-linear theory has shown that, for a panel with perfect initial form and 
starting from the configuration of infinitesimal deflections which correspond to the 
critical loads of the linear theory, buckling may evolve through one of several paths. 
While some of these paths are favorable and give loads which increase with increasing 
deformations, others cause the load to fall after buckling. This has been proved quanti- 
tatively by von Karman and Tsien’ (1) and, independently, also by the author (2). 

Experiments have shown that the effect of initial deformations is usually to cause 
buckling in an unfavorable mode. However, Welter (3) has shown from compression 
tests on curved panels, that it is possible to delay buckling by initially bending a panel 
a smaller radius than that at which it is tested. These results suggested the present 


; 


LO 


In\ estigation 


Tur GENERAL EQUATIONS FOR A CYLINDRICAL SHELL WITH INITIAL DEFORMATIONS 


Let t’ be the thickness of the shell and write t = t’/+/12(1 — »*) where » is Poisson’s 
ratio. We call state “0” a state of deformation by which the mean surface takes the 
form of a circular eylinder whose radius is R; let 2 and y be the axial and circumferential 
coordinates, w the radial displacements (positive outwards), leading from state ‘‘0” to 
the actual state “‘A”’ 

Assuming that ?’ is small compared with R and that the displacements w have the 
same order of magnitude as ?¢’, the normal and shearing stresses in state “A’’ may be 
epresented in the forms 

of : of of 


o,=E aye? 7= -E az oy , (1) 


C~- = 


; Oy r* 


*Received August 28, 1950. 
Numbers in parentheses refer to Bibliography at the end of the paper. 
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where £ is Young’s modulus and f is a stress function. Similarly the stresses in state 
“0” can be derived from another stress function f’. If the material is completely elastic, 
the stress variation from state ‘‘0” to state “A” is related to stretching by the equation 


Vif -f) +k =9, (2) 


,- 2w (Zw 1) ( af ) 
‘ Oa \ay” oR Ox dy/ ° 


The shell is held in state ‘‘0”’ by applying a certain distribution of normal pressure 
Et'p. , where pp is, in general, a funnction of x and y; the quantity 


where 


Ley 
p’ =pP-s £ (3) 


will be zero if the shell is free from transverse shearing stresses in state ‘‘0”’. In general, 
this will not be the case; the quantity Hi’p’ represents a net inherent pressure that adds 
to the plate and membrane pressures in any state of the shell. Therefore, the equilibrium 
condition in state “‘A”’, under the applied pressure Et’p yields 

ow of _ (= 1) Of 4 9 dw of 
dx Oy” (4) 


oe a. “ 
Serer oy R 


Ox” Ox dy dx dy 
If the shell is free from stress in state “0”, the functions f’ and p’ are both zero and 
Eqs. (2) and (4) degenerate to the known equations for a shell having perfect initial 


form and no initial stresses. 
If a particular state w = w* exists for which the shell is free from stress when there 


is no normal pressure, we may write 


f’ = k*, po = —fV‘w* 
where k* is the value for k when w = w*. Eqs (2) and (4) then take the respective forms 
Viftk—k*=0 (5) 
a Oe = th ~ SEZ (Ee... 4) Hf ow of 
Ox Oy \ OY R/ dx Ox Oy Ox OY 


These equations have already been presented by the writer (2). They have been used 
by Dei Poli (4) in his studies on the cylindrical panel with initial deformations. 
THE CoMPLETE CYLINDRICAL SHELL WITH INITIAL DEFORMATIONS 


We consider the case of a cylindrical shell with the particular deformations represented 
by 


Si = 0 
(6) 
p’ = p, sinésin g + pz cos 2 + pz cos 2g + p, sin 3é sin yg 
with & = rz/a, ¢ = ry/b, where a and b are constants and zR/b is any integer. The 


first of these equations states that the shell is free from stress in state ‘‘0”. The con- 
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strained sections of the shell are assumed to be so remotely situated that their influence 
on the buckling process may be neglected. 

To obtain an approximate solution to the problem of the shell under axial load, we 
first express the functions f and w in the following non-dimensional forms 


f/? = f, sinésing + f, cos 2 + f, cos 2p + f,sin 3¢sin g — ny’/Rt, 
(7) 
w/t = w, sin —sin g + w, cos 2 + w; cos 2g + w, sin 3€sin g + UW. 
The coefficient », which is positive for compressive axial loads, represents the ratio of 
the mean axial stress to the critical stress ¢, = —2Ht/R derived from the linear theory. 
The last term in the first of Eqs. (7) therefore relates to the mean axial stress. 
To determine the values of the constants f, , fo , fs , fa , Wi , We » Ws and w, , We use 
the Galerkin method and evaluate the integral expressions 
p2a 2b 2a 2b 
[vac f (v'tt+klWidy=[ dz | py. dy =0 (8) 
‘ “0 “0 “0 
where for ¥, and y, we replace successively the functions 
sin £ sin ¢g; cos 2£; cos 29; sin 3 sin ¢ 


In this way we obtain 


Df, = —Bw, + 2w,w. + 2w,ws — 2w.w, , (9) 
328’ f./a? = —88°w, + wi — 2w,u, , (10) 
320° f,/B° = wi + 9ui, (11) 
Dif, = —9B’w, — 2w,w, + 18u,u, , (12) 


Dw, = 28° uw, a B’f, — 2w, fe tess 2u.f; — 2w, fs — 2ws fi ed 2w.fs a 2u,f. > Cy ’ (13) 


168°w,/a*> = 8n8’w, + 46° fo —- Wifi t wife t wifi te, (14) 
1607w;/8" = —wif, — 9sfs +s , (15) 
D,ws, = 18n6’w, + 9B’ fe + 2, fe + 2wefi — 18wsf, — 18wifs + €4 , (16) 
where . 

dm (+t), = (+2): 


a= a/r~/Rt, B= b/rV/Rt, 


c; = a’b’p,/r't’, (@ = 1, 2, 3, 4). 
Taking into account the equation 


du. 1faw\? 1 3 3 = 
an 1 (au) = 5. — v0) = 54-9 4 (17) 
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where wu is the axial component of displacement, the mean axial negative shortening 


l aa Ou 
| — Ax (1S) 
243o O2 


t 


may be calculated. The condition of uniform shortening, requiring that e, be inde- 
pendent of ¢, is satisfied by Eqs. (7), when Eq. (11) is considered. Denoting by ¢ the 


| 


constant value of «. , we obtain 
—eR/t = 2n + (uw; + Sw, + 9w)/8a°*. (19) 


For a given initial deformation and given wave length, the constants c, are first 


calculated from Eqs. (9)-(16) by replacing for the w’s the values w] , wi , wi , w{ that 














1255 





° 
a 
l 





I i i ! | = 3 
257 0-5 1 15 Et/6, 





Fic. 1. Stress-strain diagram for the complete cylinder. 


, 


Line a (w;’ = w’ = 0). Line b (w;’ = w.’ = 1/4). Line c (w,’ = 2, we’ = —1/4). 


correspond to 7 = 0. Then successive sets of values of n, fi . fo, fs,fs, Wi, We, Ws, Ws 
satisfying Eqs. (9)-(16) are to be found. The corresponding stress, waveform and 
shortening can be deduced from Eqs. (7) and (19). 
(a) First form of deformation equation 
As a first step, we consider the case where 

a= 8B, w= w= fg = 0 


and hence disregard Eqs. (12), (15) and (16). In this case, the results obtained for a 
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special value of 8 may be generalized by considering that the parameter 6 disappears 
from the equations if the quantities 
f,/B", f2/B*, f,/B°, w,/B’, w./B’, c,/B", c2/B°, (n/B’) — (2/8*) 


are introduced 
The results of the calculations for 8 = 2 are represented in Figs. 1 and 2. In Fig. 1, 


= g/o, are plotted against e/c, for various values of wf and wi . Considering 


\ alues OI n 


+5 





j l l l l i l 
“12 -08  -O% 


Fic. 2. Displacement components wu; versus W»2 . 








Line a (w)’ = we’ = 0). Line ec (w;’ = 2, wo’ = —1/4). Line d (w,’ = 1/2, we’ = —1/16 
Line e (w,;’ = 1/32, we’ = 1/16). Line f (w;’ = 1/32, we’ = 1/4). 

Line g (w,’ = 1, we’ = —1/4). Line h (w;’ = 1/4, we’ = —1/16). 

Line i = 1/8, w.’ = —1/16 


the graph for a shell of perfect initial form (wf = wi = 0), O is the point at which buck- 


ling begins: for this point, w, = w, = 0 and 7 = 1. The lower branch of the curve through 
B represents the unfavorable path with decreasing stresses and strains accompanying 
the first phase of increasing deformations. The upper branch of the curve through A 
represents the favorable path. If Hooke’s law were obeyed and the curve followed 
the upper branch, the behaviour of the shell would be practically unaffected by buckling. 
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Figure 2 shows graphs of the deflection components w, versus w, for various initial 


values wi and w3 . These graphs are limited to positive values of w, , since the diagrams 
would be symmetrical with respect to the w, axis. The graph for a shell with wf = w} = 0, 
indicates that, for small values of w, and w. , w, = +4w, ; the positive sign corresponds 
to the unfavorable branch of the curve in Fig. 1. A third possible state of equilibrium 
would be represented in Fig. 1 by a horizontal line through 0; this corresponds to axially 
symmetrical buckling with w, = 0. 

From Fig. 2, the evolution of the buckled form of the cylinder may be traced for 
rarious initial deformations, the corresponding stress-strain diagrams being given in 
Figs. 1, 3 and 5°. The graphs of w, vs. w. in Fig. 2 are of two kinds: (a) those which 
approach that branch of the line for a cylinder of perfect initial form which has positive 
slope and, (b) those which tend to approach the —w, axis; in any case they diverge 
from line 0A. From Figs. 1, 3 and 5, it may be seen that the first kind of initial deforma- 
tion leads to an unfavorable path while the second kind indicates a rather favorable 
behaviour with the load approaching the critical load of the linear theory. It will be 
seen later that this latter possibility is fictitious. 


(b) Second type of deflection equation 


Calculations taking into account the w; component of displacement and neglecting w, 
and f, as before, have already been carried out (2) for the case of a cylindrical shell 


with perfect initial form—that is, with c, = c. = cz; = 0. The effect of the w; com- 
ponent does nat modify the results to any great extent. For example, the dotted line k 
in Fig. 5 shows the stress-strain curve for 8 = 2 and wi = 1/32, wi; = 1/16, c; = 0. 


This curve does not differ substantially from the full line e which corresponds to the 
same initial deformation, but with w; and Eq. (15) neglected. 
(c) Third type of deflection equation 


The stress and displacement functions will now be expressed in the forms 





f/? = fysinésing + fy cosésing + f, cos 2E + f, cos 2p — ny’ /Rt, 
(20) 
w/t = w, sin ésin g + w,; cosésin g + wW, cos Z2E + wW. 
The coefficients are determined from Egs. (8) in which we write successively 
vy, = ¥. = sinésing, = cosésing, = cos 2, ¥, = cos 2y. 
By eliminating the coefficients f, we obtain for a = 8: 
4 9 9 © 9 2 
. 8 38° W} 2 C; 38° Wit 2 Cr 
2nB —- 4— = —- —w+—7>4¢u4+—-F te wtoa tet 
’ f _ 8 W, lie 8 . Wry 
(21) 
38° 9 9 ] 9 9 Co 
= — —— (wi — Wr) + 5 (wi + WD) +— - 
32w, ad | Il) 8 (1 II) Wo 


The constants c, , ¢;; and c, depend upon the values w{ , wf; and w of the displacement 
components when 7 = 0. 


2Corresponding curves are marked by the same letter. 
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It should be noted that the solution does not modify if the quantities w, , wy: , W2 
are respectively replaced by wy, , w, , —W2. 


As before, calculations were made for 6 = 2. For the initial displacement com- 
ponents (wf = 1, wi; = 1/32, w) = —1/4), the stress-strain diagram is represented by 


the full line in Fig. 3. Comparing this graph with line g, corresponding to 


7F O@) 


o 
0, 
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Fig. 3. Stress-strain diagram for the complete cylinder. 
Line f (w;’ = 1/32, wry’ = 0, we’ = 1/4), or (w,’ = 0, wy’ = 1/32, we’ = —1/4). 
Line g (w;’ = 1, wy’ = 0, we’ = —1/4). 
Line g’ and vector diagrams (0), (A), (B), (C), (D) of displacement components: (w;’ = 1, wy’ = 1/32, 
we’ = —1/4). 
(wi = 1, wi = —1/4), the effect of the additional displacement term is seen to change 


the form from a favorable to an unfavorable one. 

The evolution of the buckled form will be evident from the vector diagrams (0), 
(A), (B), (C) and (D) in Fig. 3, which represent the displacements for the corresponding 
points on the full curve. In these diagrams, the vector “02” represents the component 
w> , the horizontal and vertical components of the vector “01” represent respectively 
w, and w;; . Diagram (0), corresponding to 7 = 0, indicates that w’ is negligible. How- 
ever, as one proceeds from (0) to (D), wr increases rapidly, while w, decreases. The 
units for displacements are indicated on the diagrams. 

Fig. 4 represents the deflections along a generator y = 0 or y = nb (nis any positive 
integer), for points corresponding to those on the full curve in Fig. 3.° The initial form 

’The position of the line w = 0 with respect to the curves in Fig. 4 is of no great interest: the value of 
wo in the expression for w is determined from the condition that the mean circumferential stress is zero, 





The diagrams in Fig. 4 are drawn to such a scale that the quantity V wi? + wy); is always represented 
by the distance apart of the horizontal lines. 
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in Fig. 4—graph 0—, almost symmetrical about the verticals « = —a/2, a/2, 3a/2, 


ete., corresponds to the most favorable type. This symmetry is destroyed as the load 
increases and, during the buckling, the form approaches the unfavorable type which has 
a shape similar to graph (0), but with opposite signs for the deflections. 





{ 
——N| Q 
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| 8 
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Fic. 4. Evolution of the deformation of a generator during the buckling process corresponding to line g’ 


in Fig. 3:u positive outward 


Analogous results are plotted in Fig. 5. As before, the evolution of the three de- 
formation components is represented vectorially. Starting with the initial diagram (0), 
for which (wf = 1/8, wi; = 1/32, ws = —1/16), the vector diagrams (A), (B) and 
C) show, to various scales, the deformation components, in successive phases of buck- 
ling, when —w, has the values 1/8, 1/4 and 1/2. The corresponding points in the stress- 
strain diagram fall approximately on curve e. Starting with the initial diagram (0’), 


for which (w{ 1/4, wy = 1/32, w$ = —1/16), in successive phases when —w, has 
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the values 1/2 and 1, the conditions represented by diagrams (D) and (E) are attained. 
rhe corresponding points lie practically on the same curve e. 

These results show that the behaviour of the cylinder, as calculated when the com- 
ponents W, , Wr; , We are considered, always approaches the most unfavorable path 
between the one determined neglecting the w, component and the one neglecting the 


u component. 
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Fic. 5. Stress-strain diagram for the complete cylinder. 
Line e (w,' = 1/32, wy’ = 0, we’ = 1/16), or (wy’ = 0, wy’ = 1/32, we’ = —1/16). 
Line h, i (w;’ = 1/8 + 1/4, wy’ = 0, we’ = —1/16). 
Points A, B, C and vector diagrams (0), (A), (B), (C): (w’ = 1/8, wy’ = 1/32, we’ = — 1/16). 
Points D, E and vector diagrams (0’), (D), (E): (w;’ = 1/4, wy’ = 1/32, we’ = —1/16). 


THe CurveD PANEL 


The problem of a curved panel which is compressed in a direction parallel to the straight 
sides is now considered. It is assumed that the straight sides y = 0 and y = b of the 
panel are simply supported and free from normal stresses and that the strain along these 
edges is constant (as if restrained by flexurally rigid, laterally and torsionally weak, 


ixially very stiff side-stringers). The boundary conditions for y = 0 and y = b are 
therefore 
Ow of af 2nl 
w= =—-—5; = 0; £=—2. (22) 


oy = Ox” oy” R 
Che parameter » represents here the ratio of the strain along the straight sides to the 
critical strain obtained from the linear theory for a complete cylinder of which the 


panel may be assumed to be a part. 
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For the curved edges x = 0 and x = 2a, we assume perfect clamping and hence 
Ow 
w=—=-0 (23) 
OX 


for these edges; in addition, we assume that these edges are free from shearing stress 


and write 
of 
——a— an @, (24) 
Ox OY 

The problem is considered for two extreme cases: (a) uniform longitudinal shortening 

and, (b) uniform longitudinal stress. 


(a) Panel subjected to uniform longitudinal shortening 
For the stress and displacement functions, we write respectively 


f/3nt = (f,; cosé + f. cos 2E + fs) sing — ny’ /3nRt, 


w/3rt = [w,(1 — cosé) + w, (1— cos 2£)] sin ¢, 


where & = rx/a, ¢ = ry/b. These functions satisfy all the boundary conditions given 
in Eqs. (22), (23) and (24). The condition of uniform axial shortening will be satisfied 
in an approximate form by writing* 

| (ec, — e-)sing dy = 0, (26) 
where —e, is the axial shortening calculated according to Eqs. (17) and (18) and —e 


is the edge value given by —e = 2nt/R. Thus we get 
a’ f,/B’ = —2Qw; — 8w2. (27) 


Furthermore, to determine the constants f; , fo , f; , W: , W2 , we use the Galerkin method 


and evaluate the same Eqs. (8) in which, in this case, the field of integration is defined 
by the panel contour and for y, , ¥. we write successively the functions 


(1 — cos &) sin g, (1 — cos 2€) sin ¢. 
Thus, assuming f’ = 0, we obtain 


Df, = Bw, + 8w{, + 4w uz , 





e 5.2 — 3 2 
D.f. = Bw. + 8w2 + 8w,W2 — 5 M1 
. cos sl —"" ; (28) 
2nB'w, = Dw, + 2a°(w, + wW2)/B + (8 + 16w, — 4u2)f; 
+ (32w. — 12w,)f. — 8wif3 +c ; 
52 ee ' 2 2 , . 
2nB’w, = Dw, + a (w, + w.)/28° + (8B + 16w, + 8w,) fo — Wifi — 8W2fs + C2 , 
‘By substituting an arbitrary function of ¢ instead of f; sin ¢ in the expression for f and using ez = € 


to find the function, it is possible to satisfy exactly the condition that the shortening of all generators is 


constant. 
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where 
a BY a 28\? 
1 8 + a ’ 2 28 + a 
By elimination of f, , f. and f, , from Eqs. (27) and (28) we obtain two equations in 
W, , W. and n, and hence a fourth degree equation in w, and w, . This furnishes corre- 
sponding pairs of w, and w, and these may be used to determine the ratio 7 and the 
axial shortening. 
The load sustained by the panel, expressed by 


eb 


wnt fal 


varies with the section considered. We therefore define an equivalent mean stress o as 
the stress which, when multiplied by the mean shortening and the volume of the panel, 


























\ 
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Fic. 6. Panel subjected to uniform shortening. Shortening parameter 7 versus displacement 
component wu; . 


represents the virtual work done by the forces applied to the panel for the displace- 
ments corresponding to the actual deformation. It may be shown that, according to 
the condition of uniform shortening, this gives 


1 we 
c= a7 | N dz 


and finally 


o/o. = 9 + 3fs/8° (29) 
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| | ! | l 
Fic. 7. Panel subjected to uniform shortening. Stress-s train diagram. 
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“0-4 
Fic. 8. Panel subjected to uniform shortening. Displacement components w2 versus w . 


Numerical calculations were made for a = 6 = 2 and the results are plotted in Figs. 
6, 7 and 8; corresponding points and lines in the figures are marked with the same letters. 
The curves m,m,. and m,m, in the figures apply to a panel with perfect initial form, 
and the arrows on the curves in Figs. 6 and 7 denote points where the deflection com- 


ponents w, and w, are both zero. Of the four branches of the curves starting from these 
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points, only m, gives a favorable characteristic. Favorable behaviour of the panel can 
only be expected for initial deformations whose representative points fall on or near 
the branch m, in Fig. 8. This indicates that favorable initial deformations are directed 
outwards, with maximum values near the center of the panel. 

Curves n, and n, depart from the favorable path m, and this phenomenon is similar 
to that observed for a complete shell. Curve n. departs immediately towards the large 
deflection regime with little increase in the end load. However, the initial part of the 
































1+ 
ST 

| [ l L | 
-2 -1 O 5 W 


Fic. 9. Panel subjected to uniform shortening. Shortening parameter » versus displacement 
component wv, 


n, curve corresponds to increasing loads, as shown in Fig. 7. After reaching a certain 
point, the curve in Fig. 8 suddenly veers in the direction of the m, branch. This phase 
occurs with decrease of load but, before it occurs, the load reaches such large values 
that stresses may be obtained in excess of the elastic limit of the material. For this 
reason, the portion of the n, curve in Fig. 7 which corresponds to decreasing loads is 


not shown. 


Similar results were obtained from calculations made for a = 6 = 4, the various 
graphs being given in Figs. 9, 10 and 11. From the points corresponding to w, = w, = 0, 


which are indicated by arrows in Figs. 9 and 10, the favorable path is p, . Line 7, sud- 
denly veers towards p, but, before it does so, very high loads are attained. Line g, soon 
approaches the branch p; in Fig. 11 and indicates unfavorable initial deformatioas. 

The dotted curves in Figs. 6 and 9 correspond to w, = f. = 0, when the second and 


fourth of Eqs. (28) are left out. 
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Fic. 10. Panel subjected to uniform shortening. Stress-strain diagram. 
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Fic. 11. Panel subjected to uniform shortening. Displacement components wz versus wv; . 
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(b) Panel subjected to uniform longitudinal stress 

The new boundary condition, replacing the one expressed by Eq. (26), is 
af : 

— = —2nt/R for z=0 and 2 = 2a. 

oY 


The parameter 7 represents here the ratio of the applied uniform stress to the critical 
stress obtained from the linear theory for a complete cylindrical shell. The new con- 
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Fic. 12. Panel subjected to uniform longitudinal stress. Stress parameter 7 versus displacement 
component wv . 
dition, as well as those expressed by Eqs. (22), (23) and (24) are found to be satisfied 


by writing 


f/3rt? = [f,(1 — cos é) + f2(1 — cos 2€)] sin g — ny’ /3nRt, 


(30) 
w/3rt = [w,(1 — cos £&) + w.(l — cos 2€)] sin ¢. 
Proceeding as before, we obtain 
Dif, + 2 ro (fi + fo) = —B’w, — 12wi — 16w; + 4u,w, , 
“we ; . 2 i % 
D.fo + 33 (fi + fo) = —Bw. — 12w, — 9 — 8w,v, , 
(31) 
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> ce td "< ¥ 
2n8w, = Dw, + 2 = (uw, + w.) — (8 + 24u, — 4.) f, — (82w. — 4u,)f2 +, , 
'e) 
> ck 2 9 ' ) 2 2 
2nfw. = Dw, + 553 (Wi + we) — (8 + 24u. + 8w,)f. — (8uv2-wW)fi te, 
28 
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Fic. 13. Panel subjected to uniform le 


mgitudinal stress. Displacement components We Versus WW) . 
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Calculations were made for a = 8 = 2 and the results plotted in Figs. 12 and 13. 
Stress-strain diagrams were not plotted since the end displacement varies with the 
particular generator chosen 
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The conclusions drawn for a panel with uniform shortening also apply to a panel 
with uniform longitudinal stress. The arrows in Fig. 12 show the points corresponding 
to zero values of w, and w, for a panel without initial deformations. Branches 7, , r2 , 
r, and r, starting from these points correspond to the branches starting from the origin 
in Fig. 13. 

Curves 7, and s, in Fig. 13 are practically coincident. The dotted curve in Fig. 12 
is obtained equating w, and f, to zero and disregarding the second and fourth of Eqs. (31). 


CONCLUSIONS 


a) The cylindrical shell under axial compression 

An approximate analysis has been made of the behaviour of an axially loaded cylindrical 
shell using a four-term expression for the stress function and also for the displacement 
function. However, numerical calculations were generally made using three terms in 

















Fic. 14. Lines of equal deflection of wave surfaces, according to eq. (10), for w; = wy = 0, w, = 3.2ue. 


the expressions for the stress functions and two in the expressions for the displacements. 
The two terms retained are of the form sin (ax/a) sin (xy/b) and cos (2rx/a), since 
previous experience has shown (2) that these terms are of greatest importance in the 
first phase of the buckling process. 

Three modes of buckling appear possible for a cylinder with perfect initial form. 
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The most unfavorable one corresponds to the lower branch OB of the corresponding 
curve in Fig. 1. In this case, the ratio of the amplitudes of the displacement components 
is about 4. The generators of the cylinder which contain the crests of the waves in the 
buckled form take a shape similar to that given by the curve (0) in Fig. 4 but with oppo- 
site sign, that is, with flat tops directed outwards from the cylinder and sharp tops 
directed inwards. 

The curves formed by the intersection of the deflected surface with coaxial cylinders 
of various radii are shown in Fig. 14;° they refer to the case when the ratio of the ampli- 
tudes is 3.2. The axis of the cylinder is vertical in Fig. 14 and it is evident that the 
outwardly directed crests are elongated in the direction of the axis. This phenomenon 
is often noticed in buckled specimens. 

The favorable type of buckling whose stress-strain characteristics are represented 
by the branch OA of the curve in Fig. 1, has generators of the form given by curve (0) 
in Fig. 4. If this type of deformation could be ensured, the behaviour of a cylindrical 
shell under axial load would be practically unaffected by buckling, and the stress-strain 
curve would continue to follow the line representing Hooke’s law. 

The third type of buckling corresponds to an axially symmetrical configuration and 
this would be represented on the stress-strain diagram in Fig. 1 as a horizontal line 
through the point 0. 

When the second deflection component in Eq. ( 
inward crests are sharper than the outward ones), under axial loading, the deformations 
rapidly approach the most unfavorable form, no matter what the magnitudes of the 
initial deformations may be. The stress-strain curves for various initial deformations 


= 
‘ 


) is initially positive (that is, the 


are given in Fig. 1. 

When the second component in Eq. (7) is initially negative, it might seem possible 
to produce a favorable mode of buckling. However, a more accurate analysis has shown 
this to be impossible. If the first component is larger than the second, the magnitude 
of the latter is further decreased by axial loading until it reverses its sign and gives rise 
to the unfavorable mode of buckling. This is indicated by the curves c and d in Fig. 2. 
If the first component is the smaller, the second component could increase under load 
and lead to axially symmetrical buckling. However, this could only occur if planes 
normal to the axis of the cylinder containing the maxima of the first component, also 
contained the maxima of the second component. In reality this perfect coincidence 
never occurs. Then, as is shown by the analysis based on Eqs. (20), a rapid relative 
shift of the two components occurs. When the position of one component has changed 
by half a wavelength with respect to the other, the deflected form is practically inverted 
as indicated by the various curves in Fig. 4. The stress-strain curve for this case is still 
unfavorable. 

It may therefore be concluded that, irrespective of the initial deformations existing 
in a cylindrical shell, buckling will always follow an unfavorable mode, with the stress- 
strain curve rapidly approaching that pertaining to the lower branch for a shell of per- 
fect initial form. With cylinders having initial deformations of wavelength about 2x +/Rt, 


it is impossible to improve the stress-strain characteristics. 


5The radial distance from the intersection to the internally tangent cylinder is given on each curve as 
a fraction of the difference in radii of the externally and internally tangent cylinders. The dotted curve is 
the locus of points where the tangent plane is parallel to the axis of the cylinder. 
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For wavelengths greater than 27+/Rt, it may be shown that, at buckling, the stress- 
strain curves start from a point for which n > 1 but that they fall rapidly to low minima. 
Such initial deformations are therefore dangerous. 

The possibility of using initial deformations to improve the stress-strain char- 
acteristics of cylinders is only feasible for wavelengths considerably smaller than 24 +/Rt. 
Such wavelengths might stiffen the shell and discourage the buckling. This possibility 
is best investigated experimentally.° 


(b) The panel under axial compression 


The compressed cylindrical panel was analyzed assuming, as to radial displacements, 
that the curved sides were clamped and that the straight sides were simply supported. 
The expression chosen for the radial displacements contains two terms and Fig. 15 














A~ x 


Fic. 15. Deformation of mid-generator, for a = 6b; (a) according to eq. (32), (b) according to eq. (45). 


shows the forms assumed by the middle generator for three ratios of the maximum 
values of these terms. 

Calculations were made for two special cases: (a) panel subjected to constant axial 
shortening and, (b) panel subjected to constant axial stress. The conclusions drawn 
from the two cases are the same. As for a complete cylindrical shell, the deformations 
tend towards the unfavorable modes but it is possible for the critical load of the linear 
theory to be exceeded before the beneficial effects of certain initial deformations are 
nullified. The approximate analysis given in the Appendix shows that favorable initial 
deformations may be obtained by initially bending a panel to a radius which is smaller 
than that at which it is tested. This may be seen from the stress-strain curves n, and 
q, in Figs. 7 and 11 respectively, which correspond to the curve w,/w, = 1/4 in Fig. 
15(a); Welter’s experimental results also indicate that a continuous and regular stress- 
strain curve is possible. With proper experimental verification, the foregoing analysis 
will serve to indicate how the behaviour of a compressed panel may be improved. 


6For short cylinders, the barrel shape might be beneficial. 
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APPENDIX 


Initial deformations caused by clamping a panel in cylindrical grips. 

In Welter’s tests (3) some panels were bent to various radii and then tested in grips 
of smaller radius. In order to estimate the initial deformations likely to have occurred 
in those conditions, we consider the following schematic problem. 


A rectangular panel whose sides are x = 0, x = 2a, y = O and y = J, is initially 


bent to a shape w* h sin ¢ and is then constrained such that the following boundary 


conditions are imposed 


ow . Ow : 
w - 0 for z = 0, 2a: w= —; = 0 for y = 0, D. 
Ox oy” ‘ 


Making use of Eq. (5) and assuming small displacements, we write 
Vi(w — w*) = O. 
According to the boundary conditions, the solution takes the form: 


hsin ¢ 
* (sinh xp + xo cosh x_)(cosh x» — cosh x) 
Xo 


sinh xo cosh x, 4 
+ (x sinh xy — xo sinh xo) sinh xo], 
where 

x = r(zr — a) ‘b, Xo = ra b. 


This equation gives an estimate of the initial deformations occurring when a panel, 


initially bent to a cylindrical surface of small curvature, is clamped in grips, also of 
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small curvature, but with a camber different by / from that of the panel. The deflections 
obtained for a b, for the generator y = 6/2 are plotted in Fig. 15(b). 
The condition encountered in Welter’s tests were, of course, somewhat different 


from the above. lor example, the grips in the tests gave a large camber to the panels 
and this would require the use of Eqs. (5) in their complete form. Rotation of the straight 
edges was also restricted by the toothed guides. For these and other reasons, it is obvious 


that the above is not a rigorous solution. However it will give an estimate of the initial 


deformations which may be expected. 
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ON THE PLANE PROBLEM OF A PERFECT PLASTIC BODY’ 


BY 
HILDA GEIRINGER 
Wheaton College 


In the equilibrium problem of the isotropic ideal plane plastic body it is assumed 
that the state of stress is restricted to a one-dimensional variety F(c, , 2) = 0 where 
a 1, 2) denote the principal stresses; moreover, the two equilibrium conditions hold. 
This problem may be derived in several ways from the complete three-dimensional 
problem of the perfect plastic body; in particular the assumption a; = ¢, = 0,0 ---/dz = 
0 leads to it without the need for consideration of any three-dimensional relation be- 
tween strain-rates and stresses. However, as long as F(o, , o2) is not specified the 
relation of the plane problem to the three-dimensional problem need not be considered. 

Important particular assumptions for the yield condition, F(c,; , a2) = 0, are the 
v. Mises condition of constant shear energy and the Saint Venant condition of constant 
maximum shearing stress. For the first, F is given by 


2 2 72 
Oo; — 0102 + o, = 4K’, (1) 
the second states that | 7 |max lo, — o,|/2 = K, if o, and o, have opposite signs, 
while otherwise either o, = 2K, or o. = 2K. Hence 
| G; ~ G2 | = 2x. if 01:02 <= 0, 
(2) 


4K — |o,+ 02], if o,\0, = 0. 


These forms of F are derived from corresponding three-dimensional conditions by 
means of the afore-mentioned assumption of “plane stress”, ¢, = 0, 0d ---/dz = 0. 

The problem most widely investigated is that with the yield condition (¢, — o2)? 
4K’ as in the first Eq. (2). It is distinguished by its comparative mathematical pra Aas 
moreover, this same yield condition appears in “plane strain”, as both v. Mises’ and 
Saint Venant’s condition. 

The general plane problem defined by the two equilibrium equations and some yield 
condition has been considered by Sokolovsky [6a, b] for the two conditions (1) and (2); 
and for general yield conditions by v. Mises [4a], Geiringer [1], and the Brown Uni- 
versity group, (see P. G. Hodge [2], Chap. VIII). As in the case of “plane strain” one 
can—in several ways—derive from the basic equations a linear system by a change of 
variables, somehow analogous to the Chaplygin transformation in the theory of com- 
pressible fluids; the plane of the new independent variables called “‘stress graph” by 
v. Mises, corresponds then to the hodograph plane, the yield condition to the “adiabatic’ 
condition. It is intended to follow up this analogy in several directions. 

First we establish linear differential equations of second order for certain magnitudes 
which may serve as coordinates in the physical plane; such an equation is amenable to 
well-known analytic methods. The coefficients depend, of course, on the yield condition; 


1Received Oct. 13, 1950. The main results of this paper have been communicated in a talk given at 
the Second Symposium on Plasticity, Providence, April 1949, and later briefly summarized in three 
abstracts of papers presented by title at the October 1949 meeting of the A.M.S. [1]. 
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we consider particularly conditions (1) and (2) and the new “parabola condition”’ 
by v. Mises. Next the characteristics are investigated, in general, 


recently proposed ta 
here most of our results while definitely 


and for the afore-m«¢ ntioned yield conditions; 
1 


more simple are essentially identical with corresponding results of Sokolovsky [6a], 


Mise S ta], Hodge Z Finally, sim ple wave solutions are considered. The stress dis- 


v. 
tributions are determined and discussed; cross characteristics and lines of principal stress 
are computed. In each of the indicated directions one can go much further, this paper 
thus being of an introductory character; solutions of concrete problems have not been 
included. 

1. The stress graph. Following v. Mises we express the plasticity condition in terms 


of an appropriate parameter, 8: 
J ii ts) = U, J a,(S), 0 = 0; 8). 


Let us assume from now on that o, = o, . It is often convenient to take for s the average 


stress (o, + o,)/ 2K. Denoting by @ the angle between the positive x direction and the 


£ direction, i.e., the (positive) first principal direction, we take s and @ as new variables 
which determine the stress tensor at each point. The plane in which s and @ are ,polar 
coordinates is called ess graph by v. Mises. Denote by 0 ---/0& and @ dn the 
directional derivatives in the first and second principal direction, respectively, then 
equilibrium conditions are equivalent to the system (3): 
Oc 06 Oo» 00 
= (f— “ 6:1) => ~ (ee Cas Ge (4 
EY 2 an On dg 
With the use of (3 ind the abbreviations. da,/ds = ot S$), ete., 
ids : ” 
tT, — G;)/0, = JI (g2 — 6))/a2 = Qs), (2) 
(4) yields the reducible equations 
Os ‘ 06 Os 00 ; 
— = f(s) —, — = 9(8) ; (6) 
d& On On 0& 


At every point where the Jacobian j = 0(s, @)/d(é, ) is different from zero in (6) the 
dependent and independent variables may be interchanged, to yield v. Mises’ linear 
equations 

On 

—s ( 


Os 


~] 


0 iui 0g 
: m f(s) —, — 
00 . Os 06 


= g(s) 


2. Linear differential equations of second order. In Eq. (6) & and 7 are not coordi- 
nates; it is, however, possible to derive from (7) in various ways systems of equations for 


quantities which may serve as coordinates in the physical plane. If, e.g., 
X xcos 6+ y sin 8, Y = ycos 6 — xsin 8, (8) 
we find the system of linear equations of first order: 


ox = ¢(s) oY ron ay os fie aX me (9) 


og Os 00 
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From (9) we obtain the equations of second order: 


ax a7x = 
a f9 “—_—* —-X + 3s (fgt+f-9), 
(10) 
a°Y a°Y _— = 
aes += (’f+f—g). 


These equations’ are hyperbolic, or elliptic, according to whether fg 2 0 (see (5)) or: 


ats We > 9; (11) 
they are parabolic if either fg or its reciprocal value equals zero. If a solution X(s, 6) 
of (10') has been found, the corresponding Y(s, @) follows from (9°) and the coordinates 
x, y follow from (8) in terms of s and @ without further integration; thus s and @ are 
defined in terms of x and y. 
Next, we satisfy the first Eq. (9) by a function ¥(s, 6) such that 


— an fw se 
ds __g(s)k(s)’ 060 ——sik(s)’ (12) 
where k satisfies the condition g(s)k’(s) + k(s) = 0. (12’) 
Substituting y into the second Eq. (9), we find 
OY FW _ Wee ge 
a6" fg ag de (fg f g)- (13) 
Also with 
ns ol FF a f(gh'(s) — his) = 0, (14) 


ds f(s)h(s)’ a0 h(s)’ 
the equation 


We _ - Fe _ Wey 
99 9 of = 95 OF Fath (15) 


follows in the same manner. 

For finding integrals of (10) or (13), or (15) two principal methods are available: 
(a) Bergman’s operator method may be adapted to our problem; (b) the classical way 
of expansion in series seems promising. In these ways regular solutions as well as solu- 
tions with appropriate singularities may be obtained. Because of the linearity of (10) 
solutions may be superimposed. On the other hand, the difficulties encountered in actual 
boundary value problems, where the boundary conditions are given in the physical 
plane, are well known. 

3. Particular yield conditions. We consider v. Mises’ “quadratic” condition (1), the 
the “hexagonal” condition of Saint Venant, and v. Mises’ new “parabola condition’’. 

1Equations (10) are simpler than Eq. (19) of v. Mises [4a] because X and Y rather than z and y have 


been chosen as dependent variables. 
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In the o, — o2-plane (Fig. 1) the Mises limit (1) is represented by the ellipse 
s= i (o1 + Go) (-2 sss +2) 
4 IK "1 2/3 \ eo=S— 2? = @}) 


=~ STS Bain OTS 
o, = K{ s oh (3 _*) | 3 = K|s a (454) | (16) 


where the upper signs apply to the right, the lower to the left part of the ellipse. Then 
a 4K?3-—s (¢2 — 0;)° (4 — 3’)? a, 
= = me Tlinngyn am . (17) 


0102 >— 2 
3 4-8” 





fe —= , o . 2 
dg O20; s-s 

It is seen that fg > O for s’ < 3, fg < Ofors’® > 3 while fg — +@ as s* — 3. Since 
we assume o, = o; we are concerned merely with the left part of the ellipse; the hyper- 



































bolic case corresponds to the interval — 7/3 < s < + V/3;tos = + V/3 ands = —+~/3 
6,22K ad 
ai Pa ? 3 | fi 
| ~ / 
! i 
—_ Pp ! oP 
s=- oe 7 / 
- 7 
- 4 
9,2 2K al Jey? / i’ 
a 4 12? XY C—— 
s=-2 s=-a 


Fic. 1. Yield conditions. 


correspond the “sonic points” 0, = 20, = 4K/+/3 and o, = 20, = —4K/+/3. In 
the physical plane the “sonic line’ separates the domain of hyperbolic solutions from 
the elliptic domain; its image in the stress graph is the circle about the origin with 
radius s = +/3. The representation of the “hyperbolic part’’, for a. > «, is: 


o, = K|s _ (4 —*) ‘| o, = Ks + (t=) ’ ls|< V3. 16’) 


The complete Saint-Venant limit (2) is represented by a hexagon inscribed in the 
ellipse (16). With s = (o,; + o.)/2K, the representation of the “hyperbolic part” of 
this limit is 

o, = K(s — 1), o, = K(s + 1), je] <1 (18) 
and 
oe, = K’, fg = 4. (19) 
The vertical and horizontal lines of the hexagon correspond to parabolic domains, where 
either o; or oc» are constant. 

The ‘parabola limit’? consists of two branches of parabolas passing through four 

corners of the Saint Venant hexagon: 


ok. a? nts)" ’ 
K os x | ( 2K P (20) 
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where a = 1 + +/2, the plus sign applying to the left, the minus to the right branch. 
The parametric representation for ¢. > o, is (with our parameter) 


a ae oe 
o, = Ks 9, (@ s°), = s+ oi (a s*), ls|<a (21) 


and 
f=a-s, g=at+s, fg=a-—s’. (22) 
It is seen that, with the only exception of the points s = --a, where fg = 0, the problem 
is everywhere hyperbolic. 
As an application let us consider the second order equations of the preceding section 
for the parabola limit, transforming them at the same time to standard forms. With 


“ = (a — s’)”, 8 = asinu, u = arcsin s/a 
and putting X(s, @) = X(u, 6), and A(u) = (1 + 2sin u)/cos u, Eq. (10) becomes 
ax @xX, = ~ 
FY ue + X-A(u) = 0, | u | < 2° 








Next with p(u) = [cos u(1 — sin u)]"’”, where Ap — 2p’ = 0, it assumes the form 


Z #Z, ,1+2sinu 


00 aw “ 4 cos’ u 





= 0. (10a) 


This is an equation in standard form which may be treated in various manners. 
Also, with ¥(s, 6) = W(u, 0), Eq. (13) takes the form 


vy _ ay, ay 1 


a haw" (18a) 


? 
and putting yi u, 0) = r(u)R(u, 6) where 


1 + sin u\'” 
r(u) = (1+ simu , 2r’ cosu — r = 0, 
cos u 


we obtain 


eR oR 1— 2sinu _ 
00 +2 4cosu ©. (186) 


Note the difference between (10a) and (13b). 
As a second application, we give an example of particular solutions of Eqs. (13). Set 


¥(s, 6) = P(s)Q(0). (23) 
One obtains in the usual way, denoting by \ a constant, and using (22): 
(a? — s°)P” — (a+ s)P’ — »\P = 0. 
Introducing the new variables u and U by 


v-au-p, @-1 @P_ iu 
a , ds 2a du’ ds’ 4a” du” 
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we find 


PU dU _ 


u(l — w) 7 u U = @. (24 
) du du ) 

which is a hypergeometric equation. Solutions of (24) are well known; with Q = C sin 

n(@ — 6), we may thus find solutions of the form (23), which may be superimposed. 


4. Characteristics. The characteristics of the reducible equations (6) are easily found. 
If ¢ denotes the angle of a characteristic direction with the &direction, we find by standard 
methods: 

: f dy f\'” 
tang = ~-, tang =>; =x : 
g dé 9 


In the hyperbolic domain, f and g have the same sign; thus, at each point of this 
domain, there are two distinct directions making equal angles with the édirection. We 














—JX 0! 








Fic. 2. Characteristics in physical plane and in stress plane. 


denote by C*, (C_) the characteristics making the angles +¢, (—¢) with the &direction, 
and by 
a=6@+4¢, a’ =9O-¢g (26) 


the angles of C” and C” with the z-axis (Fig. 2). Then, with dy/dz for the direction 
coefficient of a characteristic, we have 


1/2 


di tan 6 + (f/g) ans 
Y = ten (0+ ¢) = “4, (26’) 
dx 1] += tan @-(f/g) 

the upper (lower) sign applying to C* (to C’). If a solution s = s(zy), 0 = O(x, y) is 


introduced into (26’), this becomes an ordinary differential equation for the determination 


of the characteristic curves. 


Note also that, according to (3), (5), and (25) with F; = dF /de; , 
: ; 
9 Oo do» F, -; 
tan’ 9 = 9G = SO = > (25 
ea, day F, } 


depending merely on the yield condition: F(c; , 2) = 0. For the “quadratic’’ condition, 
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for example, this gives tan” ¢ = (a0. — 20,)/(20, — o,), showing that along the “sonic 
line’ ¢ equals 0 or 7/2. 

From the general definition of characteristics it follows that a relation exists between 
the derivatives of the dependent variables along a characteristic direction. To find this, 
. a ° . “ . 1/2 _¢ ° 
multiply the first term in the first Eq. (6) by cos g, the second by (g/f)'’~ sin ¢; in the 

‘ , : . . . 2 1 
second Eq. (6) multiply the first term by sin ¢g, the second by (f/g)'” cos ¢ and add. 
Thus, 

Of . : , (00 06 . 
cosg + =—sNng = (f9)' sf cos g + = sin ¢]. 

On 0& On 
Denoting by @ ---/dl differentiation in a characteristic direction, we obtain the desired 
relations: 

2 06 


—(f9) , (27) 


Os , 1/2 00 Os 
—— i ——- 
al al al 
These relations between the derivatives of s and @ along C* or C” do not depend on the 
validity of the transformation (7). 

It follows from equation (10) or (15) that for the fixed characteristics Tin the stress 


graph: 


‘de\* do \ 1/2 
( ) =(f)"', z= +(f). (28) 
ds ds 
Writing: 
[ (ig? ds = Fos) (28’) 
with an adequate lower limit), we obtain the equations of the two families of fixed 
characteristics, T° and T 
@ — F(s) = const, 6+ F(s) = const. (28’’) 


Equation (28) may also be derived from (7) in the same way as (24) was derived from 
(6). Then (27) shows that (as always) the characteristics in the two planes which are 
mapped onto each other (here the x, y plane and the s, 6 plane) correspond to each 
other, here C* to [°, C” to T-. We see from (28) that [ and I make equal angles 
+y and —y with the radius vector and tan y = sd6/ds (Fig. 2). 

We consider here in general the hyperbolic case; let us, however, just say a word 
about the parabolic case where f/g = o3/a{ is either zero or infinite. If 2 = O(o; = 0) 
the angle ¢ = 0, (¢ = 7/2) and the one existing characteristic direction coincides with 
the direction of first (second) principal stress. In either case y = 0, d@/ds = 0. It follows 
from (28) that in a parabolic region the characteristics T in the stress graph are radii 
through O’, and from (26’) by integration that the characteristics in the physical plane 
form likewise a set of straight lines: y = x tan 6 + G(@) if of = 0. 

5. Characteristics. Continuation. The I’-curves in the stress graph have been investi- 
gated by v. Mises for the three yield conditions and sketched with s and @ as polar 
coordinates. The formulas are: 


1In the hyperbolic case g/f is nowhere zero or infinite. 
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Quadratic condition (v. Mises) 


do_ (3 -8)" 





ds ~~ 4-8 
(29) 
rc oy /o 1 Ss om 
+ 6 = arctan [s(3 — s*)"'”*] — 5 arctan E (3 — 8’) | + const, 
and for the C- characteristics: 
(12 — 3¢°)'*7 — 
tang = +— a» Nin (30) 


2(3 — s*)'” 


At the “‘sonic points” tan ¢ is zero or infinite and y = 0. The characteristics (29) are 
congruent curves (like the epicycloids in our analogy). In rectangular coordinates, s and 8, 
a I’ line is monotonously increasing e.g. from — 7/4 to 1/4, showing central symmetry, 
inflection points at s = 0,s = + ~/2, and dé/ds = Oats = + V3 (Fig. 3). 


8 
4 


oS 
S 








Fic. 3. 6 = F(s). 


Saint Venant condition. Here d@/ds = +1/2. The I-characteristics in the hyperbolic 
area in the stress graph are ordinary spirals +6 = s/2 + const, or straight lines in the 
rectangular coordinates s, @ (in the parabolic area the unique set consists of radii through 
0’); also tan g = +1,a = 6+ 7/4, a’ = 6 — 2/4; it is thus seen that the characteristics 
in the zy-plane coincide with the shear lines, a result well known in the classical theory. 


Parabola condition. 


dé 2 2, -1/2 - § 
— = +(q? — 3”), @ = -tarcsin — + const. (31) 
ds a 
These are sine-lines s = -ta sin (@ — 4) if interpreted in rectangular coordinates and 
circles through the origin in polar coordinates. Here, 


1/2 
° oll 
tan g = +( ) ‘ 32 
¢g — 9" (32) 

We finish this section by a remark on the approximate solution of the classical initial 
value problems in the hyperbolic domain. In the Cauchy problem, on a smooth arc of 
curve, K, values of s and @ are prescribed (with continuous second derivatives) so that 
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K has nowhere a characteristic direction, i.e. that (26’) does not hold. These data assure 
uniqueness of solution in the characteristic quadrangle determined by both, C* and C’, 
at both endpoints of K. There exists certainly a solution in the neighborhood of K which on 
K takes on the prescribed values. It can be obtained approximately by using the fixed 
net of I'-characteristics (adjacent to the image K’ of K) in order to find the (s, 6) values 
for corresponding lattice points in the physical plane. In general, we may then be able 
to continue in this way until the whole characteristic quadrangle has been filled. The 
procedure may be refined in various ways. 

For the characteristic initial value problem of reducible equations, like (6), the data 
which can be arbitrarily prescribed are not the same as in the better known linear case. 
(This difference is often not appreciated, see [2], p. 197.) In the reducible case any arbi- 
trarily chosen curve C : y = y(x) may be a characteristic, but along this C neither s 
nor 6 may be arbitrarily prescribed, since along C one relation (26’) and one (28) (or 
rather (27)) must hold. More generally, for a curve y = y(x) to be a characteristic, the 
three functions y(x), s(x), 6(x) must satisfy two relations: one (26’), and one (28). 

For two intersecting characteristics, C°, C” the values of s and @ at the point of inter- 
section must satisfy both relations (26’); the relations (28), which merely connect dé 


and ds leave the values at one point unrestricted (or, in other words, the constants in 


Envelope 


Fig. 4. Family of straight characteristics. 


(28’’) may be determined so that the four equations (26’) and (28”’) are satisfied at the 
point of intersection). If in this way two intersecting characteristics are given, the 
solution in the corresponding quadrangle is uniquely determined, and it may be found 
step by step in a way analogous to that indicated above. 

6. Simple waves. A simple wave region is a domain D in the physical plane mapped 
onto one single characteristics, say I , of the s, @ plane such that to each point on 
[> corresponds a straight line in D carrying constant values of s and @ (and consequently 
of o; , o> , 0). In such a “forward wave” the straight lines are the C” characteristics, and 
@ — F(s) is constant for all (2, y) in D; all characteristics of the other family, the “cross 
characteristics”, C~, are mapped onto Ty . The Jacobian, 7, of this mapping vanishes, 
of course. Similarly we define a “backward wave’. Naturally, not all the above noted 
properties are needed for the definition of a simple wave. The notion includes as particular 
case that of constant stress tensor in D. 

Here we start with the problem of finding solutions of Eqs. (6) such that the state of 
stress is constant along each straight line of an arbitrary set of lines. Denote by 6 the angle 
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of such a straight line, Z, with the 2 direction (Fig. 4). Since o; , «2 , and 6 remain constant 


along L, (4) yields 


do, ; dé 
sin (0 — 6 = (¢. — o;) cos (@ — 6) — 
In ( ) di (G2 a,) cos (0 6) da’ 

(33) 
do» ; dé 
cos (6 — 6) —2 = (6. — 8 — $)—. 
Os (6 0) dé (G2 o;) sin (6 0) dé 

By division, we obtain 
' do» — F,(s) 
tan” (@ — 6) = = ——— 2. 
‘ do, F,(s) (34) 


as in (25). Next we multiply the first Eq. (33) by cos (@ — 4), the second by sin (@ — 4) 
and add: 


d(o, + a2) i a (35) 
a : igegk -<auem 00 


dé ~ gin (6 — 6) cos (6 — 4) 


From (35) and (34) we then find easily our Eq. (28). In fact, 


‘ d(o Talice oy x ; , (do.do,)'”” 
dé =- 1 r “Gite —- ¢€ cete— @ = ta 
o> — 0; ts. > 
dé (aia) een a 
oe oe. (35’) 
ds Os ~ GO; 


Thus, if a condition (3) and a I’ are chosen we find to each 6 the s and @ by (28) and 
(34). If, however, we want to find cross characteristics, C*, and lines of principal stress 
we have to consider not only 6 but the whole set of straight characteristics. Let it be 
given in the form 

y = Bx + his), (36) 


where 6 = tan 6 tan (@ — g). The differential equation of the cross characteristics 
is dy/dx = tan (@ + ¢) and since @ and ¢ both depend in a given way on s, and s on 6 


we have 


dy = tan (6+ ¢) = k(6). (37 
dx 


From (36) we derive dy = x dB + B dx + h’d B. Substituting this in (37), we obtain 
dx My h’(B) 
oir (38) 
dg Ks) — B 


Thus by means of quadratures x is obtained as a function of 8; this together with (36) 
provides a parametric representation of the cross characteristics. For the principal 
stress-lines the procedure is the same except that dy/dx = tan 6 is used instead of (37) 
while the relation between s and 8 is the same as before. In case of a centered wave 
h(p) = 0. 


It is well known that simple waves are adjacent to regions of constant state; this is 
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deduced immediately from the fact that a region of constant s, 6 is mappeed onto a 
single point in the stress graph. 

A simple wave can transform any state s = s, , @ = @, into another state s = 8, , 
6 = 6, provided either 6 + F(s) or @ — F(s) (where F(s) is defined in (28’)) has the same 
value for the first and second state. By combining a forward and a backward wave and 
inserting a uniform state in between, any s, , 6. can be reached; in each wave the envelope 
of the straight characteristics can still be chosen in various ways. 

7. Examples of simple waves. Consider the quadratic yield condition (1) and the 
forward wave, image of 6 = F(s). Using (28’) and (29): 


; — l 8 — 
6 = Fs arctan [s(3 — s°)"'”*] — 5 aretan E (3 — s°) a (39) 
Here the C are straight lines, 6 = 6 — ¢, and as a, > o, , tan ¢ is given by (30) with 
+- sign. Hence 
8 8 (> Read 38°)'? 3% 


] 
6 6 — ~ = arctan = 7s — = arctan — a. = an -~—— 3 
. (3-—s)'” 2 2(3 — s°)'” 2(3 — s°)’ 





To simplify, we introduce ¢ = +s(3 — s*)~'’” and obtain 


] t 1 9 > 
6 arctan ¢ — 5 arctan 57 arctan E (4+ °)'? - |, 


which reduces to 
6 = arctan ¢t — > t = tan (s + r) (40) 
or, reintroducing 
s= V3sin (s + 4] 6 = F(s). (41) 


Here as s goes from — 7/3, to 0, to + +/3, (¢ goes from — ~, to 0, to +), 6 goes from 

7/4, to 0, to +7/4 and 6 from —32/4 to +7/4 (see Fig. 3). Hence in a “complete” 
wave, 6 varies by 180°. In actual problems a solution is given in such parts of the x, y 
plane only where the straight characteristics do not intersect; thus parts of complete 
waves may appear as solutions. 

The image of the I, characteristic 6 + F(s) = 0 is a backward wave with rectilinear 
C* lines and 6 = 6+ ¢ = —[F(s) — g]. Hence this I” wave is the reflected [~ wave 
considered above and 6 = 2/4 — arctan /, etc. We see that in a foward wave as 6 in- 
creases, the mean pressure as well as @ increase in a monotonous way. In a backward 


wave the opposite is true. 


Next, we compute the cross characteristics for the backward wave @ = —F(s), 
6 r/4 — arctan ¢. With notations as in (26), (36), and (37), we obtain: 
3 t 5 tan (* retan ¢) Scat l Sal 
8 = tan 6 = tan |— — arcta 5 ik rer 
4 1+t 1+8 





> 


_f+t+32 B+B+2 _ k(8) 
?—t+2 2° +e6@+1 °°" 


tan (6 — ¢) = tana’ 
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We may then continue as explained above. For the particular case of a centered wave, 
where h(8) = 0, 8 = y/z, the differential equation of the cross characteristics becomes 


2 a 9,7 
SL 2222 (42) 
dx 2y +ay+ 2x 


tana’ = 
which can be integrated easily. One finds, using polar coordinates, rs (Fig. 5) 
: ar P { P 
r” cos (; - i) = const. (43) 


Let us conclude this consideration of the quadratic condition by deriving a simple 
characterization of the stress tensor which holds if at each point of the wave we use a 


y 
A 











Fria. 5. Quadratic limit. Complete centered wave. Straight characteristics and cross characteristics. .~ 


coordinate system x’, y’ where z’ is the direction of the straight line through this point, 
y’ the perpendicular direction. The stress tensor is then given by o,- , a, , Tz", . Denote 


by »; , v2 the two principal directions at a point and let o,, = o, o,, = 6, Tey" = T} 
then (vy,;) = 6 — 6 = ¢g and we find 
o = 9; cos ¢ + o2 Sin’ ¢, tT = (a2 — o,) SiN ¢g COS ¢. (44) 


Using primes to denote derivatives, we find from (1), (16’), (34): 


cos’ ¢ = a//2K, sin* g = o;/2K, 
1 | , , 1 d 2 4 3 d 
= OK (o10; + 0202) = 4K ds (0; + o>) — 4K ds (0,02) 


(45) 


— ey d 4K* 9? va 2K = 2K Te 
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Also, 


1 1d 
o = OK (0103 + o201) = OK ds (7192) a Qe. (46) 


Hence under quadratic yield condition the following simple relations hold for the stress 
tensor all over in a simple wave region 


o+7=-=K’, a = 2. (47) 
Parabola limit. Using (31) and s/a = t, we consider the image of the [™ characteristic 


6 = —arcsin ¢ = arccost — . (48) 














Fic. 6a. Parabola limit. Straight characteristics and cross characteristics. 
Fig. 6b. Parabola limit. Lines of principal stress. 


From 


St A" se (* 8) 
ten y = (2=4) "Vtg Ue (49) 


we conclude that 





and thus: 





= 2(5-1), s = a cos (0 + 5) = a cos 3. (50) 
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Here the maximum interval for @ equals 180° and for 6 it equals 270°; as 6 goes from 
1/2 through 0 to —7/2, 6 goes from x through 0 to —7/2. 

We compute now the cross characteristics for a centered wave. From (50), using polar 
coordinates 7 and 6 we obtain (Tig. 6a) 


dé 2 : 
j = —tan 2p = —tan (z + a) = —tan - et (51) 
adi 


. 26+ n\°" - 
r = 7r,.\sin ——— ; (52) 


We finally compute for the same example the principal stress lines. For the é-lines, 
the lines that make the angle 6 with the z-direction, we have, in polar coordinates 
(Fig. 6b) 


r dé 


. ‘ 6 T ae 
= —tan(o-— 6) = —tan (: + , (53) 
dr 3 


26+7|° 
r = To) sin ——.— (54) 
) 


96+ ~ 
f = 1,| cos ——— : (55) 
6 
These examples may suffice. 


Simple waves are widely used in the solution of actual problems. 


or 


and for the 7-lines 
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—NOTES— 
NOTE ON THE TORSIONAL RIGIDITY OF SEMI-CIRCULAR BARS* 
By P. L. SHENG (Taipeh, Taiwan) 


In such cases, the general stress function for torsion for the cross-section of a circular 


sectol reduces to 


" 6a” = woul t) 008 
eu = | -r0 + cos 26) +9 YY (-1) na(t) sent | 


2 yn a] n(n® — 4) 


where r, 6 are the polar coordinates, and a the radius of the semi-circle. It is easy to 
verify that Y = 0 at the boundary. 
The applied torque 7’ is given by the well-known formula 


T = 2Ga [ vas = 2Ga [| Wr dr dé, 
“A VA 


where f, or ff, are taken over the cross-sectional area, G denotes the shearing modulus 
of elasticity, and a the angle of twist per unit length of the bar. 
Hence 


w/2 a 
T= Gol [ [ r'(1 + cos 26) dr dé 


¥—-#/2 “0 


Neat DNF HTH) naa] 
+ > . f | = cos n6 dr dé 


2 
T antas. 00 MN — & Jes Jo 


T 


~ Gol Zot 3" oF 1 | 
~ Gal ‘* T mm ie (n — 2)n*(n + 2)* J° 


Using the well-known formula 





we can easily derive by the split 





— sa -4 (etx 5 4 2 | 
(n — 2)n?(n +2)? 32L(n—2)n ' n(n+2) nn” (n+ 2)’ 


that 
= l 1 ss 
p ————_——— 3 = s|/1-sF I. 
nmi ge.ees (Nn — 2)n'(n + 2) 8 16 
*Received September 30, 1950. 


1Cf., e.g., S. Timoshenko’s Theory of elasticity, §79, noticing that in Timoshenko’s book the stress 
function ¢ is¢ = Ga, and 77, a, and (r, 6) in the present paper are written as M,, 0, and (r, y) respec- 
tively. 
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T = Gaa'| ® = 4] 
2 ef 


exactly, and the numerical expression is 


By substitution, 


T = 0.297,556,782 Gaa* 
which differ slightly from St. Venant’s result of 0.296 or 0.2966.” 


*The value 0.296 is from the table in 8. Timoshenko’s Theory of elasticity, p. 250, 1st ed., 2nd and 8th 
impression (1933 and 1934), McGraw-Hill; while from the datum given in I. Todhunter and K. Pearson’s 
A history of the theory of elasticity and of the strength of materials, vol. II, part I, p. 193, we have 


T = 0.3776M = 0.3776 X Ga X (xa*)/2 XK a?/2 = 0.2966Gaa‘. 


A NOTE ON MY PAPER 


ON STEADY LAMINAR TWO-DIMENSIONAL JETS IN COMPRESSIBLE 
VISCOUS GASES FAR BEHIND THE SLIT* 


QUARTERLY OF APPLIED MATHEMATICS, 7, 313-323 (1949) 
By M. Z. KRZYWOBLOCKI (University of Illinois) 


Determination of the constant of integration for the temperature distribution (p. 
317, eq. 25) from the condition that that total flux of enthalpy across jet is alike at all 
cross-sections restricts solution to small Mach numbers (if the comparison cross-section 
is close to the slit, as pointed out by A. H. Shapiro, Appl. Mech. Rev. IIT (1950) p. 415, 
No. 2718) or to high Mach numbers (if the comparison section is far from the slit). To 
take into account all the relative cases, that constant may be determined from the 
condition that the total flux (enthalpy plus kinetic energy) is alike at all cross-sections: 


4 


2 | (Jc,T, + uj/2)u;(po + pi) dy = const. 





*Received Feb. 22, 1951. 


ON THE LEAST EIGENVALUE OF HILL’S EQUATION* 
By C. R. PUTNAM (The Institute for Advanced Study, Princeton) 
The differential equation 
ze’ ++ fil<e = 0, (1) 


in which \ is a real parameter and f(t), for —» < ¢t <, is a real-valued, continuous, 
periodic function (#0), arises in problems dealing with the propagation of waves in 


*Received November 10, 1950. 
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periodic media; cf., e.g., [1], [5], [6]. For a proper choice of units on the ¢-axis it may be 
supposed that f(¢) has period 1 and hence possesses a Fourier series 


fO~ Dae", (-. = &). (2) 
There exists a sequence of finite intervals (the intervals of stability)J, : A, S A S 3 
where A, < A° < A,,, and k = 1, 2, --+ , such that (1) possesses, or fails to possess, a 
solution x (#0) which is bounded on — © < ¢ <@ according as \ does, or does not, 
belong to the (closed) set S = >> J, ; ef. [6], p. 16. It is known ({8], [3]) that the set S 
is identical with the invariant spectrum (Weyl [7], p. 251) associated with the differ- 
ential equation (1), on either the full line — © < ¢ <o@ ora half-line, say,0O St <@. 
Let \, be denoted by yu so that u is the least point of the set S. Alternatively, u can 
be defined by the requirement that (1) be oscillatory, that is, every solution of (1) 
should possess an infinity of zeros on 0 S t <© (or equivalently, in the present case, 
on —o < ¢<o) whenever) > ug, and non-oscillatory whenever \ < u. (This twofold 
characterization of u holds also in the general case of (1) in which f need not be periodic; 
cf. [2].) 
Wintner ([9], p. 116, and [10]) has obtained the following estimates of u in terms of 
the Fourier coefficients of the function f(t) defined by (2): 
—o —2>  |c, Psu —m. (3) 


| 
n=1 


The present note will be devoted to the problem of obtaining other such estimates. A 
formula for » involving the Fourier coefficients of the given function f(t) and those of 
arbitrary periodic functions, subject to certain specified conditions, will be deduced 
(cf. (12) below); furthermore, as a corollary of this formula, upper bounds for y, involving 
not only cy) , as in (3), but also arbitrary coefficients c, , will be obtained. Specifically, 
it will be shown that 

L = nN’ = Ge + Rey), N = :. 2, eee. (4) 
where §2(c) denotes the real part of a complex number c. In case f(t) satisfies 


—f) =fti+t+o), 
for some real number ¢, as, e.g., is the case if f = sin 2zt, it is clear from the properties 
of » that (4) can be refined to 
pS WN — & — | Rew) |, N = 1,2, °--. 

2. It follows from [4] (ef. the Remark on p. 636) that » satisfies 

uw = lim {elb. | (x’? — fx’) a/ | x ar}}, (5) 

Too \ JT T 

where x(t) belongs to the class of functions, 2; , which, on the half-line T S$ t <o, 
are real-valued, continuous, and have piecewise continuous first derivatives (with 
respect to any finite subinterval of 7 < t <@) and, furthermore, satisfy 


(T)=0, O<f #dt<o, | x?dt<o. 6) 
J 


“7 
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Since f(t) is periodic, it is clear that the expression {---} occurring in (5) is independent 
of 7. Hence (5) remains valid if the limit sign is removed and the expression [-- -] of 
(5) is evaluated only for functions of class Q, . 

Let x(t) denote any function of class 2, . Clearly, it is possible to define a function 


y(t) in Q , such that y(¢) = 0 for sufficiently large t, and, in addition, is such that the 
expression [---] of (5) (for 7’ = 0) evaluated for y differs from the corresponding ex- 


pression for x by less than an arbitrarily preassigned positive number. (In fact, since 
fo x dt <~, there exists a sequence of points ¢, such that t, —~© and xz(t,) — 0, as 
n—o., Let z,(t), where n = 1, 2,--- andt, < ¢ S$ t, + 1, be a sequence of continuous 
functions, with piecewise continuous first derivatives, such that z,(t,) = x(t,), Z,(t, + 1) = 
0, and 

2’ dt > 0, | 2? dt o, as t, ©. 


in Utn 


Ify,=xonO0Ositst,,y, =z, ont, Stst,+1landy, =Ofort,+1<t<o, 
it is clear that y, satisfies the conditions claimed for the function y above, provided ¢, is 
sufficiently large.) The corners of y can be “smoothed out” so that, in addition to the 
properties required above, y has a continuous first derivative on 0 S t <o. 


It follows from the above discussion that u can be defined by 


nQ »Q 
g.l.b. | (x — fax’) a/ x a], (7) 


where x belongs to Ig , the set of (real-valued) functions x(t) (40) which possess con- 
tinuous first derivatives on 0 S ¢t S Q and satisfy 2(0) = 2(Q) = 0, and Q is, for con- 
venience, an arbitrary (variable) positive integer. Any function x of Tg , and its de- 
rivative x’, can be uniformly approximated on 0 S$ t S Q by a sequence of trigonometric 
polynomials of Ig , and their derivatives, respectively. Hence, u can be defined by 


(7) where, now, x(t) is any function of the type 
x(t) = i ee ti., = G,, Zz a, = 0, (8) 


and where N, P and Q denote arbitrary positive integers. The remainder of this section 
will be devoted to obtaining a transcription of (7) in terms of the Fourier coefficients 
of the given f, defined by (2), and the Fourier coefficients of the variable function z, 
defined by (8). 

Let P and Q denote relatively prime positive integers and note that (2) and the 


first relation of (8) can be rewritten as 

, Y 2rint/¢ 2rint/Q 

Jj ~ > Ca si and x(t) = > Af =; (9) 
where 

¢. = Co, and a. = Ae forn = 0, +1, ---, +N. (10) 

(By the uniqueness theorem for Fourier series, C, = 0 if Q does not divide n and A, = 0 
if P does not divide n.) Since x? = >> B,e*"'"’°, where B, = >, A; A,-x , it follows 
from the Parseval relation that 


ee 


fea=-ed 147, | 2 a=4r/Q Xn] A.0 


70 
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-Q 


| fx dit = Q z= C.B, ‘ 


“0 


Thus, the expression [---] of (7) is equal to 


[ ix?/Q’) Yn? | A. 7 —Q z C,( x Awa) |/0 +1 A. f. (11) 


From the definition (10) of C, and A, , it is clear that, in the numerator of (11), the 
n occurring in the first term may be replaced by Pn, while the n and k of the second 
term may be replaced by Qn and Pk respectively, so that the summaticn occurring in 
the second term of the numerator of (11) becomes : Ca $3 A ppAgn—p,x). Since P 
and Q are relatively prime, Zs A ppAgn-px = O unless nis a multiple of P. But Agp,-p, = 
A pren-x) ; hence, from (10), (11) and (7), there follows 


( 
u = g.Lbd| rP*/Q’) > n? | a, |? — Ee | 2 avte,-) | te rt (12) 
k 


n 


In the expression (12), the arbitrary positive integers P and Q and the finite set of com- 


plex numbers a_y , @-y.;, *** , @y are subject to 


(P, Q) = 1, Ain = Un and ps a, = 0. (13) 


—n 
n=—N 


(It is to be noticed that each summation of (12) extends over only a finite range.) The 
formula (12), subject to (13), will be used in the next section to deduce (4). 


3. Choose a, = 1/27, a, = —1/27, a, = Oif n ¥ +1; let Q = 1 and let P denote 
an arbitrary positive integer. (This selection corresponds to the choice x = sin 2rP¢t in 
(8).) The expression {---} of (12) becomes 2°(2P)* — co + R(c2r) so that the equation 
(4) is proved for an arbitrary even positive integer N = 22. Suppose now that N of 

1) is odd and choose P = N and Q = 2; let the a,’s be defined as before. Proceeding as 
above one sees that relation (12) implies » S 2°N* — co + R(cy) and the proof of (4) 
is complete. 

Appendix 


I. For computational purposes, one may easily verify that the formula (12) can 


be modified to 


" 
w= —t) + g.lb| @r'P* Q”) ps n’ | a, |” 


n=1 


~an( S( Zanes))]/a+2 E101) 


For appropriate choices of the sequences {a;} and the pairs of integers P and Q, subject 
to (13), various upper bounds for yu, in addition to those of (4), may be readily obtained 
from the above formula and (13). Non-trivial lower bounds for yu are not as easily ob- 


g.l.b.” in the 


a 


tained, in this manner, corresponding to the presence of the symbol 


formulas for yp. 
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II. The question has been raised by Wintner [10] whether the constant 2 occurring 
as the coefficient in (3) is the least value of a for which 


p> —m—a Dd) Ic, |’ (14) 
n=1 

holds for an arbitrary periodic function f(é) defined by (2). Although this question will 

remain unanswered, it can easily be shown, as a consequence of (4), that a = 1/47”. 

For, suppose (14) holds for all f(t) defined by (2); then, by (4), —a D°2., |c, | 


lA 


ax N° + R(cy) holds for N = 1, 2, --- . If cy is real, it follows that r°N* + cy + acz = ( 


hence, by a consideration of the discriminant of this last quadratic expression, |! 
4anr°N* < 0. For N = 1, this implies a = 1/42°, which was to be shown. 


— 
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A MODIFICATION OF SOUTHWELL’S METHOD* 
By W. H. INGRAM (New York) 

J. L. Synge’ has given a geometrical interpretation of Southwell’s method of solution 
of the problem Az = b when A = (a,;) is symmetric and }> >> a,;2,2; is a positive 
definite form. A modification of the method having application to the more general 
case in which 2A ,Az is a positive definite form makes use of the ellipsoids of the Gauss- 
Seidel process. 

For any vector 2, there is an error e defined by the equation 


Az —- b=e, (1) 


therefore 


(cAy — b)W(Azx — b) = eWe; (2) 


*Received Jan. 2, 1951. 
1J, L. Synge, A geometrical interpretation of the relaxation method, Q. Appl. Math., 2, p. 87 (1944). 
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W is a diagonal weighting matrix for weighting the relative importance of the equations 
of the set (1). An orthogonal transformation 7' exists such that 
TrA;,WAT = D=[d,,d2,--- , d,) 
is a diagonal matrix of positive elements d; . The substitution 
t=Ty «=yTr, 

into (2) gives 

dy; + doy> + dsys + +++ — 28,4, — 2Bsy2 — 2Bsys — --- = eWe — bW, (3) 


an equation in which 8 = T,A,Wb. It is seen that (3) is the equation of a family of 
hyperellipsoids with the properties 
(a) a common center at (6,/d; , B2/d2, -++) 
(b) a common orientation, 
(c) common principal axes’ ratios, 
(d) parameter e. 
To complete the squares, one adds Zz. B;/d; to both sides of (3): 
> 63/d; = BD"'8 = bATD"'T;zAzb = bWb. 
It is seen that the ellipsoids converge to a point as e — 0 and that their common center, 
given by e = 0, is the solution of (1). 
The solution may be approximated in the same way as for the family of ellipsoids 
employed by Southwell and by the same process. 


THE LAPLACIAN AND MEAN VALUES* 
By R. M. REDHEFFER anp R. STEINBERG (University of California, Los Angeles) 


Introduction. The Laplacian V’*f represents deviation from the average in the following 
sense [1]. Let f, be the mean value of f over a cube of edge 2, 


8? f. = ill f(x + u, y +0, 2+ w) du dv dw. (1) 
Then 
lim 6(fo = f/t - V's, (2) 


provided f is sufficiently smooth. This result is interesting in that it is independent of 
the co-ordinate system. Also it sheds a certain light on the wave equation. Thus, if the 
restoring force at a point in a medium is proportional to the deviation from the average, 
in some sense, then one might expect the equation of motion to be a V* f = bf., , where 
a measures stiffness, b inertia. 








*Received April 9, 1951. 








316 NOTES (Vol. IX, No. 3 

\ harmonic function could be characterized as one such that the limit (2) is zero, 
in other words, such that its value at a point is close to its average over a suitable neigh- 
borhood. This seems interesting too, since a harmonic function can be characterized by 
an exact equality of that sort over the surface of a sphere. 

A result of the type (2) is true in one or two dimensions, and also when the average 
is taken over the surface, edges or vertices of the cube, rather than over its volume. 
From (2) for the boundary of a square in two dimensions one can deduce the fact that 
the line integral of a harmonic function in two dimensions is independent of path. Again, 
(2) over the vertices of a square or cube is related to methods of numerical integration 
[2], where one wishes to find f,by computing f or its derivatives at a point. 

The Problem. What we want to do here is to characterize the point sets for which a 
result like (2) holds. It is seen from [1] that symmetry in the coordinate planes through 
(x, y, 2) is sufficient, but it is not necessary. For a solid set, / is taken as a triple integral, 
divided by the volume of the set; for an area one uses a double integral divided by the 
area; for a space curve one uses a line integral divided by the length; and for a finite 
point set one takes the sum of the values, divided by the number of the points. Concise- 
ness requires that the result be given here in a form that applies to all these cases at 
once, and to any number of dimensions. Accordingly we denote the number, length, 
surface, or volume by a general measure m: 

I. Let. S be a bounded point set on which is defined a measure m such that 0 < m(S) = 
M <o. If the vector P represents a fixed point of S, and the vector Q a variable point 
of S, the following statements are equivalent: 

(a) The set S has mass J/, its center of muss is at P?, and its moment of inertia about 
every axis through P is J. 

b) f 1{/P+ 4«(Q P)| din(Q) — Mf(P) ~ (t/2)°T VW f(P) ast > 0 for an arbitrary 

function f having continuous third derivatives near P. 

We do not dwell on the proof, which is rather trivial. There is no restriction in taking 
P = 0. Then, Taylor’s theorem with remainder about P = 0 (in Cartesian co-ordinates 
with Q = (zx, y, z)) shows (b) to be equivalent to the mass and center of mass results of 
(a) and 


| yz = | zr = | xz = 0. (4) 


In each integration Q ranges over S.) 


The identity 2? = [(2° + 2°) + (a7 + y°) — (y° + 2)]/2 shows that (3) are equivalent to 
| (y? + 2’) = | (2? + 2”) = | (2? + y’) = I. (5) 
Equations (4) and (5) are equivalent to the identity in /, m, n: 


| [((mz — ny)? + (nx — lz)? + (ly — mz)’) = (PF + mM’ 4+ n’)I; 


e 


and this is the moment of inertia condition of (a). 
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\s for the speed with which the limit in (b) is approached, we get, again using Taylor’s 


expansion, the following: 


II. In I(b), the error is O(¢°). If f has continuous fourth derivatives near P, the error 
can be made O(/') by a further restriction on S. Whatever set S is chosen, the error 
cannot be made o(¢') even if the class of functions is restricted to polynomials. 


REFERENCES 
1. L. Hopf, Differential equations of physics, Dover Publications, 1948, p. 62. 
2. Garrett Birkhoff and David Young, Numerical quadrature of analytic and harmonic functions, J. of 


Math. and Physies 29, 217-221 (1950). 


A NOTE ON ASYMPTOTIC STABILITY* 
By H. A. ANTOSIEWICZ (Montana State College) 


1. In this note we shall develop a stability criterion for a vector differential equation 


of the form 


dx 
= A(é)x, (1) 


where the elements of the matrix A(t) = (a,;(), 7,7 = 1, 2, --+ , are real continuous 
and uniformly bounded functions for all positive 4 = t . 

A. Wintner** recently established the following criterion: Let \,(¢) be the greatest, 
and .(/) the least characteristic value of the matrix 3[A(4) + A’()], and let || x() 
denote the Euclidean length of the vector x(é). If [° \,(0) dt <@, f° d.(0) dt <@, then 

r(t) >« ~ 0ast—o- for every non-trivial solution x(¢) of (1). 

It is to be noted that the condition of integrability of 4, (4, A2(¢) over (¢) , ©) implies 

f° [trace A(t)| dt < ©. Furthermore, this condition automatically excludes the important 


ease A(t) = const. unless A(t) = const. is skew-symmetric. 
In the following we shall establish a stability criterion which is free of the above 
objection, i.e. which will also apply to the general case A(t) = const. We shall consider 


a condition to be satisfied by the matrix A(é) which will suffice to insure that || x(d) || 
of every non-trivial solution 2(¢) of (1) tends to zero as t >. According to Liapounofff, 


the trivial solution x(/) = 0 is then said to be asymptotically stable. 
2. Consider a funetion V(x, t) which is defined and continuous for all x and ¢ in R: 
x,|S<e¢t=T7= 1,2, --- n). If for equation (1) there exists in R a function V (2, ¢) 


which is of fixed sign and admits of an infinitely small upper bound, and for which 
dV /dt by virtue of (1) is opposite in sign to V (a, 2) in R, then the trivial solution a(t) = 0 


of (1) is asymptotically stable. Liapounoff proved that the existence of such a function 


*Received April 16, 1951. 
**A Wintner, On free vibrations with amplitudinal limits, Quart. Applied Math. 8, 102-10. (1950). 
tA. Liapounoff, Probléme général de la stabilité du mouvement, Ann. Math. Studies, No. 17, 1949. 








318 NOTES [Vol. IX, No. 3 


V(z, t) is sufficient for asymptotic stability; it is, however, not necessary as was shown 


by J. Malkin.* 
We shall make use of Malkin’s results to establish the following theorem: 


Let \,(t) be the greatest, and d.(t) the least characteristic value of the matrix 3{[A(t) + 


A'(t)]. If f° (7) dr — —©, J 2(7) dr — —@ ast >, then || x(t) || ~ 0 as t@ 
for every non-trivial solution x(t) of (1), 7.e. the trivial solution x(t) = 0 ts asymptotically 
stable. 

Note that now {* [trace A (é)]dt diverges. 

3. First, we transform (1) into diagonal form. Let 2, , x2, --- 2, be a base of solutions 
of (1), and use this base to construct an orthogonal matrix C(t). If y = C™'(é)z, then 
(1) reduces to 

Y~1 
oY = By, Bit)=C7AC+ C (2) 
where the matrix B(i) = (6;;(), 7,7 = 1, 2, --- n, is diagonal, i.e. 6;;(4) = 0 for all 
it > j. If y:(2), yo(t), «++ yp(t) is that base of solutions of (2) for which y,(é) = J*, the 
i-th column vector of the identity matrix J, then || z;(¢) || = || y.(0) || as is easily verified. 
Evidently, C(t) and C~*(t) have bounded elements and | C(é) | = | C~*(t) | = 1; hence 


stability properties are preserved in both directions. 

Observing that C~’(#) = C’(é) by construction, we find by differentiating the identity 
C(t)C7*(t) = I that (dC~*/dt)C is skew-symmetric. Therefore B(t) + B’(t) = C™’[A() + 
A’(t)]C, and thus the characteristic values of 3[B(t) + B’(t)] are identical with those 


of 3[A(#) + A’(é)]. Hence it is sufficient to prove our theorem for the reduced equation 


> 


(2). We shall show that there exists a function V(y, ¢) which satisfies Liapounoff’s 
criterion for asymptotic stability. 

Consider the diagonal elements b;;(é), 7 = 1, 2, --- n, of the matrix B(t). Since 
(dC~*/dt)C is skew-symmetric, trace (dC~*/dt)C = 0, and thus 
b(t) = (C7'),AC* = (C')'AC* = 4(C")'[A(O) + A’(OIC’. (3) 


All diagonal elements of B(t) are quadratic forms in the components of the column vectors 
= 1. These quadratic forms 


C* of the matrix C(t) for which we evidently have || C’ | 
|| = 1 (compact set); if 


attain their maximum and minimum on the unit sphere || C’ 
i, (t) is the greatest, \,(¢) the least characteristic value of 3[B(t) + B’()], then A,(® is 


the maximum, \.(/) the minimum. From (3) we then obtain 


i () = b(t) = ro(0) (4) 
whence for all ¢ = ¢ 
exp (fi, Ai(r) dr) = exp (fi, b::(r) dr) = exp (Jt, A2(7) dr). (5) 
By hypothesis {‘ \,(7) dr ~ —” ast—~, k = 1, 2, and thus 
y(t) = exp (Ji, b::(r) dr) + 0 as t—>o, (6) 


As Malkin has shown, (6) involves for all ¢ = t, 


nt d " 
0 | ody ¥° @ 


*J, Malkin, Certain questions on the theory of the stability of motion in the sense of Liapounoff, American 





Math. Soc., Translation No. 20, 1950. 
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and (6) and (7) together, in turn, imply f° [g;(0)]? dt < ©. Hence the functions 


. 


VO) = lelOl* | lede)F de (8) 


exist for all £ = %& and are uniformly bounded; in fact, a* < y,(é) < 6° where a and b 
are certain constants. 
Now consider the function 


V(y, 2) = Wldyi + v(dys +--+ + vA(Dye - 


It evidently satisfies Liapounoff’s criterion for asymptotic stability; it is a positive 
definite quadratic form, admitting of an infinitely small upper bound, and its derivative, 


by virtue of (2), becomes 


dV 


ap = Ute t es ty) + WY, 


where W(y, 2) is a quadratic form whose coefficients depend upon those elements b, ; (¢) 
of B(t) for which 7 < 7, 7,7 = 1, 2, --- n. Since these elements can always be made 
sufficiently small by a transformation with constant coefficients (which will not affect 
stability properties) the derivative dV/dt will be a negative definite quadratic form. 
Hence the trivial solution y(t) = O of (2) is asymptotically stable, and therefore the 
trivial solution x(t) 0 of (1) is asymptotically stable. This establishes our theorem. 
O. Perron* was the first to prove directly that the conditions 
dr 


r(t) $C | — 
a ) = 1») Y y I g(r) 


IIA 


C, 


| 


are necessary and sufficient for the trivial solution x(t) = 0 of (1) to be asymptotically 
stable. 


*O. Perron, Die Stabilitaetsfrage bei Di fferentialgleichungen, Math. Zeitschrift 32, 703-728 (1930). 


CONDITIONS SATISFIED BY THE EXPANSION AND VORTICITY OF 
A VISCOUS FLUID IN A FIXED CONTAINER* 


By J. L. SYNGE (Dublin Institute for Advanced Studies) 


1. Introduction. In plane motion of a viscous fluid inside a fixed container, the 
expansion 6 and the vorticity w cannot be arbitrarily assigned. A necessary and sufficient 
condition’ for the consistency of given @ and w with vanishing velocity on the walls is 


[ ov + 2wV) dS = 0, (1.1) 





*Received April 17, 1951. This paper was written while the author was on leave of absence at the 
Institute for Fluid Dynamics and Applied Mathematics, University of Maryland. 

1J, L. Synge, Quarterly of Applied Mathematics, 8, 107-108 (1950). The condition with 6 = 0 was 
originally due to G. Hamel, Géttinger Nachr. Math.-Phys. KI. 1911, 261-270. 
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where U’, V is any pair of conjugate harmonic functions and the integration is taken 
over the region occupied by the fluid. 

The purpose of the present paper is to extend this result to three dimensions. In 
space we have expansion @ and a vorticity vector w; (suffixes range 1, 2, 3 with the 
summation convention), and the theorem which will be proved may be stated as follows: 
Given expansion 6 and vorticity w,; are consistent with vanishing velocity on the walls if, and 
only 7f, 


| (OP + w.Q,) dV = 0, (1.2) 
where the integration is taken through the fluid, Q; being any (arbitrary) solution of thé 


partial diffe ential equations 
AQ; = Qh.xi (1.3) 


and P satisfying 
Po, = $6€5j:Q:.; - (1.4) 


Here ¢,;, is the usual permutation symbol, the comma denotes partial differentiation 
(Y., = 0Y/dér;,), and A the Laplace operator, so that AQ; = Q;,,, . In vector notation 
(1.3) and (1.4) read V°Q = VV -Q and VP = 3V X Q. 

This theorem will be considered only for simply connected regions. For such regions, 
(1.3) are precisely the integrability conditions of (1.4), so that, given any solution of 
(1.3), P exists satisfying (1.4), unique to within an additive constant. 

2. Necessity of condition (1.2). Expansion 6 and vorticity w; are connected with 
velocity u; by 


u = Q, € u = 2w j (2.1) 


The problem of finding a motion with given 6 and w, in a region V, bounded by a fixed 
surface B to which the fluid adheres, is the problem of solving the partial differential 


equations (2.1) for u, with the boundary condition 


u; = 0 on B. (2.2) 


The theorem stated above asserts that (1.2) is a necessary and sufficient condition on 
6 and w, for the existence of this solution. 

The necessity of the condition is easy to prove. We assume the existence of a solution 
of (2.1) and (2.2). 

By virtue of (2.1), it follows that for any P and Q; at all (not subject to any condi- 
tions save those of smoothness) we have 


> 


[ (oP + 0.0.) dV + | uP. — desinQe.s) AV 


e 


= (uP). + wiQ; — (heijuiQs).; + F essus,;Qi] dV (2.3) 


© 


— | (nu;P — 3 nje;;u;Q,) dB, 


* 
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the last integral being taken over the boundary B, on which n, is the unit normal, drawn 
outward. 

Note that the above follows from (2.1) only. We now bring in (2.2). This makes the 
last integral in (2.3) vanish. If we then subject P and Q, to (1.3) and (1.4), the second 
integral in the first line vanishes, and we are left with the equation (1.2); the necessity 
of (1.2) is thus established. : 

3. A lemma. The proof of sufficiency is harder. It rests on a lemma, for which the 
proof offered here is not mathematically rigorous, resting as it does on the assumption 
that a certain minimum is attained. A precise mathematical proof would of course have 
to specify the requisite smoothness of the bounding surface B and of the tangential 
component assigned on it (see immediately below). 


Lemma: Given the tangential component of a vector Q; on the boundary B of a region V, 


then Q, exists satisfying (1.3) and this boundary condition. 


To prove this (or at least make it plausible), consider the integral 


1Q) = | iinQu.ieirnQe AV. (3.1) 
In vector notation, the integrand is (V X Q)° and cannot be negative. Thus for all 
Q, satisfying the stated boundary condition, 7(Q) is bounded below. We assume that 
the minimum is actually attained by some vector field; let Q; be it. 
Then, if ¢ is any constant, and g; any vector field with zero tangential component 
on B, it follows that 


I(Q) < 1(Q + eq), (3.2) 


and hence by the usual procedure associated with Dirichlet’s principle, 


” 


| €:itQe.stirds.r dV = 0 (3.3) 
for all such g; . This may be transformed into 


+ 


f e:iQs.i¢iregan, dB — | €:i1Qe.ic€irege AV = 0. (3.4) 
But €;,.9.m, is the tangential component of g; , turned through a right angle, and so 


vanishes. Further 
€csntire = 6;, Oks ts 0; Ox, ; (3.5) 


and so we get 


e 


| (aQ. — Qi.:)% dV = 0. (3.6) 


Since g; is arbitrary except for the boundary condition it follows that the minimising 
(); satisfies (1.3). 

As pointed out already, P then exists satisfying (1.4). 

4. The sufficiency of (1.2). We assume that @ and w; satisfy (1.2) for all Q; and P 
satisfying (1.3) and (1.4). We have to prove the existence of a solution u; of (2.1) and 


(2.2) 
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Choose a particular solution of (1.3), (1.4): P = 1, Q,; = 0. Then (1.2) gives 
edV =0. (4.1) 


Now choose P 0, Q, W,; , the gradient of any scalar field W; these satisfy 
(1.3), (1.4) without restriction on W. Then (1.2) gives 


| wW..aV = 0, (4.2) 
or 
| win,W dB — | w.WaV =0. (4.3) 


Hence, in view of the arbitrariness of W, 
W; ; 0 in V, w,n; = 0 on B., (4.4) 
We now try to solve (2.1) and (2.2). This we do in two steps, first solving the system 
= 6, €;jata,; = ow; , (4.5) 


with the boundary condition 


un, = 0 on B. (4.6) 


There is a well known procedure for this; we express u; in terms of a scalar potential 
and a vector potential, and obtain for these certain Poisson integrals, yielding particular 
solutions of (4.5). The problem of solving (4.5) and (4.6) is thus reduced to a Neumann 
problem, and it is known that (4.1) and (4.4) are sufficient conditions for the existence 
of a solution to (4.5) and (4.6). 

The second step is to prove that the tangential component of this u, vanishes, the 
normal component being already zero by (4.6). This is easy. Equation (2.3) is valid, 
since it depends only on (2.1), i.e. (4.5). But the first two integrals in (2.3) vanish by 
(1.2) and (1.4) respectively. Also the first part of the last integrand vanishes by (4.6). 
So we are left with 


7 


| n,€;;,u,Q, dB = 0. (4.7) 


< 


But ¢,;,7;Q, is the tangential component of Q,; , turned through a right angle, and 
this, as we saw in the lemma, may be chosen arbitrarily. From this it follows at once 
that the tangential component of u, must vanish on B. 

Thus the solution of (4.5) and (4.6) is in fact the solution of (2.1) and (2.2). This 
completes the proof of the sufficiency of the condition (1.2), subject of course to the 


assumption that the minimum of (3.1) is in fact attained.* 

*Added in proof, June 22, 1951: In recent papers which the author had not seen when the present 
paper was written, C. Truesdell, Comptes Rendus Ac. Sci. Paris, 232, 1277, 1396 (1951), has given the 
condition (1.2) with consideration of its sufficiency, and this condition has also been given by F. H. 


van den Dungen, Q. Appl. Math. 9, 203 (1951), but the question of sufficiency has not been fully con- 


sidered by him. 
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CONFORMAL MAPPINGS FOR CERTAIN DOUBLY CONNECTED DOMAINS* 
By J. A. McFADDEN (University of Michigan)** 


1. Introduction. In this paper conformal mappings are derived for some doubly 
connected domains, the inner boundary being formed by slots and the outer boundary 
being the unit circle. All these mapping functions utilize elliptic functions; those for 
the more complicated domains involve elliptic functions of two different moduli. The 
domains appear, for example, in the linearized theory of supersonic flow past conical 
bodies. 

2. Mapping of domains bounded by a slot and the unit circle. Consider first an 
annular domain in the ¢,-plane, r; < |¢,| < 1, and the corresponding rectangular 
domain in the z,-plane, —2K, < 2, < 2K,,0 < y, < K{, obtained by a logarithmic 
transformation from the ¢,-plane. (See Fig. 1. Circular domains will be designated by 
Greek letters, ¢; = & + 7m; , and rectangular domains by Roman letters, z; = x; + 7y; .) 


= (2K,/m7) log ¢, (0 < arg ft; < 2m) (1) 


The radius 7, of the inner circle in the ¢,-plane is given by the expression 


r, = exp (—7K{/2K,). (2) 
The cirele | ¢, r, can be collapsed onto a horizontal slot, —(1 — k{)/ki < # < 
(1 k{)/k, , n° = 0, (see Fig. 1), with the circle | ¢; | = 1 mapping onto the circle 
| ¢ |, by means of the transformation 
., _. en (2, 3 ky) + tky sn (2, ; ky) (3) 
: dn (2, ; k,) . ‘ 
where sn (2, ; k,), en : k,), and dn (z, ; k,) are the Jacobian elliptic functions with 


argument z, and modulus k, , and where kf = (1 — k{)'”’. (See, for example, reference 
(1].) The constants K, and K{ in equation (2) may then be identified as the complete 
elliptic integrals of the first kind having moduli k, and kf , respectively. An alternative 
form of equation (3) is the following: 
—(1 + ¢”’)/2¢’ = sn (z, — K, ; ki) 

(3a) 


(—2K, <2, < 2K, : 0<y, < Kj) 


In this form we can understand the reason for the appearance of the Jacobian elliptic 
functions. The left member of (3a) maps the two boundaries in the ¢’-plane onto the 
real axis of an intermediate plane. The right member of (3a) maps the rectangular 
boundary in the z,-plane onto the same real axis, since the sn-function is real on all 
four sides of the rectangle. 

Note: This transformation is related to the transformation ¢’ = —(k))'”? sn (zo ; ko) 
by a Gauss transformation of elliptic functions (see, for example, reference [2]) and a 
translation of the coordinate axes. In the aerodynamic application mentioned above, 
the form (3) is found to be much simpler. The related form is discussed by Holzmiiller [3]. 

*Received April 6, 1951. This paper is a chapter of a dissertation submitted in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy at the University of Michigan, December, 1950. 
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The transformation (3) and equivalent combinations of functions have been used 
by previous investigators in the treatment of the linearized supersonic flow past a delta 
wing. (See, for example, reference [4].) We shall use it here as a basie transformation, 


from which more complicated transformations can be derived. 
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Fic. 1. Conformal Mapping of Domains Bounded by a Slot and the Unit Circle. 
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In equation (3) we may replace ¢’ by 7¢ and z, by (z, + K,). The resulting transforma- 
r, onto a vertical slot € = 0, —(1 — kf)/k, < 9 < 


tion collapses the circle | ¢ é 


I k')/k, (see Fig. 1). The transformation function is 
¢ = en (z, ;k,) + 728n (2, ; k,). (4) 


An alternative form of equation (4) is the following: 


(1 — ¢°)/2¢ = sn (2, ;s k,) 
(4a) 
(=—2h, <2) < 2a, 0 <9; < KD 


Che transformation (1), followed by either (3) or (4), leaves the included segments 
of both axes invariant as well as the unit circle. (See Fig. 1.) 

3. Mapping of domains bounded by a cross and the unit circle. Because of the invari- 
ance of the axes, the horizontal slot may be opened into a circle with symmetric hori- 
zontal fins in the ¢,-plane (see Fig. 2). We apply the inverse of transformations (1) 
and (2), but with a smaller modulus k, < k, . That is, let 
en (z; ; ks) + tks sn (23 ; ks) . 
ee (5) 

an (s, * Ks) 


~S 


¢ = 


’ 


where 
z, = (2K;/mi) log ¢; , (0 < arg f; < 27m). (6) 
The corresponding domain in the z,-plane is a slotted rectangle (see Fig. 2). 
The radius of the inner circle in the ¢,-plane is r, , given by the equation 
r, = exp (—7K3/2K;). (7) 
An alternative form of equation (5) can be written, similar to equation (3a). 


—(1 + ¢”)/2¢’ = sn (23 — Kz ; ks) 
(5a) 


(—-2k. <a < Bh: 0 < yy; < K4) 


If we combine equations (3a) and (5a), we may map directly from the z,-plane to 
the z,-plane. The resulting transformation may be written in four forms: (Related 


functions have been given by Kronsbein [5].) 


sn (2, — K, ;k,) = sn (@, — Ks ; ks) 
en (2, — K, ;k,) = m (@; — Kz ; ks) 
dn (2, — K, ;k,) = +[1 — kisn’ (2, — Ks ;3ks)]'” (8) 
+[1 — kj sn’ (z, — K, ;k,)]'” = dn (z; — Ks ; ks) 


(—2K, < 2, < 2K, ;0 < y, < Ki; — 2K; < 2; < 2K; ;0 < ys < K3) 


= ? 


If we argue in similar manner with equation (4a), we may construct a domain 
bounded internally by a circle with symmetric vertical fins in the {.-plane. (See Fig. 2.) 
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Fic. 2, Conformal Mapping of Domains Bounded by a Cross and the Unit Circle. 
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Again a new modulus k, < k, is used. The resulting transformation may be written in 
four forms: 
sn (z, ;k,) = sn (22 ; ke) 


en (z, ;k,) = en (22 ; ke) 
dn (z, ;k:) = +[1 — kjsn’ (2 ; k.)]'” (9) 
+[1 — kjsn? (z, ;k,)]' = dn (2 ; ke) 
(—2K, < 2, < 2K, ;0 < y, < Ki; — 2K, < 2 < 2K, ;0 < y, < Kj) 


where 
2. = (2K,/xi) log & , (O < arg {2 < 2rz). (10) 
The corresponding domain in the z.-plane is again a slotted rectangle. (See Fig. 2.) 
The radius of the inner circle in the {.-plane is r, , given by the equation 
T, = exp (—7K}/2K,). (11) 
A transformation similar to equation (3) will collapse the circle | ¢,| = r, onto the 
real axis in the ¢,-plane, leaving the vertical fins on the imaginary axis. (See Fig. 2.) 
The transformation equation is the following: 





cn (2 ; ke tki sn (2. ; ke) 
fp = EEE a ite (12) 
dn (z2 ; ke) 
A transformation similar to equation (4) will collapse the circle | ¢; | = 7; onto the 


imaginary axis in the ¢,-plane, leaving the horizontal fins on the real axis. (See Fig. 2.) 
The transformation equation is the following: 
t. = en (z; ; ks) + 2 sn (2s ; ks) (13) 
The two definitions of ¢, in equations (12) and (13) are equivalent if and only if 
the moduli are related by the equation 
ki = kzk3 . (14) 
[This statement can easily be proved in the form (17).] Then the inner boundary in the 
¢,-plane is a cross consisting of the horizontal slot, —(1 — k3)/k. < & < (1 — ki)/ke, 
n, = 0, and the vertical slot, &£, = 0, —(1 — k3)/ks < m < (1 — k3)/ks . 
By a series of steps we have mapped the circle | ¢; | = 7, (Fig. 1) onto a cross, the 
unit circle remaining invariant (as well as the coordinate axes). 
It is convenient to introduce a shape factor e for the cross, defined by the relation 
e = ki/(ki)'? = (ki)? /ks . (15) 
e has the range (k{)'’* < « < (k{)”'”. Then the moduli k, and k; and the respective 


complementary moduli are given in terms of ¢ and k, by the equations 


(1 — ki)”, ks _ e(ki)'”, 


k, 
(16) 


- 2 
/€. 


kz = (1 — ki/é)™, ks = (ki) 
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If « = 1, then ki = ki = (k{)'” and the vertical fins have the same length as the 
horizontal fins. If « = (k{)~'’*, then k, = 0, k; = k, , r2 = 0, and equations (8) become 
an identity. Then the cross degenerates to a vertical slot. If « = (k{)'””, then k. = k, , 
k; = 0, r; = 0, and equations (9) become an identity. Then the cross degenerates to 


a horizontal slot. In general, if « > 1, the vertical fins are longer than the horizontal 
fins. If e < 1, the horizontal fins are longer. (Fig. 2 is drawn for e < 1.) 
The combination of equations (9) and (12) or of equations (8) and (13) gives the 


resulting equation 


en (z, :k,) + ze(ki)’”” sn (2, : ky) 


l 


wre ; (17 
en (2, ;k,) — ie(k:)'”” sn (2, ; ky) } 


Although the function ¢; is an elliptic function, the function ¢, is not. It is doubly periodic, 
but it has branch points. In particular, there are branch points at the images of the 


point ¢, = 0, the intersection of the cross. These points are defined by the equations 
fees sea ok, ; 
z= t(5K, +7 + iK{, + 9 + 7iK{. (18) 
[These points are the zeros of the functions dn (z, ; k,) and dn (z; — K; ; k3). See Fig. 2.] 
Ife< 1,7 >0O0.Ife > 1,7 < 0. If « = (k{)”’, then s = 3K, and the oe points 
coalesce to form simple zeros at the points z, = +A, + 7A{, the zeros of the function 
3). If e = (kf) “~, then r = —3A, and the branch points coalesce to form simple zeros 

at the points z, = 7K{ , 2K, + 7K{, as in the function (4). 


Let us derive the mapping from the ¢,-plane onto the ¢,-plane by a distinctly different 
method. First, we square the variable ¢, , making the two circles doubly covered, the 
inner one being of radius 7; . Let the appropriate modulus be k, , so that the corresponding 


complete elliptic integrals are related to the radius by the equation 
r; = exp (—7K{/2K,). (19) 


Second, apply the transformation (3) with modulus k, in the two-sheeted plane, pre- 
serving the doubly covered unit circle and creating a quadruply covered horizontal slot. 
Third, extract the square root. The final boundaries are the simply covered unit circle 
and a cross with four fins of equal length. The resulting transformation may be written 


as follows: (A related function has been given by Kronsbein [5].) 


en ((4K,/ri] log ¢ ; ki) + aki sn ((4K,/r7] log £1 ; ki) (20) 
dn ([4K,/mi] log ¢, ; k,) " 


fel 


(C4) = 


Comparison of equations (2) and (19) reveals that the moduli k, and k, are related 
by a Landen transformation. (See, for example, reference [6].) If we apply the Landen 
transformation, whereby 
k= 11 -k)/14+ £), ki = 2(k)'?/(1 + ki), 
(21) 


K, = 401 + kDK, , Ki =(1+k)K, 


and 
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sn (2, ; k,) en (2, ; ki) 
dn (2, ; k,) 





sn ({1 + kijz, 34) = (1 + &) 


’ 





a i ng eh) — LTO +k) snt & 5 hy) 2 
en ({1 + kijJa, 3 *) = iG 143 ; (22) 

| 1 — (1 — ki) sn’? @ 5h) 

dn ({1 + kz, :k,) = — 1 
In ({1 + £7] k) ist kD : 

then equation (20) may be written 
| » ek */T9\1/2 - d- 
(i) = en (2: 5 ki) + (hi) sn (2 5 hs) (23) 


en (z, ;k,) — i(k)’ sn (2, ; &,)’ 


where 2, is given as in (1). 
We observe that equation (23) is indeed the mapping function for a cross of equal 
lengths, for it is a special case of equation (17), namely e = 1. 
The cross of equal lengths can easily be mapped into the cross of unequal lengths in 
Fig. 2. We apply a linear fractional transformation which preserves the unit circle in 
the two-sheeted plane. If the transformation 


(ci)? — (e — D/e + D 
—n 7 3 24 
i —[e- De + DIK (24) 


is combined with equation (23), then equation (17) is the result. 





2 
4 


¢ 
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EFFECT OF A RIGID ELLIPTIC DISK ON THE 
STRESS DISTRIBUTION IN AN ORTHOTROPIC PLATE* 


By A. J. OWENS anp C. B. SMITH, (University of Florida) 


A thin orthotropic plate of uniform thickness will possess two perpendicular axes 
of symmetry in the plane of the plate. An infinite rectangular plate of this type con- 
taining a rigid elliptic disk with major and minor axes coinciding with the axes of sym- 
metry is discussed. A uniform tension is assumed to act along two opposite edges of the 
plate and a mathematical analysis of the stress distribution is given. It is assumed the 
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strains are small and remain within the limits of perfect elasticity. The solution obtained 
is applied to a plain-sawn Sitka spruce plate. 

Choose as the origin of the coordinate system the center of the ellipse. The boundary 
conditions may be stated mathematically using the notation of A. E. H. Love’ as 


SS a Y, | .. = 0,) 


( (1) 
ee eee x = ©: 


and u(x, y) = v(x, y) = 0 on the circumference of the disk given by the parametric 
equations x = a cos 6, y = b sin 6. The displacements in the x and y directions are 


u(x, y) and v(x, y), respectively. S is the uniform tension applied at two edges of the 
plate. 

The components of stress and strain are connected by the following relations if the 
x and y axes are taken as the axes of elastic symmetry of the orthotropic plate: 


Ou l oo 
OS ee 9 
ee £. ae E, Vy (2) 
Ov —— .§ : 
— dy EE, td E, Vy, (3) 
ey = 4 Ba ,, (4) 


Oy Ox Hes 
It is desirable to find a stress function F(x, y)° such that 
oF : OF . OF 


‘ A, * ee, 
ne . Ox Oy 





(5) 


XxX, = " ,= ;, 
; oy” . Ox” 


For the problem of a thin orthotropic plate in a state of plane stress F must satisfy the 
differential equation‘ 


o*F . OF oF 

— 2K sss +@ si = 6 

a* +o Ox On On ' (6) 

where 

: 2 . 

K = 4(E,E,)'” (+ - a) (7) 

Mery E, 
n = &y, and e = (E,/E,)". (8) 


A suitable stress function is 


r= RA 
\2 


71 


iF (Z, — Wi)? + vi In (Z, + rv) | 


3 Be ) 
- c: 2 E (Z. — W2)* + v2 In(Z, + w.) | (9) 
2Y¥2 LS 


Sn° 
Qe" 


a 


1A. E. H. Love, A treatise on the mathematical theory of elasticity, Dover Publications, New York, 1944. 
2H, W. March, Stress-strain relations in wood and plywood considered as orthotropic materials, Forest 
Products Laboratory Rep. No. R1503, p. 2 (1944). 
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where the symbol R{ | denotes the real part, 

A=A,+7A,, B= B, + 7B,, (10) 

Z, = x + ian, Z, = x + 7Bn, n = €&, (11) 

W,=(Z4i-%7)'", W.= (4-1), (12) 

71 = a’ — a’e’d’, v2 = a” — B’e’b’, (13) 


a = {K + (K? — 1)'“}", 


and 


2 


B = {K — (K’ — 1)'"7}"”. 
Using the transformation 


K = cosh 4, 
(14) 


a =e” and B=e°, 


In order that the stresses may be single valued W, and W, are to be assigned values 


so that the inequalities 





Z.+W,l2n and |Z, + W.| 2 v2 (15) 
are satisfied. 
It follows that 
, Oe. og 4 a. ee 
Y= 3" Mine +W,) + Wid + a} (16) 
: ,OF ae A B’eB } . - 
= “4 : = —_—_——O —_—_—_—_—_——— S 
An = ¢ oH wi +W)*Wi.+Wos a7) 
and 
__ _ OF _ pf iaeA wie} 
X, = —tor an ~ “\Wid, + Wy) + Wie, + Was” (18) 


It can be shown that the exterior boundary conditions are satisfied so long as A and 


B are finite. 
Substitution of the stresses given by (16) and (17) into (2) and integration gives 


e ate’. Oye |] 1 Bre’ | Gey |} __—1 } SE, 
oh al +E dx + "i ba al z* hs + Wilf tz, + 





where u(y) is an arbitrary function which can be shown to be identically zero. A similar 
expression for v may be found by substituting the stresses in (3) and integrating. From 
the conditions vu = v = 0 on the boundary of the disk four equations are obtained which 


uniquely determine A, , A. , B, and B,. 





3G. B. Airy, Rep. Brit. Assoc. Adv. Sci., 1862. 
‘C. Bassel Smith, Q. App. Math. 6, 452-456 (1949). 
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It follows that 


A, = Be = @: 
, = —(a + acb)[Saa(l + €B'oy2) + Sheoy(€B + ox) (19) 
eT al€a? + o,,)(1 + €8’0,.) — BIL + €a’a,,)(€B° + oy)’ , 


and 


—— (a + Beb)[ SaB(1 7 €a' Oye) + Sheoy.(ea* + 2y)]__ (20) 
a(ea + a,,)(1 + €Bo,,) — BIL + €a'e,,)(€B + o.,) 

The stress function is now determined and it can be shown that it contains as a 
special case the stress function for an isotropic plate. Let a = 6 and set « = 1. Then 
from (6) it is seen that FE, = EF, . This is a necessary condition but not a sufficient condi- 
tion that the plate in question be an isotropic plate. For the isotropic case, it is now 
sufficient that @ — 0; since by (14) K — 1 and (6) reduces to the biharmonic equation 
for the isotropic case. Parts of the stress function become indeterminate but they may 
be evaluated by successive applications of L’ Hospitals rule. 
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Fic. 1. Variation of the shear stress component X, at points along the boundary of a_ rigid circular 
disk of radius a with center at the origin for a plain-sawn plate of Sitka spruce and for an isotropic 


slate.—Ordinates: Ratio of shear stress component X, at points on the boundary of a rigid circular 
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disk to S the normal stress at infinity. 


The stress function for a plate in a state of plane strain was also determined and it 
was shown that it also contained as a special case the stress function for an isotropic 
plate. It is interesting to note that using a different method Professor I. 8. Sokolnikoff? 
obtained the same stress function for an isotropic plate. 

Consider now a large plain-sawn rectangular plate of Sitka spruce containing a rigid 
elliptic disk at its center. It is possible to apply the results to a finite plate since the 


5]. S. Sokolnikoff, Mathematical theory of elasticity (mimeographed lecture notes, Brown University, 


1941) Chap. V. 
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stress concentration is localized near the rigid disk. The z- and y- axes are parallel and 
perpendicular to the grain, respectively. 

Shear stresses are probably most important in producing failure in a wooden plate. 
The shear stress given by (18) evaluated on the boundary of the rigid disk is 


y ees aeA, sin 9 cosO 
ids ere a’e’b” cos’ 6 + a’ sin’ 6 
(21) 
BeB, sin @ cos 0 


B’eb” cos’ 6 + a’ sin’ 0 


In Figure 1 the curve for Sitka spruce was plotted utilizing formulas (19), (20) and 
(21). The isotropic curve was computed for ¢ = 0.3. 

Throughout the Sitka spruce plate containing a rigid circular disk | Y, | is less than 
0.03. The maximum value of X, is 1.23S at the point (a, 0). 

From these results, it is probable that, for S sufficiently large, failure in the plate will 
occur approximately along the lines y = -+ta sin 773° = +0.976a with the crack be- 


ginning at the edge of the disk. 


BOOK REVIEWS 


The mathematical theory of plasticity. By R. Hill. Oxford at the Clarendon Press, 1950. 
ix + 356 pp. $7.00. 


Time, temperature, Bauschinger, hysteresis, and size effects are all explicitly ruled out and the major 
but not exclusive emphasis is on ideal plasticity. However, a wide range of topics is discussed in this 
interesting and invaluable treatise. Mathematical proofs, experimental evidence and the practical evalua- 
tion of theory are kept in excellent balance from the study of stress-strain relations, variational principles, 
and the questions of uniqueness, to the solving of practical metal forming problems. A few relatively 
simple elastic-plastic solutions are included for prismatic bars, and for thick cylindrical tubes and spherical 
shells. However, most of the chapters deal with plastic-rigid techniques and solutions developed by 
the author and E. H. Lee. An extensive discussion is given of the slip line fields for plane problems. The 
necessity for complete solutions is stated strongly and repeated warning is given against the error of 
thinking in terms of static determinacy only and not considering velocity conditions as well. A number 
of interesting and truly amazing solutions are given in detail for which the configuration remains geo- 
metrically similar as the plastic deformation proceeds. Among the miscellaneous subjects covered are: 
machining, hardness tests, notched bars, normal and oblique necking, earing and anisotropy. Tensor 
notation is used throughout without apology. 

As is entirely proper for a book on the mathematical theory of plasticity, the physics of metals is 
covered only by reference to treatises on the subject. A very brief Appendix on suffix notation, the sum- 
mation convention, and hyperbolic differential equations is considered sufficient for the reader who is 
assumed to be familiar with the elementary theory of elasticity. 

The reviewer regretted to detect an apparent desire on the part of the author to have written a truly 
definitive text at a stage in the development of plasticity when too many basically new facts are being 
discovered. As the author himself is in the forefront of many of these developments he could have done a 
greater service by indicating more clearly the needed fundamentally new directions rather than glossing 
over our present shortcomings and often embarrassing lack of knowledge. He then might not have 
allowed himself to get so worried by Southwell and Allen’s elastic-plastic solution for a notched bar that 
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he unjustifiably loses faith in his rigid plastic solutions. Also, there would be no need to-keep a distinction 
between the true yield function and the plastic potential in the derivation of general stress-strain relations, 
nor to place so much faith in the path independence of the dissipated work. 


D. C. DrucKER 


Response of physical systems. By John D. Trimmer. John Wiley & Sons, Inc., New York 
and Chapman «& Hall, Limited, London, 1950. ix + 268 pp. $5.00. 


This interesting book was written from material used in a course in “Instrumentation” primarily 
for physicists and engineers. The book is interesting because it manages to cover in concise form, on a 
level suitable for seniors, the essentials of system response to sinusoidal and pulse driving. Both mechan- 
ical and electrical systems are dealt with and system stability is treated in sufficient detail to give a good 
idea of the problems. Measuring instruments are discussed as physical systems and measurements as the 
response of such systems together with measurement errors. The Chapter headings are: A Pattern for 
Systems, Physical Systems, First Order Systems, Second Order Systems, Sinusoidal Forcing of Linear 
Systems, Higher Order Systems, Measuring Instruments, Feedback Systems, Parametric Forcing, 
Distributed Systems, Nonlinear Systems. Appendices: The Laplace Transform in the Study of Linear 
Systems, Problems and Review Questions. 
R. TRUELL 


Classical mechanics. By H. C. Corben and Philip Stehle. John Wiley & Sons, Inc., New 
York and Chapman & Hall, Ltd., London, 1950. xvii + 388 pp. $6.50. 


As the authors point out in the Preface “the parts of classical mechanics which are of present-day 
interest to the physicist are not those which were of paramount interest in the nineteenth century.” The 
authors therefore set out to write a modern book which would reflect this change in emphasis in an old 
field. In this reviewer’s opinion, they have succeeded to a remarkable extent in producing such a text. 
The space available for this review does not permit a detailed discussion of the contents, and the following 
list of chapter headings (with some significant section headings added in parentheses) must suffice. 
Kinematics of particles. The laws of motion. Conservative systems with one degree of freedom. Miscel- 
laneous theorems on systems of particles Lagrange’s equations of motion (Charged particle in an electro- 
magnetic field). Applications of Lagrange’s equations (Rutherford scattering. The cylindrical magnetron). 
Linear vector spaces. Small oscillations of conservative systems (Application to molecules. The molecule 
X3Y). Rigid bodies. Hamiltonian theory. Contact transformations. Contact transformations which 
simplify the equations of motion. Poisson brackets (Poisson brackets in quantum mechanics). Infinitesi- 
mal contact transformations. Further development of transformation theory. Miscellaneous generaliza- 
tions and analogies (Analogy with thermodynamics). Introduction to special relativity theory. The motion 


of particles in high energy accelerators. 
W. PRAGER 


Methods of mathematical physics. By Harold Jeffreys and Bertha Swirles Jeffreys. Uni- 
versity Press, Cambridge, 1950. vii + 708 pp. $15.00. 


The first edition of this treatise appeared in 1946 [see this Quarterly 5, 366(1947)]. The fact that a 
second edition should have become necessary so soon clearly indicates that the work has become a 
standard reference for mathematical physicists and applied mathematicians. While the material has been 
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rearranged in some chapters and minor changes have been made throughout the book, the major changes 
are as follows. A section on block matrices has been added, the chapter on multiple integrals has been 
thoroughly revised, relaxation methods have been treated in somewhat greater detail, the presentation of 
the theory of inverse functions has been revised, and the treatment of multipole radiation has been ex- 


tended. Also, some new problems have been added. 
W. PRaGER 


Foundations of aerodynamics. By A. M. Kuethe and J. D. Schetzer. John Wiley & Sons, 
Inc., New York and Chapman & Hall, Ltd., London, 1950. ix + 374 pp. $5.75. 


This is a textbook for a general course in aerodynamics, intended for graduate and senior-under- 
graduate students whose training has been in engineering, mathematics, or physics, but who need not 
previously have studied aerodynamics. There are three major parts of the volume, devoted substantially 
to incompressible-inviscid, compressible-inviscid, and incompressible-viscous flows. It is a major attrac- 
tion of the textbook, however, that real gases are discussed in the opening chapter. Here, in addition, the 
plan of the whole book is laid out for the student and the admirable objective of his study is stated: 
“to provide a background of sound concepts for finding approximate solutions of problems in the flow of 
a compressible, viscous, inhomogeneous gas.” 

The chapters based on incompressible-flow theory include, among other matters, classical two-dimen- 
sional thin-airfoil theory and Prandtl wing theory. The compressible-fluid chapters include channel flow, 
shock and expansion waves, Prandtl-Glauert rule, and linearized supersonic airfoil theory. Except for a 
descriptive section, the last-mentioned is confined to plane flow and the equivalent sweptback infinite 
wings. 

The stress terms in the Navier-Stokes equations are derived in an Appendix, in lucid style. Two 
chapters concern mostly laminar boundary-layer flow; two more are devoted to turbulence, and another 
to boundary-layer transition, including descriptions of boundary-layer instabilities of two types. These 
chapters seem particularly valuable, and are probably unique in the field of fluid-mechanics textbooks. 

The closing chapters present material on viscous-compressible flow and some descriptive information 
on wing properties, such as transonic flow and stalling. The viscous-compressible paragraphs are soundly 
physical, concerning especially energy balance and aerodynamic heating. There is also some information 
on interactions between shock wave and boundary layer, and a regrettably short paragraph on “‘the thick- 
ness of a shock wave.” 

There are problems for the first 13 of the 18 chapters, which seem, to the reviewer, to demand rather 
little originality from the student. There are gratifying exceptions to this generalization, however. 

Throughout, the authors have emphasized the physical picture that their equations describe. More- 
over, the reader is frequently reminded of the confirmation, and occasional conflict, of experimental 
results. Thus, although it will undoubtedly teach most seniors and young graduate students a good deal 
of applied mathematics, it is really not a book in that field, but one in the field of physics, as it should be. 

In view of the purpose and the scope of the volume, it is probably not justifiable to complain about 
sins of omission; rather, the authors should be commended on their skill in having selected their topics 
well. It may suffice to remark that no single volume of 374 pages can be adequate for full-dress presenta- 
tions of theoretical aerodynamics at the graduate level. In its own field—which is one that should be 
vigorously encouraged—this is certainly the best textbook, and a badly-needed one. 


W. R. Sears 


Theory of probability. By M. E. Munroe. First Edition. McGraw-Hill Book Company , 
New York, Toronto, London, 1951. viii + 211 pp. $4.50. 


This is an elementary textbook designed to introduce students who have had a first year course in 
calculus to the theory of probability. The author makes frequent references to Cramér’s and Uspensky’s 
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books, and the general purpose of the book seems to be to provide a book more accessible to less advanced 
students than these books are. As a matter of fact, the present book is more in the spirit of probability 
theory than is Uspensky’s book. Even a version of the strong law of large numbers is discussed. However, 
strong convergence is defined as convergence in the space of infinite sequences. This is rather sophisticated 
for students with a year of calculus and an alternate limit definition might have been more intelligible. 
The book has many illustrative problems and exercises for the student. 


J. Wo.Lrow1tTz 


Vector and tensor analysis. By Harry Lass. McGraw-Hill Book Company, Inc., New 
York, Toronto and London, 1950. xi + 347 pp. $4.50. 


This book contains a thorough treatment of vector analysis, with extensive applications. It also 
contains a short but good introduction to tensors. There is a large number of examples and problems. 
The style is clear and concise, and there is a sharp segregation between general theory and the many 
applications. This renders the book suitable for the standard undergraduate course in vector analysis, as 
well as for more advanced work dealing with an introducton to mathematical physics. 

There are nine chapters. Chapters I, II and IV deal with the algebra, differential calculus and integral 
calculus of vectors. These chapters, comprising about one quarter of the book, contain the theory appear- 
ing in the standard undergraduate course in vector analysis. Chapters III, V, VI and VII, which comprise 
about one half of the book, deal respectively with differential geometry, electricity, mechanics, and 
hydrodynamics and elasticity. Of course, in a book of this nature the treatments of these topics can be an 
introduction only, but the treatment is systematic, and particularly in the case of electricity is rather ex- 
tensive. Chapters VIII and IX, which comprise the remaining one quarter of the book, contain a brief 
introduction to the theory of tensors, as well as brief applications of tensors to geometry, mechanics and 


elasticity. 


G. E. Hay 


Elasticity. Proceedings of the Third Symposium in Applied Mathematics of the American 
Mathematical Society, Volume III. McGraw-Hill Book Company, Inc., New York, 
Toronto, London, 1950. v + 233 pp. $6.00. 


This volume contains the papers presented at the Third Symposium on Applied Mathematics of 
the American Mathematical Society, held at the University of Michigan on June 14-16, 1949. Seventeen 
papers are included, two of these being presented in the form of abstracts. 

The subject matter of thirteen of the papers is elasticity, that of the remainder, plasticity. The 
papers in elasticity include studies in the plane problems of anisotropic materials, finite deflections of 
axisymmetric shells and asymptotic integration in shell theory, the bending of plates, general theory of 
finite elastic deformation, the beam of circular cross-section under concentrated loading parallel to its 
axis, approximate determination of the torsional rigidity of beams, torsion of the axially symmetric 
shaft, and the boundary layer edge effect in elastic plates. The plasticity contributions are on the stress 
analysis of indeterminate elastic-plastic structures, the use of characteristics in obtaining graphic solu- 
tions of plane problems, and the solution of problems in plane flow with stress discontinuties. 


WiuuraM H. Pei 
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SUGGESTIONS CONCERNING THE PREPARATION OF 
MANUSCRIPTS FOR THE QUARTERLY OF 
APPLIED MATHEMATICS 
The editors will appreciate the authors’ cooperation in taking note of the following directions for the preparation of 


manuscripts. These directions have been drawn up with a view toward eliminating unnecessary correspondence, 
avoiding the return of papers for changes, and reducing the charges made for “‘author’s corrections.” 





Manuscripts: Papers should be submitted in original typewriting on one side only of white paper sheets and be 
double or triple spaced with wide margins. Marginal instructions to the printer should be written in pencil to 


distinguish them clearly from the body of the text. 
The papers should be submitted in final form. Only typographical errors may be corrected in proofs; composi- 


tion charges for all major deviations from the manuscript will be passed on to the author. 

Titles: The title should be brief but express adequately the subject of the paper. The name and initials of the 
author should be written as he prefers; all titles and degrees or honors will be omitted. The name of the organization 
with which the author is associated should be given in a separate line to follow his name. 

Mathematical Work: As far as possible, formulas should be pgeinen: Greek letters and other symbols not 
available on the typewriter should be carefully inserted in ink. Manuscripts containing pencilled material other 


than marginal instructions to the printer will not be accepted. F 
The difference between capital and lower-case letters should be clearly shown; care should be taken to avoid 


confusion between zero (0) and the letter O, between the numeral one (1), the letter / and the prime (’), between 


alpha and a, kappa and k, mu amd u, nu and », eta and n. 
The level of subscripts, exponents, subscripts to subscripts and exponents in exponents should be clearly 


indicated. 
Dots, bars, and other markings to be set above letters should be strictly avoided because they require costly 
hand-composition; in their stead markings (such as primes or indices) which follow the letter should be used. 


Square roots should be written with the exponent } rather than with the sign +/. 
Complicated exponents and subscripts should be avoided. Any complicated expression that recurs fre- 


quently should be represented by a special symbol. i a 
For exponentials with lengthy or complicated exponents the symbol exp should be used, particularly if such 


exponentials appear in the body of the text. Thus, 
exp [(a? + b*)"/*] is preferable to e***+**)*/* 
Fractions in the body of the text and fractions occurring in the numerators or denominators of fractions should 
be written with the solidus. Thus, 





ez 

cos — 

cos (xx/2b) . 2b 
—————— is preferable to 

cos (ra/2b) 7a 

cos — 

2b 


In many instances the use of negative exponents permits saving of space. Thus, 
, , sin u 
J ué sin u du is preferable to / —— du. 
u 


Whereas the intended grouping of symbols in handwritten formulas can be made clear by slight variations in 
spacing, this procedure is not acceptable in printed formulas, To avoid misunderstanding, the order of symbols 


should therefore be carefully considered. Thus, 
(a + bz) cos t is preferable to cos t (a + bz). 


In handwritten formulas the size of parentheses, brackets and braces can vary more widely than in print. 
Particular attention should therefore be paid to the proper use of parentheses, brackets and braces. Thus, 


{{a + (6 + cxr)"] cos ky}? is preferable to ((a + (b + cx)") cos ky)’. 


Cuts: Drawings should be made with black India ink on white paper or tracing cloth. It is recommended to submit 
drawings of at least double the desired size of the cut. The width of the lines of such drawings and the size of the 
lettering must allow for the necessary reduction. Drawings which are unsuitable for reproduction will be returned 
to the author for redrawing. Legends accompanying the drawings should be written on a separate sheet. 
Bibliography: References should be grouped together in a Bibliography at the end of the manuscript. References 
to the Bibliography should be made by numerals between square brackets. 

The following examples show the desired arrangements: (for books—S. Timoshenko, Strength of materials, 
vol. 2, Macmillan and Co., London, 1931, p. 237; for periodicals—Lord Rayleigh, On the flow of viscous liquids, 
especially in three dimensions, Phil. Mag. (5) 36, 354-372 (1893). Note that the number of the series is not separated 
by commas from the name of the periodical or the number of the volume. 

Authors’ initials should precede their names rather than follow it. 

In quoted titles of books or papers, capital letters should be used only where the language requires this. Thus, 
On the flow of viscous fluids is preferable to On the Flow of Viscous Fluids, but the corresponding German title would 
have to be rendered as Uber die Strémung zdher Flissigkeiten 

Titles of books or papers should be quoted in the original language (with an English translation added in 

parentheses, if this seems desirable), but only English abbreviations should be used for bibliographical details like 
ed., vol., no., chap., p 
Footnotes: As far as possible, footnotes should be avoided. Footnotes containing mathematica] formulas are not 
acceptable. 
Abbreviations: Much space can be saved by the use of standard abbreviations like Eq., Eqs., Fig., Sec., Art., etc. 
These should be used, however, only if they are followed by a reference number. Thus, “Eq. (25)”’ is acceptable 
but not ‘“‘the preceding Eq.” Moreover, if any one of these terms occurs as the first word of a sentence, it should 
be spelled out. 


Special abbreviations should be avoided. Thus “boundary conditions” should always be spelled out and not 
be abbreviated as “‘b.c.,” even if this special abbreviation is defined somewhere in the text. 





The Eighth International Congress for Applied Mechanics 


will be held at the University of Istanbul in Istanbul, Turkey, from August 20 to 28, 1952. 

Membership. Membership in the Congress will be open to all qualified persons, whether they are 
able to be present in person or not. For regular members of the Congress the fee is $10: these persons will 
be entitled to participate in all the scientific and social features of the Congress and may, if they wish, 
contribute papers. Members of families, accompanying Congress members and not participating in the 
scientific meetings, may become associate members for a fee of $5: they will be entitled to all the privi- 
leges of membership. 

Contribution of papers. Each member of the Congress may present at most two papers, and the 
time allotted for each paper will not exceed fifteen minutes. Facilities for projection of formulae and 
graphs will be available in order to save time; members wishing to make use of these should prepare 
slides beforehand not exceeding 8 x 8 em (3 x 3’). Abstracts for contributed papers should not exceed 
400 words in length and must be submitted on blanks which may be secured from the Secretary of the 
Congress. Abstracts must be in the hands of the Organizing Committee not later than June Ist, 1952. 
The Organizing Committee will endeavour to accept papers arriving after this date but will be under no 
obligation to do so. Members of the Congress who will not be able to attend in person may send only one 
paper which must be read at the Congress by a member of the Congress designated by the author. The 
official languages of the Congress are English, French, German, and Italian. 

The Sections for the presentation of short contributed papers will be as follows: 


SECTION I: Elasticity—Plasticity—Rheology 

SECTION 11: Fluid Mechanics (Aerodynamics—Hydrodynamics) 

SECTION 111: Mechanics of Solids (Ballistics—Vibrations—Friction—Lubrication) 
SECTION Iv: Statistical Mechanics—Thermodynamics—Heat Transfer 

SECTION v: Mathematics of Physics and Mechanics—Methods of Computation 


The Proceedings of the Congress will contain in extenso the speeches, lectures, major addresses and 
short contributed papers read at the Congress. If any of these are not available, the abstracts will be 
printed in their place. Members who wish to receive the Proceedings of the Congress will be required to 
pay a fee which will be announced at a later date. 

Five minutes will be allowed for the free discussion of the contents of each paper read at the Congress. 
The discussions will not be printed in the Proceedings. 

Major addresses and lectures. The Organizing Committee intends to invite a few outstanding 
persons to deliver stated addresses. There will also be a lecture by a specialist in each section of the 
Congress, mainly in order to present research in fields in which important advances have been made in 
the recent years. 

Entertainment. ‘There will be many interesting entertainment features, amongst them a reception 
and a banquet which al] members of the Congress are entitled to attend. Excursions by boat or bus are 
also being considered. Finally, members will be given the opportunity of visiting the main historical 
monuments, and architectural and natural beauties of Istanbul, under the guidance of competent persons. 
In connection with the Congress, Messrs. Henry Van der Zee and Co. of New York and Istanbul, acting 
as agents for the American Express Co. Inc. of New York, will provide for organized tours to Troy, 
Bergama, Ephesus, Bursa, Izmir, and Ankara. 

Travel. The Organizing Committee has designated as the official travel agency of the Congress the 
American Express Company, Inc. of New York, from which all information concerning travel and 
accommodation can be obtained. 

The Organizing Committee is providing for a 25% reduction on the one way fare and a 50% reduction 
on the round trip fare on all boats belonging to the Turkish State Lines, operating in the Mediterranean, 
and on trains of the Turkish State Railways. It is also expected that transportation by sea from the 
United States will be made available for approximately $450 for the round trip. 

Accommodation. First class hotel rooms are available for $3.50 upwards, whilst second class rooms 
are from $2 to $3.50. Free accommodation will be provided in school dormitories. 

Meals. Excellent meals in first class restaurants cost from $1.50 upwards. Also meals will be pro- 
vided by the University cafeteria for approximately $0.65 upwards. 

Information. All requests for information concerning the Congress should be addressed to: 

The Eighth International Congress for Applied Mechanics, P. O. Box 245, Istanbul (Turkey) 


Chairman of the Organizing Committee: Secretariat of the Congress: 
Cahit Arf (General Secretary), 
Giacomo Saban (Acting Secretary). 


Kerim Erim. 























